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Introduction. It is the purpose of this work to derive functions which represent the 
flow of a perfect fluid inside a cylindrical wall. The velocity perpendicular to this wall 
must be zero everywhere. Some of the problems discussed will have the axis of the 
cylinder as an axis of symmetry, but several problems are discussed where the flow is 
not symmetric about an axis. 

In Sec. 1 a certain set of boundary conditions is shown to lead to a series involving 
Bessel functions. This solution is shown in Sec. 2 to be the result of a point sink at the 
origin and a circular ring source lying in a plane perpendicular to the axis of the enclosing 
cylinder and through the origin. This discussion enables the stream function for a source 
alone and a ring alone to be found. The problem of the combined ring source and point 
sink in a uniform stream through the cylinder is next discussed. The stream function 
for a point source in the cylinder not on its axis is then found in Sec. 3.2. 

In Sec. 3.4 a new result is obtained which shows that the potential due to a sink not 
on the axis can be decomposed into a sum of terms, each term being the potential due 
to a ring source circular in shape but with non-constant strength. It would seem that in 
many ways these particular ring sources offer the simplest possible configurations with 
which to begin the study of fields which do not possess axial symmetry. 

In See. 4 the field of a source and sink on the axis in a stream inside the cylinder is 
discussed and the semi-axes of the resulting closed stream-surface are computed. 

In Sec. 5 the field of the ring with strength M cos @ is worked out in detail, this being 
the least complicated of the ring sources of Sec. 3.4. 

1. Let ¥(p, 2) be a stream function in cylindrical coordinates for a field which is 
symmetric about the z-axis. Then ¥(p, z) is a solution of the equation 


ay _ lay, av _ 4 
Op pop 27 


which has solutions of the type 


¥(p,2) = >» An exp (—jim2/@) pI s(jime/Q) , 


where j,,, is the mth positive zero of the Bessel function of the first order J,(x). It follows 


*Received February 8, 1960; revised manuscript received June 23, 1960. This work was supported 
by the Ordnance Research Laboratory, Garfield Thomas Water Tunnel, under U. S. Navy Contract 
NOrd 16597. 
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that ¥(a, z) = 0 for all z since J,(j,,,) = 0, and ¥(0, z) = 0 for all z since p is a factor in 


each term. There is an expansion theorem that if f(t) be a function arbitrarily defined 


in the interval (0, 1) and 
al 


l t/*4(t) dt 


exists, and if 
al 


a, = {2/[JAj,.)]°} Uf(Od sGimt) at; 


e 


then 
2: Bad (jrat) 
m=~1 
converges and its sum is (1/2) {f(z + 0) + f(a — 0)} at any internal point z of an 
interval (a, b) such that 0 < a < b < 1 and f(¢) has limited total fluctuation in (a, b).* 
Consider the function F(p) 
ae \M, 0 < p< ka, & <I 
F(p) =< 
0, ha <p = 4, 
then F(p)/p has limited total fluctuation in any interval 0 < « < b < 1, and 


1 k k 


[ ¢?RO/d at = [ Me? at = 2Mt”?| = 2Mk” 


@ 


exists (where t = p/a). 
Therefore, the theorem applies to this function and 


ll 


al 
2/laJiiim)] | F(at)J \(iml) at, 


~0 


Am 


Il 


2M/{aJi(irm)] | Ji(jimt) dt, 


= 2M[1 — Joljimk)V/Girmd iim); 


since 


a 


| J(u) du = —ZJ,(u). 
Now with A,, = a,, the solution for Y becomes 
¥(p,2z) = p = {2M[1 — Jo(jimk))/[QjimI 2(j1m))}JFi(jime/@) exp [(—Jjim | 2 |)/a] (1.1) 
and 


. 2M[1 — Joljwk)] 7: 
+ 
ed Aji md 2(jim) JsGime/a) 


¥(p, 0) 


;' aa (M, O<p<ka 
= pF(p)/p = F(p) = 4” r ) 


| 


0, ka<p<1 


*G. N. Watson, Bessel functions, Cambridge Univ. Press, 1922, Sec. 18.24. 
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It should be noted that J2(jim) = J3(jim) So that only values of J, and J, are needed. 
Hence we have a flow such that y = 0 consists of the cylinder p = a and the part of the 
plane z = 0 between p = ka and p = a, and the z axis. ¥ = M consists of the part of 
» = 0 inside p = ka. The plot of the values of the function shows that the origin is a 


sink into which the fluid flows and there is a ring source of radius ka in the plane z = 0. 


(Fig. 1* was made for k = 1/3, M = z.) 


1.0)P 











a a ee a ca. 


“10-50 O .50 1.0 z 


Fia. 1 
2. If we take 


oe « 
"7 - _) p + V(p, 2), 


Ld ° ° . . . . 
where ¥(p, 2) is given by Eq. (1.1), we will have a stream flowing around the ring source 
. . 9. . . 
and simple sink at its center. If q, , g, are the velocity components, we have 


= Loy _ - 7 a/ (; e) ium 9 
q, = —™ >» A,, exp (—jimz/a)J o jim 7) a (2.1) 
_lay¥__j5 > (—j..2/q) im (i 2) a 
q, = a p> A,, exp (—jimz/a) - J, jim a) (2.2) 
If p = 0, J,(0) = 1, so on the z axis 


. = Ache ‘ 
qg(0,z) = U — > —— exp (—jimz/a). 


m=1 


The stagnation point on the axis will be given by the root of the equation 


aU >> 1 — Jol(jimk) ; 
—— = —, *~—— exp (—j1.2/a) (2.3) 
2M m=1 J (Jim) 
when solved for u and z = ua. It will be noted that if a stagnation point can occur at a 
place not on the z-axis g, and g, both are zero. So from (2.2) a value can be assigned for 
z and the resulting equation for p can be solved, giving the coordinates of the stagnation 
point. Then these values can be inserted in (2.1) and the value of U/M can be calcu- 
*Caleulations were done by the Ordnance Research Laboratory computers. See Appendix I for 


values of y, and Fig. 1. 
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lated, giving the velocity of the stream at a great distance from the origin. It was found 
that z = 0.20 gives an equation 
f(p) = 1.0547, (3.8817) + 2.6207, (7.1056 ) 

+ 2.825J, (10.1725p) + 1.9487, (13.3237) 

+ 0.972J, (16.4706p) + 0.4547, (19.6159p) 

+ 0.278J, (22.7601 p) + 0.226J/, (25.9037 p) 

+ 0.167J, (29.0468) + 0.1057, (32.1897 p) 

+ 0.0587, (35.3323) + --- = 0. 


This gives the following table of values: 
a) 


ry 0 0.3 0.4 0.49 0.5 O07 

f(p 0 2.450 0.668 0.043 —(0.020 —0.074 
therefore, z = 0.200, p = 0.497 is a stagnation point. This is approximately z 0.20, 
p = 0.50. The value of U/2M = 0.649 from (2.2); or U = 1.298m. Equation (2.3) gives 
a stagnation point at — 0.65 on the z axis. By direct calculation from (2.0) the values 


of ¥/M can be calculated.* 


PIO te ee 




















3. Since the problem of the sink in a cylinder where the flow is along the cylindrical 
wall is a Neumann problem, it is necessary for @/dn to vanish over the cylinder. The 
velocity potential satisfies the differential equation: 


a® _ 1 0o® 1 1 o® 4 eo 0. 


Op hall pop p dg Oz 


@=e ” cossgP(p), 


ele 2,7 4 eP = 6, 

dp pdp = p 

, dF +p dl + (k p” —s)P = 0. 
dp 


P dp 


Therefore, P = J,(kp). 


*See Appendix II for a table of values and Fig. 2. 
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&’= >) Dd Axe cos spJ (kp) 


s=0 


08/dp = > > A,.e"*"(cos sg)kJ'(kp). 


s=0 


/ 


Therefore, if d&/dp = 0 for p = a, J’(ka) = 0. If h,.,, is mth root of Jf = 0,k = h, ,,/a. 
3.1 For any axially symmetric field s = 0, J{(ka) = — J,(ka) = 0; so 


@= > A, exp (—jin | 2 |/a)Jo(j:.0/a), 
n=1 


08/dz = — > A,(jin/a) exp (—jin | 2 |/a) J o(jinp/a), 
n=1 
O® [ 0; p =0, simple pole] 
ee | po <i ae, 
0; O0O<p<a 


| | (2) pJ(j;,p/a) d0dp = —2xM if M 
e~) v¥W Ox 0 


is the strength of the sink, because 


(2%) = —M io). 


Oz 


Hence, the integral = .WJ,(0) {$fo" d 6 dp = total flux from the sink on one side of 
z=0 = —-2rA, [ (j:,p/a) J o(j:-p/a)J o(j,,0/a) dp, 


since every integral vanishes except for 


n=r= —2rA,,y,,aJ6(j:,)*, 


A, = 2M/[j,,aJo(n,)); 


hence, 


@= D> (2M/[j,,aJi(jin)]} exp (din | 2 |/@So(jine/a) 
n=l 

is the velocity potential of some field with a sink at the origin with normal velocity 

O on p a. 

The stream function is 


Y= DY (2M /LiaI2isn)]} oJ (irne/a) exp (— inn | 2 |/a), (3.1) 


n=l 


*G. N. Watson, Bessel functions, Cambridge Univ. Press, 1922, Sec. 5.11 (11). 
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which obviously gives ¥(a, z) = 0 (since J,(j:,) = 0). Also note that this is a part of 
Eq. (1.1). However, this ecuation gives a velocity 0 at z = ©, so it is necessary to add a 
stream flowing from right to left; i.e., a term Mp*/a*, then the resulting field is that of 
a sink in a cylinder. 

If a stream of velocity U from right to left is added (Fig. 3): 


hill — ite iedee Mp U 
wy = (| z |/z)(2Mp/a) = {jind ond} J 1CGinp/@) exp (—jin | 2 |/a) + £ +—p. 
ae Pare, a 2 


EAA SALA PAE ASAE EAL E LES 


oe 

+ 

e— 

*+— VELOCITY 
ee rae 


Se y 








~~“ 


_ EEO 
-_— 
—— 


WVSESSVAVANUSANS ANN NANS YS 


> 


Fic. 3 








Assume a = 1. Then for z > 0 


v = 2Mp D> finJojun)} J i(jinoy EXP (— jin | 2 |) + (a1 + DN 


for z < 0, 


eo , —_— U 
Y = —2Mp D> {ism aiin)} J (ino) EXP (—in | 2 |) + (J - M)o; 


l dy O® = r . 5 27: . 4 | ‘ r 
a’ teak “ae > {[2MJojino)]/ Join) } exp (—din | 2 |p + 2M — U, 
U>2M 
on 
p=0, q7 =O if D Jilin) exp (-is ° ) —~ OM I. 
Since 
Js(jin) exp (—i., z |) jin+t . 1; 
_> . lee ak _— 1.n+ l,n Z 
Joli nti) CXP (—)1., 2 |) jin exp J parry ) 


= (1 +2) exp (—za | z |) & exp (—7 | z }) 


as n — o, the series converges if | z| > 0. 


The graph of U/M as a function of z is made from the following Table which was 
computed from the formula just derived. 


Li) 


2 l ] 1/3 


U/M | 0.150 | 1.376 3.769 
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The graph of U/M = 2 + 2 >> can be compared with the function (U/M)'” which 
gives the position of the stagnation point for a source in a free stream and it is evident 
that a stronger stream is needed inside the cylinder to push the stagnation point to the 
position it has in a free stream. The curves are plotted in Fig. 4. 








The 
Of 
Uv 7} 
M L 
5+ 
3 CYLINDER 
If FREE 
1 ly | 
4% z 
Fia. 4 


3.2. In the case where the sink of strength VW is at (b, 0, 0) the same method can 
be used. Let 


po = y > ye exp (—Ren lz | a)d (Rene a) COS Spy, (3.2) 


e=0 n=1 


where h,,, is the nth positive root of J/(z) = 0. Now 





(8@/02z),-0 = —(h,,/a)P(p, ¢, 0). 
90° 60° 
Z=0.10 / 
©MAX@=.6321 
AT P = 1142 / 
® <05615 





NUbaoON ®wO 


>. 








o° 
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Hence 


[ i’ (22) cos sgJ ,(h.np/a) p dp dp 
0 v0 « s=0 


= (1 ¥ 5s0)4(hsn/) Asn: (a’, 2)[h.n —_ 5h pee po | 
where 





0, a> st 
90° 6 
Z=.15 
© MAX @ =.3472 
AT P =.1378 





NWO N® © - 











But (06/0z),-» is zero everywhere except at the point (b, 0, 0), so the integral must be 
equal to 
J (h,,/a)(—27M). 
The integral must be interpreted as a Stieltjes Integral of a discontinuous function. 
Therefore 
—2rM J (h,,b/2) = —7A,,a(hi, — 8 (hen) J i(hen(L + 8,0). 
Therefore 
2Mh,A2 — 840) 


. ~~ 2 2») 72 t ( ij : 
= (iin — FF Uhan) on) 


3.3. In the special case where we demand symmetry with respect to the axis we get 
a set of terms 


6 = oma > FinlTo(jin)] exp (—jin | 2 [/a) So(jind/a) J o(jine a} 4+ 2Mz 











eee 
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90° 60° 
Z=0.20 


© MAX ® =.2152 
AT P =.1675 


















0° 


which still satisfy the differential equation and from (1.1) and (2.1) this must be the 
potential of a ring sink with 2M for its strength and b for its radius, inside a cylinder of 
radius a after the addition of the stream term. The stream function for this will be 


Y= 2Ma*'S D7 jralJoiw]* exp (— jnt/a)FoGrnb/a)J (jrno/a)} + Mame. 


3.4. A ring source with variable strength [//(2rb)] cos s’¢) at the point (b, go , 0) 
gives for the potential of the whole ring 


ale oo 


©, = M(2r)~' cos a> b A,,, exp (—h,, |z|/a)J,(h.,p/a) cos s(y — o) | dgo 5 
n=1 


/0 \s=0 


COS 8’¢o COS 8(9 — go) dgo = 


[" [0; s’ #8, 


le COS Sy; s’ = 8; 


ee) 


&, = (m/2) cos p> A,, exp (—hin | 2 |/a) J (Rene, a}. 
n=1 
Since these are exactly the terms of the original series Eq. (3.2.1), which contain cos 
sy, we have the theorem that the potential @ for a point source at (b, 0, 0) of strength 
M on a cylinder consists of the sum of the potential ®, for s = 0, 1, 2, --- where ®, is the 
potential of a ring source of variable strength in the cylinder and passes through the 
point source and has a strength [1//(27b)] cos sg at the point (b, go , 0) of this ring. 
The ring for s = 1 will be referred to as a ring of the first order. The ring of zero 
order is just a ring of constant strength since cos sg. = 1 if s = 0. 
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In Sec. 5 below, the field of the first order ring will be worked out in detail. 

4. Let there be a source of strength 7 at p = 0, 2 = b and a sink of strength M at 
p = 0,2 = —b and a stream flowing from right to left inside a cylindrical wall with 
axis along the z axis. If z > b the streaming motions at infinity due to source and sink 
cancel each other and these terms may be omitted. The same thing occurs for z < — b. 
But between the two singularities — b < z < b there will be a streaming motion from 


source to sink. 
Forz> b 


yp = (U/2)p? — doe2M pj iA[Jo(jin)] J s(jine) {ExP (—jin | 2 — b |) — exp(—ji, [2 +5 |)} 


1 
= (U/2)p? — Dd 2M pj) [Jo(iin)] Ji Ginp)(Sinh j,,b) exp (—jir2)- 
n=1 
If y = 0, we will have the dividing stream surface so this surface has the equation 


4 U/2 — >» 2M [Joljin)] + [271Gine)/Cinp) (sinh j,,b) exp (—ji.2)7 = 0. 


A 
If p + 0, the part of this stream surface, not coincident with the axis p = 0, has for its 
equation 


U/(4M) = YS [Jolsine) + Join) ][JoGine)] “(sinh j,,b) exp (—jin2) 


n=1 


since (2/2) J,(z) = Jo(z) + J2(z). Now the points on the axis p = 0 belonging to this 


surface are the stagnation points, so 


U/(4M) = >> [Jolin] (sinh j,,.b) exp (—j,,2) 

















90° 60° 
| 2 +0.50 FY 
Fs | 
© MAX ® = .04584 / | 
AT P = .4315 F, 
fo | 
04024 Fé | 
Loe: 245 | 
2.035 | 
P=.030 | 
30° 
aN 
e” 
\ \O 
IT; @-° 
Y, J }o° 
— ) 
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1961] 


is an equation for the z coordinate of the stagnation point. This gives the value U/M = 
0.2828 for z = 1.50, and U/M = 2.6120 for z = 1.00. 
For-—b<z<b, 


y = (U/2)p° + > 2M pjinlJi(jin)) “J o(jinp) [exp (—jinb)]2 cosh j,,2 + 2Mp’. 


Now ¥ = 2M, z = 0 gives the position where the boundary of the closed dividing stream 
surface crosses the plane z = 0. Since p ¥ O the equation 


2M + (U/2) + 2) 2M[Jolirw)]*[Joliinp) + Ja(jinp)] exp (—jind) = 2M p™ 


gives the intercept. This gives p = 0.955 for the first case and p = 0.742 for the second. 


90° 60° 


| 
2 =1.0 / 


] = .018585 










4 


NUDAaAON®DwWO 


e™ 











5. In the case of the most simple ring with variable strength discussed in (3.4), we 


have 


S = 3 cosy D> 2Mhigd (hinb/a)(h2, — 1)7"[Jx(Ian)] 2 Gin /a) exp (—Iin | z |/a)- 


n=] 


To obtain some numerical values we will take b = a/10, a = 1; so that the ring is one 
tenth of the radius of the cylindrical wall and in a plane perpendicular to the axis of 
the cylinder with its center on the axis. It is necessary to have the values h,, which are 


zeros of J{(z); since 2J{(z) = zJ9(z) — J,(z) and the value z = 0 is not one in which we 
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Fic. 10 


are interested; these values can be found much more readily from the second form, i.e., 


zJ,(z) = 0. The values found are as follows 
n l 2 3 4 5 6 
ee Neaneeane ae Sem - a aa 
hin | 1.8412 5.3314 | 8.5363 11.7060 14.8636 18.0155 
n 7 S | 9 10 ll 12 
hin 21.1644 24.3113 | 27.4571 30.6019 33.7462 | 36.8900 
| 








In computing @ for z = 0.10, it was necessary to use 12 terms of the series to get con- 
sistent values, a smaller number sufficing for larger values of z. 
It will be noted that the strength of the source along the ring varies from a maximum 


M at ¢ = 0, to 0 at gy = 2/2, toa minimum — M at gy = zo, and back through 0 at 
g = 3nr/2 to M at o = 22x. Evidently 6 = + © on the upper half of the 
ring and 6 = — - on the lower half; and any surface @ = c must pass through the points 


(1/10, x/2, 0) and (1/10, 37/2, 0). If c is large and positive, we must obtain a sausage- 
shaped surface enclosing the upper half of the ring; and if c is large and negative, one 
enclosing the lower half of the ring. The diametral plane g = 1/2, 9 = 37/2 is the surface 
® = 0. For small values of @ there is a closed curve on the cylinder p = a which is symmet- 
ric to the plane z = 0 and to the plane & = 0, in which the equipotential surface meets the 
cylinder orthogonally; and then passes down through the two points of strength 0 on 
the ring without cutting the plane surface = 0 at any other points. In calculating the 
values of ® the values z = 1.0, 0.5, 0.2, 0.15, 0.1 were used, and values of p from 0 to 
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1 for about 20 values, so chosen as to give the more interesting situations. Graphs of the 
equipotentials are given in Figs. 5-9 for each section perpendicular to the axis at the 
constant values of z in the above series. A sketch of the equipotentials is given in Fig. 
10. The streamlines are the curves which are orthogonal to the equipotential surfaces. 
Since the ring is a curved line-source, the streamlines from any point on the ring must 
all leave the ring in the plane normal to the ring at that point, and these streamlines 
form a closed surface which meets the ring again in the point symmetric to the starting 
point below the surface = 0. These closed surfaces are nested and close down to point 
limits at the two points of zero strength. They also have as an outer limit, the surface 
of the bounding cylindrical wall. The streamlines which leave a point on the ring source 


are shown in lig. ‘. 


Acknowledgment. The author wishes to express his thanks to Dr. George F. 
Wislicenus for directing his attention to the problem and for many helpful discussions 
as the work progressed. The numerical work would have been almost impossible without 
the indefatigable aid of the Ordnance Research Laboratory’s computing group, par- 
ticularly Mr. Paul C. Sweetland and Mrs. Barbara R. Kautz. For the figure of the 
equipotential surfaces of the first order ring, which realizes the form of this solution so 
completely, the author is deeply indebted to Mrs. Joan P. Lampman. 


Appendix I 


The stream function 


Wp, 2) = 2ep D> [exp (—juz)][1 — Foljn/3) ix (Jolin) 2d Gn) 
m=1 


must be evaluated for 0 < z,0 < p < 1, where J,(z) and J,(z) are the Bessel functions 
of order zero and one, and j,, is the mth positive zero of J,(z). Intervals of 0.1 were chosen 
for p and the series was evaluated for z = .25 to z = 2.00 with intervals of 0.25. This table 
and that in App. II were computed by the Computing Section of the Ordnance Research 


Laboratory. 











Fia. 11 
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TABLE 1—y 
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z/p a 2 3 A 5 .6 ia 8 9 
0 1 1 1 0 0 0 0 0 0 
.25 . 1687 . 4502 . 5639 .4847 .3437 . 2262 .1445 .0867 0415 
50 .0245 .0814 . 1360 . 1646 . 1637 . 1426 1115 .0763 .0393 
40 .0059 0212 .0396 .0550 .0631 .0628 .0550 .0408 0221 

1.00 .0018 0067 .0133 .0195 .0239 .0253 .0255 .O181 0101 

1.25 .0006 0024 .0047 .0072 0091 .0099 .0094 .0074 0042 

1.50 .0002 0009 .0018 .0027 .0035 .0038 .0036 .0029 .0016 

1.75 0001 0003 .0007 .0010 .0013 .0015 .0014 .0012 0006 

2.00 0001 .0003 .0004 .0005 .0006 .0005 0004 0002 

Appendix II 
Values of ¥/M Computed from Eq. (2.0) 

p/z — .30 20 —.15 Q— O+ 10 15 . 20 30 
0 0.000 000 0.000 —1.000 1.000 0.000 0.000 0.000 0.000 
] 036 —0.994 1.006 0.147 0.049 
2 067 —0.052 —0.975 1.025 0.187 0.004 0.114 
Oo —0.102 —0.941 1.059 0.219 
} 0.171 020 0.000 +0.105 0.105 0.237 0.208 0.188 0.091 
5 0.182 104 0.082 0.162 0.162 0.21 0.207 0.220 0.182 
6 0.283 198 0.204 0.234 0.234 0.259 0.263 0.283 
7 0.260 0.318 0.318 0.264 0.243 
8 0.415 0.415 
9 0.525 0.525 

1.0 0.649 0.649 0.649 0.649 0.649 0.649 0.649 0.649 0.649 

A few other values are of interest: For p = 0—, Z = 1/3-—, ¥ ~ 1.071, and for p = 04+, Z = 1/34, 
y = 0.071, y has a period 2 and lines for which y = a or a + 2 etc. are the same line, as illustrated by 
the line through the stagnation point y = 0.22 ory = —1.78. 


Appendix III 


Figures 5-9 are cross sections of the field of the first order ring at distances z = 0.1, 
0.15, 0.2, 0.5, 1.0 from the plane of the ring, the curves being the equipotentials. These 
were made from computations done by the ORL Computing Section. 

Figure 10 is a sketch of the general form of the equipotentials in one octant of the 
cylinder. This is the work of Mrs. Joan Lampman. 

Figure 11 is a sketch of the stream lines which emanate from one point of the ring 
source in the case of the first order ring. 
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ON THE TIME DERIVATIVE OF TENSORS IN MECHANICS OF CONTINUA* 


By 
P. M. NAGHDI (University of California, Berkeley, Calif.) 
AND 
W. L. WAINWRIGHT (U.S. Naval Postgraduate School, Monterey, Calif.) 


1. Introduction. In describing the motion of a continuum, not only have spatial 
coordinates been employed, but sometimes the continuum at time ¢ is referred to a 
convected or a fixed reference frame. Since tensor fields, such as stress and strain, 
which arise in the mechanics of continua generally have different transformation laws, 
the time derivative of these tensors depends not only upon the particular reference 
frame chosen, but is intimately tied in with the constitutive equations of the medium. 
Moreover, while the covariant derivative and the time derivative of a tensor commute 
in a material coordinate system, such an interchangeability of order of derivatives is 
not permissible in spatial or even in a convected coordinate system. 

The question of the time derivative of tensors, especially with reference to the stress 
tensor, i.e. the stress rate, has aroused much interest recently and has been discussed 
and investigated by Truesdell [1, 2], Noll [3], Oldroyd [4], Thomas [5], Cotter and Rivlin 
[6], and Prager [7]; the last, with the limitation to rectangular Cartesian coordinates, 
contains an account of the various definitions for the stress rate proposed in [1, 3, 4, 5, 6], 
as well as in others, including the earlier work of Jaumann [8]'. Although the expressions 
for the stress rate deduced in [1-6] differ from one another, their underlying objective is 
the requirement that the constitutive equations of the medium involving the stress 
rate must remain invariant under rigid motions. In addition, Prager [7] has pointed 
out the desirability for a definition of the stress rate, the vanishing of which will render 
the invariants of the stress tensor stationary. Also, mention should be made of a very 
recent paper by Sedov [9], dealing with the time derivative of tensors, which will be 
referred to again. 

The main purpose of the present investigation is to set forth a single general ex- 
pression for the time derivative (or the time rate) of tensors, valid in all curvilinear 
coordinates, and applicable to any tensor fields (representing the variables of state in 
mechanics of continua) whose transformation law is known a priori. Indeed, the time 
derivative of a tensor as derived here is necessarily dependent on the transformation 
connections (which oecur in the transformation law) of the tensor in question, and in 
particular obeys the same transformation law as the tensor itself. Furthermore, if the 
rate of a tensor (such as the stress rate), is defined in the context of the present paper, 
then (for a suitable representation of stress) the following conditions are fulfilled: (a) 
the stress rate vanishes when the medium executes rigid motion alone and when the 
stress is referred to a coordinate system participating in this motion, and (b) the rate 
of invariants of the stress tensor is stationary when the stress rate vanishes. This latter 


*Received August 24, 1960; revised manuscript received October 13, 1960. The results presented 
in this paper were obtained in the course of research sponsored by the Office of Naval Research under 
Contract Nonr-222 (69), Project NR-064-436, with the University of California, Berkeley. 

1As pointed out in [7], Noll’s result [3] is identical to that given much earlier by Jaumann [8]. 
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requirement is equivalent to that stipulated by Prager [7] when the coordinate system 
is rectangular Cartesian and the time derivative is taken in the sense of Jaumann. 

Following some preliminary background in Sec. 2, which includes certain results 
needed in the subsequent analysis, the expression for the time derivative of tensors 
mentioned above is deduced in Sec. 3 [Eq. (3.13)], where its various parts are con- 
veniently arranged in terms of symbols reminiscent of Christoffel symbols, and is 
applied to strain as well as stress. The reduction of the latter (i.e., the stress rate) to 
special cases and its comparison with the previous works, discussed in Sec. 4, are con- 
fined, for economy of space, to those in [2, 3, 4, 8, 9]. Finally, in analogy with the time 
derivative (3.13), an expression for a coordinate derivative of tensors is deduced in Sec. 
5, which has the properties that (i) when applied to a tensor, it has the same trans- 
formation Jaw as the tensor itself, and (ii) it commutes with the time derivative operator 
(3.13) in all coordinate systems. Aside from its possible future utility, this coordinate 
derivative of a tensor [Eq. (5.17)] is included here because of the manner in which its 
development parallels that of (3.13). 

2. Preliminaries. Let a generic point of a continuum 7 occupy initially 
t = t)) a position P with the material coordinates X*, and at time ¢ a position P’ with 
coordinates X’*; both positions, which may also be identified by position vectors R(X) 
and r(X’, ¢), respectively, are measured relative to the same arbitrary but fixed curvilinear 
reference frame X. It is clear that while X* and X’* refer to the same generic point of 
the continuum at two different times, they represent two distinct positions in space 


at time 


separated by a displacement u such that 

r=R+u. (2.1) 
If we also introduce an arbitrary curvilinear reference frame x and denote the coordinates 
of P’ in this frame by x’, then the spatial coordinates x‘ and X’* refer to the same generic 
point of 7 at time ¢, and the motion of the continuum relative to the fixed frame of 
reference X may be represented by 


Ox’ 


i i rE S a <é 
a’ = g(X",2); ax2|> 0, (2.2) 
or by transformations of the type 
ai = 2'(X", 1 
pare oon (2.3) 


x4 = X'“(x?, d), 


with the restriction that their Jacobians be positive. When only the mapping (2.2) and 
its inverse are admitted, then the kinematics of continua may be presented as in [10], 
a description which is also adopted here”. Except for some modifications, the notation 
employed is similar to that in [10]; here, however, covariant differentiation is designated 
by a stroke (|), while comma is reserved for partial differentiation. 

Because of (2.3), tensor transformations from xz‘ to X’* have the transformation 
connections 
ox’’ 
Ox" 


dx" 
ax’ ’ 


2As will become apparent later, the coordinates X’4 are introduced mainly for added clarity, al- 
though it will also offer some advantages. 


i 


LA = 





ree 
A"; = 


, 
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whereas the tensor transformation laws between x‘ and X“ may involve X4 , r4 [see 
Eqs. (2.7)] or both, as well as their inverses. In what follows, it is convenient to select 
any one of the three sets (X“, t), (X’“, t), and (z*, t) as independent variables. Thus, if 
y is representative of a tensor of any order, the partial derivatives of y with respect 
to the three space variables are related by expressions of the type 


vi= viaX4 = Via X's . (2.4) 


On the other hand, since the time ¢ is common to all three sets of variables, in taking 
partial derivative with respect to time, it may be difficult to recognize which of the 
three sets of variables is taken as independent. To avoid this ambiguity when writing 
partial derivative with respect to time, it is advisable to specify which coordinates are 
being held fixed by attaching an appropriate subscript to the partial derivative in 
question. Accordingly, the partial derivatives of y with respect to ¢ will be displayed as 


(2¢ ). . bX) + (24) (2.5a) 


v(t). (02), om 


and similar expressions for (dy/dt)x- and (dy/dt), . 
For future reference, we recall here that with the use of (2.1) and (2.4) the base vectors 
of the coordinate curves are related by 


Il 


A ryA , 
SE oe (2.6) 
g' = G4 = 2',.6", 
where 
rega—uwla, r= X44U*|,. (2.7) 


Further, if h denotes an arbitrary vector with components h‘ and H“ in the coordinate 
systems x’ and X“, respectively, then by (2.6) 


hi = r,H%, (2.8a) 
while 
hé |; = XJriH* |p . (2.8b) 
With the notation 
A A yA 
Ts = bp + L lp ’ (2.9) 


Rj = 6; —u'|;, 
where 6: is the Kronecker delta, it follows also from (2.7) that 
A Aw B 
_ T X ‘) 
ee (2.10) 
ra = Riazi, . 


It is instructive to consider here the velocity of a generic point of the continuum r. 
If we select the set. (X“, t) as independent and specify the position vector of this point 
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at time t by r = r(X%, #), and if we follow the motion of this point (by holding X“ fixed), 
then the velocity vector is simply 
or 
oh ae (2.11 
" (3) ) 


If, on the other hand, the set (X’“, ¢) is taken as independent, in which case r = r(X"", 2), 
then for an assigned value of X’* (which may be the final position in space), r = const., 
(dr/dt)x» = 0, and by (2.5), 


_ (8) _ pra, _ (2x) , 219 
y= (*). = Ys, = al os : (2.12) 


which is an expected result. 
While (2.12) represents the total (or absolute) velocity of P’ relative to the X frame, 
the total velocity of P’ relative to the x frame of reference is not (dx'/dt)x , but is given by 


Ox’ Ox’ 
ere) = (2.13 
; (22 ) (22) | ; 


which is easily verified by application of (2.5a) to x’. The first term on the right-hand 
side of (2.13) is the velocity of P’ relative to the x frame and will be denoted by 


Ox 
ee O43 
q (2) ?, 48) 


and in order to assign an interpretation to the second term on the right-hand side of 


Ox ox’ +) ™ 
<< we : (2.1% 
(% ). ( at J,’ 


where (0X’“/dt), is easily interpreted as the absolute velocity of a point fixed in the 
zx frame (i.e., a point, the z‘ coordinates of which are constant in time). Since (0X’"/dt), 
gives the components of this velocity in the X’ system, it follows from (2.15) that 
(dx'/dt)x- is the negative of this same velocity in the x frame. Introducing the notation 


Ox ; — 
( at ) Pp; (2.16) 


then the combination of (2.13) and (2.14) yields 
v=qt?D, (2.17) 


a well-known result signifying that the total velocity of P’ may be expressed as the sum 
of its velocity relative to the x frame and the velocity of the x frame itself or more 
precisely, the velocity of the point in the x frame with which P’ is instantaneously in 


(2.13), we first observe that 


contact. 
In the remainder of this section, we dispose of certain results preliminary to our 
main task. Selecting the set (2, t) as independent, then with the aid of (2.16), 


24.) a, wall (2.18 
(28. PP lh» sited 
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and by application of (2.5b) to X 4 , and then to X* we deduce 


ox ‘) yA 
1 = — ? _ mR 
( at ) X75 - (2.19) 
Recalling also that the Jacobian of transformation from x to X is given by [10, p. 63], 
J=(|T3| = |&2+ Us|, (2.20) 
then since (07}/dt)y = (8U“*|,/dt)x , by (2.1), (2.12), and (2.10) we arrive at 
a1) oJ (272) -~1,BryrA } 
= = J(T V* ‘ 2.21 
(2 t), = apa\ at ls (Ta B (2.21) 


Therefore, if in (2.8b) we replace H*|, by V“|, and introduce the result into (2.21) 


we obtain the identity 
21) 
- = Jy’ |, (2.22 
(2 xX ) 
which will be utilized in See. 3. 


3. The time derivative of tensors. For the sake of clarity, we first consider the time 
derivative of a class of tensors with transformation connections of the type r, , and 
confine attention, without loss in generality, to a mixed second order tensor which 


transforms according to 
1 Ai 
He =r Talk; ; (3.1) 


and which may be regarded to include the transformation law appropriate to velocity 
and displacement components (Sec. 2). 

Since there is no difficulty in defining the time derivative of Hj in the X frame, i.e. 
(0H} dt), , the main question confronting us is the definition of the corresponding 
quantity in the x frame or equivalently the establishment of an appropriate expression 
for the time derivative (or the time rate) of h} which will be designated by Dh;/Dt. 
Here this is achieved by seeking the tensor transformation of (@H7/dt)x in the sense of 
the original transformation (3.1), ie., Dh'/Dt is the desired time derivative provided 
it transforms according to 

Dh' en »( ate) , 3.2 

le Oe (3.2) 
To this end, we apply (@/0@)x to (3.1), multiply throughout by rir , and by virtue of 
(3.2) write the resulting expression in the form 


Dh (a) (2) | (2) | ‘ 
Dt ~ \at). + Vat) ~ Nar Bi (3.3) 


In (3.3), the coefficient of terms involving h; may be expressed as 


i or’ — 0 A. 
ri( ol ) 7 [2 (G 5) | (3.4) 


—— 28: ) 
~ 8 (28 I 
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and since by (2.5b), (2.14), and (2.18), 


0 i i i 9 - 
(2e:) = 8:59 +p |; & , (3.5) 
xX 
(3.4) becomes 
ri(%) "a { k hy co (3.6) 
“hie 4 Tr jr 


Motivated by the fact that the Christoffel symbol of the second kind has the form 


(ie 
j bf & Bik» 


\ 


we introduce the notation 


y = g'. 26: ie i k a | he Be 
()-0() (i rern, as 


and call (‘) the first time symbol. Substitution of (3.7) into (3.3) finally yields the 
J 


Dhi _ (2%) (*) (! 
ie "laws * hi — 


desired result, i.€., 
Ws } (3.8) 


which is reminiscent of that for covariant derivative of second order tensors. 

A generalization of the foregoing result is necessary, as not all tensors arising in 
mechanics of continua (especially in the constitutive equations) have transformation 
laws of the type (3.1). Notable among these are the stress o°’ and the strain e'’ which 
have transformations 


_ ” = re. 4 i. »4 “— : ’ \ 3 ‘: )) 
Eup ane XL aX. pei; ’ (3.10) 


S*” and E“*” being, respectively, the Kirchhoff tensor® and a measure of strain in the 
X frame. In order to distinguish between transformation connections of the type r4 
and x, , temporarily we use Latin indices for the former and Greek indices for the 
latter. For our present purposes, it is sufficient to consider a mixed tensor of weight n 
which transforms according to the law 
AT n_/ iv ia ‘ 
Hoa = oe ri X' a his . (3.11) 
To establish the time derivative of h‘, , i.e. Dh';,/Dt, as the tensor transformation 
of (0H%;/dt)x in the sense of (3.11), the procedure is the same as in the case of (3.8), 
except that in addition to terms involving (‘) there also arise expressions x*, (0X *;/dt)x 
J 
and (01 /dt)x for which use may be made of (2.19) and (2.22), respectively. In addition 
it is expedient to define the second time symbol by 


a a OX , a - 4 
B = ot (2X2) - 2 i EE (3.12) 


where (2.17) and (3.7) have been used. 


3In general, the stress tensor in the initial X frame has different representations. Here S4? is a 
symmetric tensor. 
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Thus, omitting the details, the result sought for the time derivative of a tensor with 
the transformation law (3.11) is given by 


Dhig eT {Has (282) 
— 8 = Jy jr?e* X 4( —*) = | —# 
Dt Olt /x Ot /x (3.13) 


+ (nas - ("nis + i . | fais + nhisv* |e , 


which includes (3.8) as a special case. With reference to the role played by the two time 


symbols (‘) and | in (3.13), it should be emphasized that the former is associated 
J 
with the transformation connection r; and the latter with 2°, . Having distinguished 


between the two transformation connections by utilizing Latin and Greek indices, in 
the remainder of the paper we dispense with the latter and use only Latin indices. 

Before proceeding further, we consider as an example, the application of (3.13) to 
the strain tensor 


e;; = 3u. |; tu; |, —u" |, a. |;). (3.14) 


Since by (3.10) the strain tensor has only transformation connections of the type 2‘, , 


then 
De;; _ de, ; 7 k . k 
Dit — ( ol ) [Ee LE fa . (3.15) 


Applying (3.13) also to u,|; , [which by virtue of (2.7) transforms by both r; and 
z*,|, and after computing (de;;/dt), and substituting the results into (3.15), simplifi- 
cation yields 

Pest = 40, |; +05 |) (3.16) 
for the rate of strain tensor. 

Although the application of (3.13) to the stress tensor will be postponed until the 
next. section, it should be emphasized that the time derivative operator defined by 
(3.13) is applicable to any tensor, whether relative or absolute, provided the trans- 
formation law for the tensor in question is known a priori. This should not be construed 
to imply that the application of (3.13) inherently provides for any special requirements 
which may be demanded of the time derivative of a particular tensor or of the con- 
stitutive equations containing it. Indeed, such requirements depend on the transfor- 
mation connections of the tensor which may even have more than one transformation 
law; this is exemplified by the representation of stress in the X frame which, in addition 
to S*”, may also be defined by 


p ge = Jr X Fe"). (3.17) 


4. The stress rate and comparisons with previous works. As pointed out in [2, 3, 5] 
and again in [7], an acceptable definition of stress rate must be one such that if the 
medium executes rigid motion only and if the stress is independent of time when referred 
to a coordinate system fixed in the medium (i.e., a coordinate system participating in 
the motion), then the stress rate must vanish when referred to this coordinate system. 
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It can be easily verified that the time derivative defined by (3.13) automatically 
satisfies the above criterion when applied to tensors which transform by x‘, alone, as 
e.g., in the case of stress with the transformation law (3.9), i.e. 


Do _ (aa° | a i “is | 
Dt — ( at ), * |e 5 H tov k. (4.1) 


For, if we let x’ be the coordinate system participating in the rigid motion of the medium, 
then by (2.14), g' = 0, and (2.17) and (3.12) become 


; 
v=?p, .| = 0, (4.2) 
, | 


respectively. Furthermore, since for a rigid motion v"|, = 0, the last term in (3.13) 
vanishes even for relative tensors. Hence, it follows at once from (2.5b) and (4.1) that 
if ¢ is a representative of a tensor field (such as o'’) which transforms by x‘, alone, then 
Do (2) (4.3) 
Dt at}, =" 
which meets the foregoing criterion, stated as condition (a) in Sec. 1. 
Moreover, if the rate of stress invariants (such as o'’ ¢;;) which arise in the theory of 
plasticity, is defined by application of (3.13), e.g., 
D (o'? ' etl De; ; (4.4) 
0 OG; —= 40 ° 
Dt Dt 
then, clearly, the vanishing of the stress rate will render the left-hand side of (4.4) 
stationary, and Prager’s requirement [condition (b) in Sec. 1] is fulfilled. 
If, instead of (3.9), the transformation law (3.17) is adopted for stress, then the 
application of (3.13) will not yield (4.1). To scrutinize this seemingly paradoxical result, 
we first observe that by combination of (3.9) and (3.17) 


> gaa = et lel (4.5) 


and that 


ere (25— (4.6a) 
Ay B a 
JX AX, Dt \ at ) 
; Dea"’ (2 2") 
: pA 8 = ; ( 6k 
Ir’2 Di ry ; t.6b) 


Next, in conformity with the requirement demanded of a suitable definition of stress 
rate, we assume at the outset that (dc'’/dt), = 0 in the spatial coordinate system (par- 
ticipating in the rigid motion of the medium) for which (4.2) holds, and proceed to 
determine the partial derivative with respect to ¢ of both >>**% and S**. Through con- 
sideration of the stress vector and its transformation relation involving the stress tensor 
and the use of (2.5a), it may be shown that 


0 » —) ( ) k CB 7 
a) =r re >, (4. 

( ot x na k . ; 4) 
with the aid of which application of (0/dt)x to (4.5) yields 


(28) =0. (4.8) 


ot 
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Hence, by the results (4.7) and (4.8), it is clear that of the two representations for stress 
in the X frame, only S*” with the transformation law (3.9) supplies a suitable stress rate*. 
The time derivative of stress with the transformation law (3.17) specifically leads to 


De" _ (i), 
Di ~~ \k/? 
which jibes with (4.6b) and (4.7). 


In the remainder of this section, we consider the various definitions of stress rate 
proposed in the recent literature all of which can be derived as special cases of (4.1). In 
this connection, it is pertinent to mention that since the stress rate (4.1) meets both 
conditions (a) and (b) discussed above, it follows that any definition of stress rate 
deduced from (4.1) will also satisfy these conditions. 

(a) Truesdell’s definition. If, as in [2], we select the x frame as one fixed in space, 
then p' = 0, q' = v', and by (3.12) 


H " { ope oe i (4.9) 


With the aid of (2.5a) and (4.9), (4.1) reduces to 


De q - . i - si ; wi . 
Dt * (*2-) +o" |e —v' po! —v' po" +o |, . (4.10) 


Observing that the first two terms on the right-hand side of (4.10) form the material 
derivative, designated in [2] by a superposed dot, then it becomes immediately apparent 
that (4.10) is identical with Truesdell’s Eq. (1.10) in [2]. 

(b) Jaumann’s definition. This definition [8], which has been rediscovered by Noll 
[3], may be deduced as a special case of (4.1) by considering the continuum to undergo 
rigid motion and referring the rate of stress to a coordinate system fixed in space. Under 
these conditions by (2.22), p’ = 0, and q’ = v' with 


v, |, =, ,0 | = 0; |. = 0, (4.11) 


1 
Ww; a(V; j 


— »,|;) being the vorticity. With the aid of (2.5a), (3.12) and (4.11), (4.1) 
reduces to 
Da’ (22:) i k k 
= {-— +o |,v —w' o,; + ow; , (4.12 
Dt at : hk kj ki ) 
the result given by Noll [3]. 
(c) Oldroyd’s definition. In order to establish the relationship between (3.13) and 
the corresponding expression in Oldroyd’s work [4], it is convenient to consider separately 
tensor fields which transform by x‘, and rj . For this purpose, we limit ourselves (with- 


out loss in generality) to 


og, wo OR (4.13) 
and 
hi = rar"H , (4.14) 


‘It is of interest to note that the transformation law for $4 is of the same type as that for strain 
given by (3.10). This is in keeping with the remark made in Sec. 1 that the choices of stress, strain and 
their rates are intimately tied in with the constitutive equations of the media. 

5In this connection, see the remarks made in [7] where other references employing Jaumann’s 


definition are cited. 
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which have, respectively, the representations 
Ss = Xiz' 50; (4.15) 


and 
Rew Foe Mi . (4.16) 
in the X coordinate system. 
Oldroyd’s approach in [4] is to obtain the transformation into the X’ system of 
{0/dt ( )], in the x frame; the latter frame (x) being taken as convected with the medium 
in which case g’ = 0, p* = v’, and by (3.7) and (3.12) 


(‘) =v'|,, H = 0. (4.17) 
J J 


In order to show the reduction of (3.13) to Oldroyd’s results for tensor fields with trans- 
formation connections of the type x‘, , we proceed to determine Do}/Dt and DS;,*/Dt, 
which (when x’ is convected with the medium), are related through 





DS; _ yr A.i Do; ; 
Di — X 4a B’ Dt . (4.18) 


Taking into account (4.17), (3.13) when applied to o*’ with the aid of (2.5b) yields 


De; _ (2) a (221) 
= (5). = (3) - (4.19) 


The definition (3.13) may also be used to calculate (DS%*/Dt) provided we identify 
zx’ in (3.13) with X’“ as utilized in this section. Denoting by V’“, Q’*, and P’* the counter- 
parts’ of v’, g’, and p’, respectively, in the X’ system which is fixed in space, it follows 
that P’* = 0, Q’* = V", and by (3.7) and (3.12) when 2° is identified with X’* 


a) ={ 
Hi 7 is “h 


Applying (3.13) to the combination of (3.7) and (4.15), and utilizing (4.20), (2.5a), and 


(4.19), we obtain 


0 , J A * ’ 1 r 4 r ! 7D r yes 
(921) X"2'» = (? se | +S) ViO+V' |p SH°+ VV"? |e So’. (4.21) 
Ot Ot /x: 


A 
BC 
A vy’ — v"* lees . 


“a (4.20) 
\ 
) 


The right-hand side of (4.21) is Oldroyd’s definition of stress rate which in [4] is denoted 
by (6S’3/6t). 

We now turn our attention to the reduction of (3.13) to Oldroyd’s result when the 
tensor in question transforms by r; . Following the procedure which led to (4.21), the 


‘It should be kept in mind that P’4 and Q’4 defined as counterparts of p‘ and q‘ in the X’ system, 
are not the transformations (X‘4 p‘) and (X‘4 q‘) which hold only if X’ and z are two different types of 
coordinate systems; in Sec. 4 of this paper these two coordinate systems are of the same type. 
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time derivatives of (4.14) and (4.16) with the aid of (4.17), (4.20), and (2.5a) may be 
written, respectively, as 


Dhi (i) + of [ht — oF |, hi (4.22) 
and 
Dae = QE) ate an 
which may be shown to be related through 
a Xz), = DH (4.24) 


Substitution of (4.22) and (4.23) into (4.24) leads to an equation which when put in 
the form 


e i + A 
(ani) X'iz'g: = (228) + H5' |e V’ — V" |p He? + VV" |p He (4.25) 
ot /, Ot /x: ; , 

is similar to (4.21) and agrees with the corresponding result in [4]. 

In identifying the prediction of (3.13) with Oldroyd’s results, no mention was made 
of the time derivative of relative tensors included in [4]. If the transformation between 
x and X involves weighting, then a term of the type v*|, will occur on both sides of (4.21) 
and (4.25) and will cancel each other. If the weighting is present in the transformation 
between x and X’, then an additional term, obtained from the determinant | X’+|, 
will be present on the right hand side of (4.21) and (4.25) and these will again be in 
agreement with those in [4] for relative tensors. 

(d) Sedov’s results. The results given in [9] are special cases of (3.13) in two respects: 
first, as in [4], the coordinate system is convected; and second, in view of the method of 
approach, the treatment is limited to tensors which transform by r, alone. On account 
of the latter, only the first time symbols arise which as in (c) above for convected co- 
ordinates are given by (4.17) and (4.20). 

Since the vanishing of the time derivative (D/Dt) of the base vectors together with 
(4.17) and (4.20) yield expressions of the type (6g,/dt), = v‘|.g, , then Sedov’s main 
results expressed in Eqs. (9) of Ref. [9] can be obtained by direct application of (3.13). 
Thus, for example if h'’ transforms by r, alone, we have 


DH’** aH’** f A 4 7,Dyzz/CB B ; rrDyyr AC 
Di = ( zw) +e oS} tie pa ase 


which are, respectively, the third and fourth Eqs. (9) in Ref. [9]. The first and second 
of Sedov’s Eqs. (9) may be obtained in a similar manner, while in order to deduce the 
last of Eqs. (9) in [9], we assume as Sedov has that the coordinate system x‘ is Cartesian 
in which case by (3.7), 


1 , 
‘] = g'- x ;= - 
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The remaining results in [9] are merely specializations of the above and are derived 
with the further assumption that the coordinate systems x* and X* coincide at some 
particular instant of time. 

Before closing this section it should be mentioned that (3.13) also may be shown to 
reduce to other definitions of stress rate, e.g., that given by Thomas [5], where a special 
Cartesian coordinate system called ‘“‘kinematically preferred” is employed. As this 
reduction is reasonably lengthy, it will not be included here. 

5. A coordinate derivative analogous to the time derivative (3.13). While the order 
of covariant derivatives of a tensor and [d(_ )/dt], is commutative in a material co- 
ordinate system, such as interchangeability of order of derivatives, though often useful, 
is not in general permissible. With the primary aim of providing a mechanism for commu- 
tation in the x frame of D(_ )/Dt and an appropriate coordinate derivative, we deduce 
in this section an expression for the latter which will be designated by double strokes 
(j|). 

As in Sec. 3, where the time derivative (3.13) is obtained as the tensor transformation 
of the total time derivative in the X frame, the coordinate derivative that we seek plays 
a similar role with respect to the covariant derivative, and again depends on the fact 
that not all tensors transform by the same law from X* to x‘. For our present purposes, 
it will suffice to consider separately the following two transformation laws, 1.e., 


Ri rH’, (5.1) 


and 


We first take the partial derivative of (5.1) and with the use of (2.9) write h', as 
h', = g';-G,H* + g'-G, .X"H* + r4H",X 5 . (5.3) 
Recognizing that the left-hand side of this equation and the first term on the right-hand 
side (with G,H* = g,h’), upon transfer, is the covariant derivative of h' and combining 


the remaining two terms in the form 


iw B A J A ) Cc 
ak ‘| 4 B + ‘ai pit | 
we arrive at 


h' ae. _ -* (5.4) 


as the tensor transformation of the covariant derivative of H*. But, more than this, 
the transformation law (5.4) also states that the derivative part of H*|, transforms 
by X * and not by r* . 

Before proceeding further, we recall [10] the vector sets 


c, = XiG, = Rig, 


(5.5) 
c’ = 2°,G* = (R")ig’ 
and the Cauchy-Green measure of deformation 
ci = c’-c’ = x',7',G*", oo 
(0.60) 


7A B 
Ci; = C; °C; = XX Gap ] 
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where (2.10) and (2.6) have been employed on the right-hand sides of (5.5). We now 
turn to (5.2), take its partial derivative, use (5.5) and rearrange the result in a manner 
similar to (5.4) to obtain 
i i k i wBQA} » 
¢o,;t+ec “Cy 50 = 2°,X‘8 |B - (5.7) 
In keeping with the original transformation law (5.2), the left-hand side of (5.7) is the 
tensor transformation of the covariant derivative S*|, , where in particular we note 
that its derivative part transforms by X * . 


IE $ |] =c'-¢,, (5.8) 
v J 


then the left-hand side of (5.7) may be designated by 


o ||; = o'; + IF lle (5.9) 


which is the appropriate tensor transformation of S*|, in the X frame. Since the coordi- 
nate derivative (5.9) reduces to the covariant derivative in the X frame, to be con- 
sistent, we employ throughout this section the notation ‘||; for the coordinate deriva- 
tive of any tensor regardless of its transformation law. In view of (5.1) and (5.4), as 
well as (5.2) and (5.7), it is understood that this coordinate derivative (1) is identical 
with the covariant derivative y'|; if ¥* transforms by r7 , and (2) it is defined by the 
right-hand side of (5.9) with o* replaced by y’, if ¥* transforms by X“ . 

In order to generalize (5.9) for the coordinate derivative of any tensor, we need to 


establish the relationship between iF Ai and the Christoffel symbol ‘ fe To 
jk ; 


bes 
this end, it will suffice to consider the transformation 


Introducing the notation 


i yA 
u' =r7r,U’,7 
the first coordinate derivative of which is 


ue ||; =u |; = 7X 5U* |e. (5.10) 


In contrast to the covariant derivative u'|;, , since in (5.10), one index transforms by r, 
and the other by X “ , the coordinate derivative u'||;, reads as 


u »,=u'd; li 


(a)... , - . 
=U Ilha ct lr “aft \|;- IF | fe lle (5.11) 


iy By CyyA | 
= raX ;X pl IBC 3» 
which is obtained by direct computation. 
By (5.5) and (5.8), as well as (2.9), 


i udu 2 ea | 
IF ll = sf ei, k oo (5.12) 
- ‘ a — (R™) Ju” |x. 
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On the other hand, if we solve for [, ml from (5.11) with the aid of (2.9), we obtain 


ell, ln} (eet 


Combination of (5.12) and (5.13), following a lengthy manipulation, yields the relation- 
y, La) J b 


ship between the two symbols IF ai and ‘ - in the form 


1 ae t ” 
if rll- - ott, ‘ sii 


which is similar to the difference between the two time symbols given by (3.12). 

The relation (5.14) may be used to advantage in obtaining the partial derivative 
of J = | T3|. Thus, with the aid of (2.10) and (5.14), 

Fue IE Ves (5.15) 
an expression similar to (2.22). 

With (5.4) and (5.7) as well as the result (5.15) at our disposal, the generalization 
of the coordinate deriv ative for any tensor field is immediate. Again, as in Sec. 3, employ- 
ing temporarily Latin indices for transformation by r¢ and Greek indices for trans- 
formation by X 4 , if H% transforms according to 

ia i Ba yA A ~ > 

J"hjs = Tl; . r X 3H, , (5.16) 

then it can be shown as in the case of (5.4) and (5.7) that the desired coordinate derivative 
has the form 


nit |, = his, 


+4! a bat, ws a 
(5.17) 


+ (Ll bs - [LJ be 


+ nhijgu™ || mk 
Dispensing again with the use of Greek indices, and considering the transformations 
of 


Dhi* 
mp ik ||  entded I I 
J*hji |\|m and J Di’ 


the truth of 


Dhit\ _ 
(75H) a = Dy Mt ln 6.18 
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can be verified. That the identity (5.18) holds is not surprising if we recall that (3.13) 
and (5.17) were deduced, respectively, as the tensor transformations of (0( )/dt)x and 
( )|a which commute in the X frame. 
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BOOK REVIEWS 


Solutions numériques des équations algébriques. By E. Durand. Volume 1. Masson et 

Cie., Paris, 1960. vii + 327 pp. $13.27. 

Whereas the present volume is concerned with the numerical solution of a single equation, a forth- 
coming second volume is to treat systems of equations. Three introductory chapters (90 pp.) present 
methods that apply to transcendental as well as algebraic equations (series development, continued 
fractions, iteration, inverse interpolation, etc.). The principal part of the work, which is concerned with 
polynominal equations begin with a chapter on computing with polynomials (20 pp.). There follow 
chapters on transformations of polynomial equations (35 pp.), location of the roots in the complex 
plane (25 pp.), methods using the powers of the roots (34 pp.), consecutive determination of the roots 
by iteration (45 pp.), and simultaneous determination of the roots by iteration (14 pp.). A brief chapter 
on double precision and floating point arithmetic (14 pp.) concludes the volume. 

Through out the book the exposition is pragmatical, making only very modest demands on the 
reader’s mathematical background. A wealth of numerical examples illustrate all important features 


of the methods. 
W. PRAGER 


Qualitative theory of differential equations. By V. V. Nemytskii and V. V. Stepanov. 
Princeton University Press, Princeton, New Jersey, 1960. viii + 523 pp. $12.50. 


This volume is a translation of an excellent two-volume Russian work on the qualitative and quanti- 
tative theory of differential equations. The first half of the book (Volume One of the original) is devoted 
to a study of the stability of solutions of linear and nonlinear differential vector equations of the form 
dz/dt = A(t)x + f(x, t), while the second half covers topological dynamics and ergodic theory. 

Both halves of the book are well-written and lucidly presented. The second part will be of most 
interest since the material of the first part is available in a number of contemporary books in the English 
language and in translations of the American Mathematical Society. 

RIcHARD BELLMAN 


Theory and solution of ordinary differential equations. By Donald Greenspan. The 

Macmillan Co., New York, 1960. viii + 148 pp. $5.50. 

This is an excellent outline for an introductory text in the theory of ordinary differential equations. 
The only serious weakness as far as the outline is concerned is the lack of matrix techniques in dealing 
with N-th order linear equations and linear systems. 

The author attempts to present in 147 pages, everything from basic concepts of point sets to 
linear equations to existence theory; from power series solutions to approximate solutions to Sturm- 


Liouville theory. It can’t be done. 
{ICHARD BELLMAN 


Mathematical tables. Volume 4: Tables of Weber parabolic cylinder functions and other 
functions of large argument. By L. Fox. Her Majesty’s Stationery Office, London, 
1960. iii + 40 pp. $1.85. 

The main table gives the Weber functions W(a, + zx), which occur in the solution of the two- 
dimensional wave equation in parabolic co-ordinates, in the range 2! = 0(0.005) 0.1, a = —10(1)10 to 
an accuracy of 8 significant figures. In addition, the exponential integral, sine and cosine integrals, 
Airy integrals, the error integral and the first six factorial functions are tabulated to 10 significant 
figures for large values of the argument z. Full interpolation facilities are provided. 


(Continued on p. 136) 
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ON THE EQUATIONS OF EXTENSIONAL MOTION 
OF CRYSTAL PLATES* 


BY 
R. D. MINDLIN anv H. L. COOPER** 


Columbia University 


1. Introduction. The Cauchy-Voigt [1] equations of extensional motion of aniso- 


tropic, elastic plates may be written as 
(L, + p®’)u + L;w = 0, (1) 
Lsu + (Ls + p®’)w = 0, 


where u and w are the components of displacement in the plane of x and z (the plane of 


the plate), 2 = 70/dt and the L, , p = 1, 3, 5, are the quadratic operators, 
a” o° a” 
L, = Tis a + 2715 a aa + Ys55 33 » 
OX Ox Oz 02 
a” a a 
Ls = ss = 2 + Dy = TF Ye 5 (2) 
On OX OZ Oz 
a a” a 
Ls = Ys = 2 + Ws + Ya) = OT Fas 2.2" 
O2 OX 02 Oz 
In (2), the y,, are the elastic stiffnesses relating the stress-resultants (7,) and the strains 
(€, according to 
T = Ype€a ; Yva = Yap ; P; q _ 1, 3, 5, (3) 
where 
Ou Ow Ou ow 
ak ee é =; ¢é = ry 
dx © Oz dz Ox 
In the case of an isotropic plate, 
V11 = Y33 > 2755 - Ta Fats Tis, * Jas = 0 
and Lamé’s [2] representation 
0g 0 0 0 
“u=—— “3 " w= = ¥ (4) 
Ox Oz Oz Ox 


reduces solutions of the displacement-equations of motion to solutions of 


yiVie + pe =0, Y55V7°*v + 2’ = 0. (5) 


In a study of the three-dimensional equations of anisotropic elasticity, G. F. 
[3] considered a special elastic symmetry which, in the present context, corresponds to 
an orthotropic plate (y:; = ys; = 0) with the additional condition 


= (vir — Ys5)(¥a3 — Yss) = (is + Ys) (6) 


*Received September 12, 1960. 
**Now at Queens College, New York, N. Y. 


Carrier 
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and showed that, if 


. = (yin = %)"" = a (Ys3 — ¥s5) ° ro ; - 
(7) 
5 a Aes aon gee — v.27 Ov 
u (Y33 Y55) az + (v1 Ys5) ae? 
solutions of P 
dy dy 2 2 2 
Y11 ax? + 33 Oz + pQ oe = 0, Y55V y + pQ y = 0 (8) 


yield, through (7), solutions of the displacement equations of motion. 

In this paper, we find the most general representation, of this type, for the Cauchy- 
Voigt equations; and we show, with the example of quartz, how to find orientations of 
crystal plates which have constants of elasticity that satisfy the conditions required for 
reduction of the equations of motion to the simpler forms. 

2. Generalization of Carrier’s representation in two dimensions. Let 

0 0 
Ne (a, = + Gis —Jx = Mix, 


Ox 0z 
(9) 


0 0 
w= (a, ax + ds33 2, = M:x, 


expecting that there will be two sets of constants a;; and, hente, two functions x. Then 


(1) become 


[(L, + p2°*)M, + L;M;)x = 0, (Lz + p2*)M; + L;M,]x = 0. (10) 


t 
Because of the presence of the term p®’, the only linear operator that can be a factor 
of the operator in brackets in the first of (10) is M, and in the second of (10) is M/;, . 
Hence, if the bracketed operators in (10) are each to be expressed as a product of a 
linear and a quadratic operator, 17, must be a factor of L;\/, and M, must be a factor 
of L;M, . Accordingly, we set 
M,M; = L;. (11) 
Upon equating coefficients of like derivatives in (11), we find 
43:93, = Y¥i5 » 433013 = Y35 » 44:33 + 43431 = Vis + Ys - (12) 
If (12) are satisfied, (10) become 
M\(L; + M3 + p2’)x = 0, M;(L; + Mi + pQ’)x = 0. (13) 
Now, we require the quadratic operators to be the same in both of (13). Hence, we set 


L,+M3;=L1,;+ Mi=N, say, (14) 


where 


3” o" Fa 
N = bigat 2bs aa, + ban 
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and the b, are constants to be determined. Then solutions of 
(N + p2’)x = 0 (15) 
will yield, through (9), solutions of the Cauchy-Voigt equations. 
Upon equating coefficients in (14), we obtain 
Yu + a3, = Y55 + ai, =b, 
Ys3 + dis = ¥ss + @33 = b; , (16) 
Vis + Q33031 = Yas + i103 = bs. 


Iquations (12) and (16) contain six relations, among the four a;; and the y,, , whose 
consistency requires two relations among the y,, . Upon elimination of the a;; from (12) 
and (16), these relations are found to be 


(Yu — Ys55)(Y33 - "Y55) = (Yis + Ys5)- o> (vis + Yss)'5 (A) 


(Yis + Yss)(Yis — 35) - Yas(¥u1 = V5) = V1s(¥33 = Y55)- 

If 4 Ys5 0, the second of A is satisfied identically and the first of A reduces 
to (6). 

Subject to conditions A, (12) and (16) may be solved for the a,; and b, . We note, 
first, that the first two of (12) and the last two of (16) are satisfied by a,; = as, if b; # 0 
and by a,, — @;; if b; = 0. Hence there are two sets of solutions: 

an = tris — 33) (vis = vss)" 
a33 = ty3sl(¥1 = Y33) (vis = v35))” 


Qi3 = a3, = +[lvis — vs)/(vi1 — Ya3))'”” (17) 


by = + (vis — ¥55)/C1u — ss) 
bs = yas + (vis — vis)/(Y¥11 — Ysa) 
bs = Nis + Ys5 
and 
Qi, = —A33 = +1[(yis — vas)/(¥11 — Yss))'”” 
Ay3 = +tys5[(¥11 — Ya3)/(vis — Y35)]'” 
Q3, = Fiyisl(vi — ¥s3)/(vis — ¥35)]'” (18) 
by = bs = ¥55 — (vis — Y35)/(¥11 — Yss) 
b, = 0 
in which either the upper or lower signs may be taken, without loss of generality. 


3. Consequences of conditions A. Conditions A take a simpler form when expressed 
in terms of the elastic compliances, s,, , defined by 


Ey = Bebe | Seg = Bey ; p,q = 1,3, 5. (19) 


Pp 








112 R. D. MINDLIN AND H. L. COOPER [Vol. XIX, No. 2 


and showed that, if 


/2 OY \v2 OW 

a (Yi1 = Yss). a = (Ys3 = Yss)" a 2 

(7) 
/ /2 fete) /2 re} 
w= (Vs3 — Yss). = + (vu = ss). ov ’ 
solutions of , 

ay Fi 2 P - 

Yi a. + 33 ~ + pXg = 0, YsV v + pQX’y = 0 (8) 


yield, through (7), solutions of the displacement equations of motion. 

In this paper, we find the most general representation, of this type, for the Cauchy- 
Voigt equations; and we show, with the example of quartz, how to find orientations of 
crystal plates which have constants of elasticity that satisfy the conditions required for 
reduction of the equations of motion to the simpler forms. 


2. Generalization of Carrier’s representation in two dimensions. Let 
fe) re) 
“u= (a, — + ais a = Mx, 


Ox Oz 
(9) 


fs) 0 
w= (as, = + ss 2, = Mx, 


expecting that there will be two sets of constants a,;; and, hence, two functions x. Then 
(1) become 
[(L, + pQ’)M, + L;M;)x = 0, ((L; + pQ*)M; + L;M,]x = 0. (10) 
Because of the presence of the term p’, the only linear operator that can be a factor 
of the operator in brackets in the first of (10) is M, and in the second of (10) is M/; . 
Hence, if the bracketed operators in (10) are each to be expressed as a product of a 
linear and a quadratic operator, 1, must be a factor of L;4/; and M; must be a factor 
of L;M, . Accordingly, we set 
MM, = L;. (11) 
Upon equating coefficients of like derivatives in (11), we find 
41:43; = Vis 5 Q33013 = Y35 » Q11:433 + 43031 = Vis + Yss - (12) 
If (12) are satisfied, (10) become 
ML, + M3 + p®)x = 0, M3(L; + Mi + p2’)x = 0. (13) 


Now, we require the quadratic operators to be the same in both of (13). Hence, we set 
L, + Mi; =1L1,+ Mi=N, say, (14) 


where 
a a 3 
N = 6 ox" + 2bs Ox dz + bs 02” 
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and the b, are constants to be determined. Then solutions of 
(N + p2’)x = 0 (15) 
will yield, through (9), solutions of the Cauchy-Voigt equations. 
Upon equating coefficients in (14), we obtain 
Yu + a3; = 55 + ai =b, 
Y33 + ais = ¥55 + 3s = b;, (16) 
Vis + 33431 = Yas + Qiidi3 = 0; . 


Equations (12) and (16) contain six relations, among the four a,; and the y,, , whose 
consistency requires two relations among the y,, . Upon elimination of the a,;; from (12) 
and (16), these relations are found to be 


9 


“Ca ¥55)(Y33 _ Yss) = (Vis + Yo)” = (vis + Y35) > (A) 


(Vis + Yss)(¥is — Y3s) = Yas(¥i1 — Yss) — Yis(¥33 — Yss)- 


If + Ys5 0, the second of A is satisfied identically and the first of A reduces 


to (6). 

Subject to conditions A, (12) and (16) may be solved for the a;; and b, . We note, 
first, that the first two of (12) and the last two of (16) are satisfied by a,; = as, if b; # 0 
and by a), — a; if b; = 0. Hence there are two sets of solutions: 


Qi = zyisllvi1 — ¥s3)/(vis — ¥33)]'” 
G33 = t+¥s5[(¥i1 — ¥ss) (vi5 = Yas)" : 


Gi3 = Gy, = +[(yis ~— Y35) (yu ed Y33) |" (17) 


b, = 1 + (vis — ¥35)/Cri — Ysa) 
bs = ¥33 + (vis — Y35)/(¥11 — Yaa) 
bs = N15 + ss 
and 
Ay, = —Q33 = +i[(vis — v35)/(v1 — Ysa)" 
Qi3 = tyssl(vi1 — Ya3)/(is — ‘Yas))'”” 
a3) = Fiyisl(yi1 — 33)/(vis — Y3s)]'”” (18) 
b, = bs = ¥55 — (vis — Yas)/(111 — Ys8) 
b = 0 


in which either the upper or lower signs may be taken, without loss of generality. 


3. Consequences of conditions A. Conditions A take a simpler form when expressed 
in terms of the elastic compliances, s,, , defined by 


6p = 8yqTe 3 Sng = Sos) =S P CQ. = 1,3, 5. (19) 
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From (19) and (3), 


Yoe = Bae 4, (20) 


where S,, is the cofactor of element s,, in the determinant A S,q|- When (20) are 


substituted into conditions A, the latter become 
81, + 833 = 2813 + 855 , $15 = $35 (A’) 
and we note that conditions A’ make 
A = 855(811833 — Sis = Bigs 


Then, with ¢ = +x/A”"”’ for (17) and y = + ix/A’” for (18), we find that (9) and (15) 
become 


dy dg Oy oy 
u = (813 — 833) =— + =) ae + 515 3 + (3,5 _ 811) a 
Ox Oz Ox Oz ? 
(21) 
0g ; . do oy Oy 
W = $15 ax + (813 — 811) 32 — (813 — S33) ae — S15 Oe 
and 
ay ; ay A ay pAde 
§33 ar? sata 2815 ar Oz + 81 oa P ‘ Ot ) 
0: Ox Oz Bss € (22 
ay , vy ay 
33 t+ aa = PSs 22 
OX Oz ot 
We now rotate to new axes x, , 24 according to 
4 = xecos¥, —zsin¥, , Z4 = rsin¥, +2 cosW, (23) 


in order to eliminate the cross term in the first of (22). The strain-stress relations, re- 
ferred to axes 2, , 24 are written as 

A A A ‘ ~ 

€ = SaTp 5 p,q = 1,3, 5, 


pa*p ?) 


where, with s = sin VW, ,c = cos WV, , 


A 4 4 fe \.2.2 ‘ 3 ‘ 3 

Si, = 81:C + 8338 + (2813 + 855)8C — 28,58C — 283,58 C, 
me . a ae i 9 2.2 9 3 9 3 
$33 = 81,8 + 833¢ + (28,3 + 85;)8'cC + 28,;8'¢ + 28358¢ , 


A 4 ‘ 2.2 2 2\2 \/.2 2 
855 = 4(8:; + 833 — 28,:3)8c + 8;;(C = iy + 4(s,5 ~ Bie = s )sc, 


> 


A 4 ee 4 4 t 2 
Sis = (81; + 833 — 855)8C + 813(8 +c) + (Sis — S35)(C 


— $x, 


835 = 2s,,8°¢ — 28338¢° + (28,3 + 855)(c’ — s°)sc + 8,;(38¢ — s*) + s3;(c' — 3s°c’), 
si; = 2s,,8c° — 283,8°c — (28,3 + 855)(c — s°)sc + 8,;(c' — 38°c’) + 855(38'c¢ — 5°). 
Now, conditions A’ make 

Sis = $13 » 835 = 855 , 81 + 833 = 8, + 833 (24) 


and if, in addition, we set 


tan 2V, = 28,;/(833 — 8),) (25) 
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we find also 


A A ‘ 
Si1 — 833 = (81, — 833) sec 2W, , (26) 


A A 
Ss = &, = 0. 


With (21) and (24)—(26), the components of displacement, referred to axes x, and 2, 


become 


0 / 1) 9 
i, = (si, _— S33) = e + (i's ae $1) a ’ 
OX, 024 


(27) 
A A dg A A Oy 
We = (813 — 81,;) —— — (813 — 833) = 
A 13 11 Oz, 13 33) Ox. 
Also, the differential equations (22) transform to 
dy P dg A.A A.A “¢ 
833 ar: + es az: = p(s1 1833 — 813813) FY ’ 
“A “A 
(28) 


ay 4 ay oy 
2.2 :* eco” 
OX 4 OZ, p ot 
If we were to convert from constants s,4 to y,4 by means of the inverse of (20), we 
should find that (28) are the same as (8) and that (27) are the same as (7), except fora 
constant factor. Thus, conditions A reduce the general anisotropic plate to an ortho- 
tropic plate satisfying Carrier’s condition (6). 
4. Slowness vector. Returning to Eqs. (1), we insert the displacements 
u = A exp [iw(tr + ¢2 — JD], w = C exp [iw(te + fg — I] (29) 
corresponding to straight-crested waves. The resulting characteristic equation, which 
must. be satisfied by the components &, ¢ of the slowness vector [4], is 
(yud + 2viséo + 556 — p)Cvss& + 2yaséo + Yast — op) 
= fyist + (Yis + Vs5)EF + Vast] = 0. (30) 
With &, ¢ as real Cartesian coordinates, (30) is, in general, the equation of a pair of four- 
lobed curves. However, if conditions A are satisfied, the discriminant of (30) is a perfect 
square and the curves are a circle, corresponding to the potential ¥, and an ellipse, 


corresponding to ¢, with principal axes oriented in accordance with (25). 
The discriminant of (30) is also a perfect square if 


II 


he 53) (V13 + ss) 2yis(Yis ~_ Y35)) (B) 
(Y33 — 55) (V1i3 + Y5s) = 2y35(¥35 = Y15)- 


In this ease the locus of the tip of the slowness vector is a pair of orthogonal ellipses 
associated with waves whose properties are described in the following section. 


5. Consequences of conditions B. In the case of conditions B, we first perform a 
rotation to new axes 


Lp = x cosV¥, — zsinW; , Ze = rsinV, +2zcosV¥z,. 
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Referred to these axes, the stress-strain relation is 
7) = You€e » 
where, with s = sin Vz ,c = cos WV;z, 
Yn = Wc + Y338° + 23 + Qys5)8'C = 4y58¢° — Ayss8°C, 
ves = Yu8* + vase" + 2Wlyis + 2yss)8'¢° + 4yiss’c + 4yasse’, 
Y55 = (vi + ¥33 — 2713)8'C° + Yss(C a 7 + Ams — Yas(C - s°)sc, 
vis = (yu + 33 — 4755)8°C° + yis(s* +c) + 2s — vas — 8°)se, 
ye, = ¥18'C — 338° + (y13 + 2¥s5)(c? — 8’)se + y15(38'¢ — 8°) + yas(¢\ — 38°C), 
ye, = ¥1180° — Ya38'¢ — (413 + 2yss)( — s*)se + yis(c’ — 38°”) + ya5(38'c — s‘). 
Now, if we set tan 2V, = (y:3 + ¥s5)/(¥3s — ¥:5) and apply conditions B, we find 
Vis = Yas = Vis + Y55 = 0. 


Hence, if uw, and w, are the components of displacement referred to axes xg and zz , the 
equations of motion become 





2 42 “2 “2 2 2 
B OUp 4 B ¢] Up ou B B re Wr B 0 Wp re) Wp (31) 
V3 33 Th as” Fee Ys 3 ~~ a7. ae t 
” OX OZB ot 4 OX B ” OZ, p dt 


and, consequently, the introduction of displacement-potentials g and y is unnecessary. 

Equations (31) have the curious property that, in the waves defined by them, the 
directions of the respective displacements are fixed: independent of the wave-normals. 
The mechanism which produces this phenomenon may be identified by noting that, 
under conditions B, 


= vsa[T1 + (vss 33) 73] (yiivss = 55755) » 
€3 = Yi T3 + (y55 )7 1] (yin¥ss ia 158735) 
Hence, since the conditions of positive-definiteness of the strain-energy-density require 
Vi1 > ©, 33 = 0, Yh > @, 1133 = 55755 A 0, 


both Poisson’s ratios, p? = — y%,/y2 and v3 = — 73/74, must be negative and in the 
ratios indicated. These are properties unlikely to be found in natural crystals. 

6. Application of conditions A to quartz. The constants, say s/’, of a plate cut from 
a crystal, at an arbitrary orientation with respect to axes of reference fixed in the crystal, 
may be expressed in terms of two angles of rotation, say 0 and 9%, and the constants, 


say s°,,, m,n = 1, --+ 6, associated with the axes of reference. When the expressions 
for the s’’ , in terms of the s°,, , 9 and ®, are inserted in conditions A’, in place of the 


constants s,, appearing there, a pair of simultaneous, transcendental equations in 0 
and ® is obtained. If there are any real roots of the pair of equations, they give the 
orientations of plates whose constants s/’ satisfy conditions A’. 
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In the case of quartz, the compliances associated with reference axes 2 , Yo , 20 , in 
accordance with the I.R.E. Standards [5], are 


0 0 0 
81; = 82 = 12.77, 82 = —1.79, 
0 > 0 0 < 
833 = 9.60, 813 = 83 = — 1.22, 
0 0 0 0 0 
Sn = 835 = 20.04, 8:14 = —8, = 385, = 4.50, 
0 0/..0 0 
Sc6 = 2(si1 = 812), 
0 0 —- —— ey oie oo 0 _ 0 = 0 —_ 
Sis = Sig = S25 = S26 = S34 = 835 = Sage = S45 = Sse = 0, 


. . . — > 2 j . 
where the numerical values, in units of 107'* (meter)’/newton, are those given by Bech- 
mann [6]. 
A rotation of coordinates about the axis of x» to axes 


2’ = Zo, y’ = yo cosO + z sin®O, 2’ = —ysinO + z% cosO 
gives formulas for the thirteen non-zero compliances s/,,, m,n = 1, «++ 6, in terms of 
the s°, and 0, as tabulated by Cady [7]. 


A further rotation, about the axis of 2’, to axes 
x’’ = 2’ cos®+ y’sin ®, y’’ = —2’sin @+ vy’ cos ®, z’=2 
gives for the six, of the twenty one, constants which we require: 
si! = sf, cos’ & + si. sin’ ® — (s/, + sho — 28/. — sf.) sin” ® cos ®, 


eff — of 
“33 “ 33 » 


sii = sf, sin? @ + s{, cos’ ®, (32) 


° 2 
all = gs, + (a, — 8;,) an ©, 


asi = sf, ain ©, 
si! = si, sin ® — (sf, — sf, — si) sin ® cos %, 
in which the s/,, are given in terms of the s°?,, and @ by Cady’s [7] formulas. 
We now suppose that all of the analysis in Secs. 2 and 3 was performed in the plane 
of x” and 2’’, so that the constants s,, , appearing there, are replaced by s//, . Then, when 


(32) are inserted in conditions A’, we find 
[s?,(3 sin? © — 4 cos’ O) + 2ksin*® O cos OQ] sin’ & 
— {s?,(4 sin‘ 0 — 3 cos’ O) + 2k sin O cos’ 9] sin & = 0, 
[2s°.(3 + cos’ Q) sin @ cosO — k sin‘ 0] sin‘ & 
+ 6[s°,(2 cos? @ — 1) sin® cos@ + ksin’ © cos’ O] sin’ ® 
— 2s°,(3 + sin? @) sin @ cosO — k cos 8 = 0, 
where k = s°, + 8, — 2s?, — s?, . These equations have the four pairs of roots 9 and 
® given in the first two rows of the accompanying table. The next six rows give the 
values of the constants s/’ computed from (32). The next row gives the values of 4 


computed from (25), with s,, replaced by s//; and the remaining rows give the values 


of the constants s,4 , referred to the axes x, and z, defined in (23). 
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The first column in the table describes the well known ‘“Z-cut”’ of quartz for which 
the Cauchy-Voigt equations reduce to isotropic form; i.e., the constants satisfy the 
conditions 


yy = gil off — Piel! of ee 
$11 = 833 855 = 2(S8i1 S13), S. =, = 


and Eqs. (4) and (5) apply. This cut has been studied in detail by Petrzilka [8], Bechmann 
[9] and Ekstein [10]. Modes of vibration of the remaining three cuts may be studied 
by means of Eqs. (27) and (28). 

As anticipated, no cut of quartz has constants of elasticity that satisfy conditions B. 


TABLE 


Orientations (QO, &, ¥) and compliances (spq, 10~"* meter?/newton) 
of the four cuts of quartz which satisfy conditions A. 


I II III IV 

e 90° —9.65° 18.79° 18.79° 
® 0 0 58.62° —58.62° 
sit 12.77 12.77 12.52 12.52 
$3 12.77 9.52 9.77 9.77 
855 29.12 23.97 37 23.27 
813 —1.79 —0.49 —0.49 —0.49 
835 0 0 —0.87 0.87 
rad 0 0 —0.87 0.87 
Vy 0 0 637° 73 .83° 
sti 12.77 12.77 12.77 9.52 
sf, 12.77 9 52 9.52 13:77 
sé 29.12 23.27 23.27 23.27 
si's —1.79 —0.49 —0.49 —0.49 
sf, 0 0 0 0 

sf 0 0 0 0 
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STRESS IN THE VICINITY OF A CRACK IN A THICK ELASTIC PLATE* 


BY 
M. LOWENGRUB 
The Department of Mathematics, Duke University 


1. Introduction. The determination of the distribution of stress in the vicinity of a 
crack plays a central part in recent theories of fracture [1], [2] and for that reason is of 
some technical importance. Much of this work is based upon an analysis of the stress 
near a circular or “‘penny-shaped”’ crack first discussed by Sneddon [3]. In that paper 
Sneddon reduced the problem of finding the stress near a penny-shaped crack in an 
infinite elastic solid to the determination of the solution of a pair of dual integral equa- 
tions. Green [4], recognizing the relation of this problem to that of the electrified disc 
problem of electrostatics (see, for instance, p. 175 of [5]) suggested a different approach 
leading to the solution of an integral equation. A solution of this crack problem using 
the ideas of Green and employing oblate spheroidal coordinates has also been given 
by Payne [6]. Still another approach involving this latter coordinate system is due to 
Sack [7]. 

All of these discussions are confined to the case of an infinite solid and only one of 
them—that due to Sneddon—appears to lend itself to generalization to the case where 
the crack is situated in a thick plate. Even then the calculations can only be carried 
out when the crack is lying in the central plane of the plate with its surfaces initially 
parallel to those of the plate. In this paper we carry out this extension of Sneddon’s 
method for the case in which the thickness of the plate is 6 times the diameter of the 
circular crack. 

As in Sneddon’s original paper it is assumed that the equations of the classical 
(infinitesimal) theory of elasticity hold. The boundary conditions appropriate to two 
types of problem are formulated in Sec. 2. In See. 3 the problem is then reduced to the 
solution of a pair of dual integral equations which are solved in Sec. 4 by a procedure 
due to Lebedev and Uflyand [8]. Exact analytical expressions are derived for the solution 
only in the case of large values of é—the solution is correct to order 6° in 6. The procedure 
to be followed in a numerical discussion of the solution is illustrated in Sec. 6 by the 
calculation of the shape of the deformed crack. 


2. The boundary value problems. We shall determine the distribution of stress in 
the vicinity of a penny-shaped crack in an elastic plate of finite thickness but infinite 
radius. The crack will be taken to lie in the central plane of the plate with its surfaces 
parallel to those of the plate. We shall assume that the deformation occurs as the result 
of the application of prescribed forces or displacements to the plane boundaries of the 
plate, these forces being distributed symmetrically about an axis which passes through 
the center of the crack and is perpendicular to the plane of the crack. If we take the 
radius of the crack to be a, the crack therefore occupies the region 0 < r < a, z = 0, 
where (r, 6, 2) are cylindrical coordinates referred to the center of the crack as origin 


*Received September 13, 1960. The results contained in this paper were obtained in the course of 
an investigation conducted under contract No. AF 18(600) 1341 with the Air Force Office of Scientific 
Research of the Air Research and Development Command. 
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and having the z-axis along the axis of symmetry of the problem. In this coordinate 
system we take the components of the displacement vector to be (u, 0, w)—the 6-com- 
ponent being identically zero because of the symmetry—and the non-vanishing com- 
ponents of the stress tensor to be ¢, , a, , ¢, , T,, - The thickness of the plate will be taken 
to be 2d. 

In the analysis of the boundary value problems we make use of a technique first 
introduced by Sneddon [3] whereby a crack problem for the thick plate — d < z < d 
is converted into a mixed boundary value problem for the elastic layer 0 < z < d. 

If we assume that the surfaces of the plate are given a uniform displacement e, then 
we have the boundary conditions 


onz = d: WwW =e, T,, = 0, yr > 0. (1) 
onz = 0: Tr: = 0, r>0; 

o, = 0, O<r<a 

w= 0, r> a. (2) 


The penultimate equation of this set is valid if there is no pressure applied to the surfaces 
of the crack, and the last one is a consequence of the prevailing symmetry. If we take the 
components of the displacement vector to be (u, 0, w + e«z/d) we see that u and w must 
be such as to lead to a displacement field satisfying the conditions 


onz = 0: w = 0, r>a; (3) 
t-, = 0, r>0; (4) 

l—vp che. 
= =< : — 1h <7 s 8: 5 
om 1 +90 — ® (‘)e. esrse ©) 
neg = d: w= 0, r > 0: (6) 
Ts = 0, F > 0. (7) 


In Eq. (3), v denotes Poisson’s ratio and E denotes Young’s modulus. 

On the other hand, if we assume that the free faces of the plate are pulled by a uniform 
tension F, we have the boundary conditions (2) again but the conditions (1) are replaced 
by: 


onz = d: oc, = —F, T,, = 0, ry 2G, (8) 


In this case, if we take the components of the displacement vector to be 


f (1 + »)(1 — 2p) (z 
4, 05 4 men 7) 


} 
we must choose the functions u, w to lead to a field satisfying the boundary conditions 
(3), (4), (7) and the conditions 

onz=0: o, =-F, O<r<a. (9) 


onz=d: o,=Q0, r > 0. (10) 


1961] STRESS IN THE VICINITY OF A CRACK 121 


3. The solution of the equilibrium equations. It is well-known ((9], p.452) that the 


expressions 
2 





oe eae ee 
,= (1 2v) ar az? (11) 
w=(1-— 2») 24 — Vo - ve! : (12) 


satisfy the equations of elastic equilibrium and the compatibility relations, provided 
that ¢ satisfies the biharmonic equation 


V‘o = 0. (13) 


Corresponding to these components of displacement we have the expressions 


_— =a E st f 9 oO ae) 
."Tea-et ad ‘¢) af (14) 

7 er 9 _- 0¢ P 
a (1 + »)(1 — 2p) \u ~) or (V'd) — or e.) (15) 


for two of the stress components. The remaining two are given by similar expressions. 
Now the function 


[ = - 14 (B + Aé€ SEF sinhé _ bs 
¢ a 3 WA + cosn é€¢ 
, ae 2v sinh — 6 +¢ 6 cosh é - 


— (B + Aéf) sinh E¢\ J o(Ep) dé, (16) 
where A and B are arbitrary functions of &, and 
p=r/a, t = 2/a, 6 = d/a (17) 
satisfies the biharmonic equation (13) and it is readily verified that it also leads to an 
expression for r,, which satisfies the boundary conditions (4) and (7). 


4. The dual integral equations. We shall consider now the first of the two boundary 
value problems discussed in Sec. 2. The expression (16) satisfies the boundary con- 
ditions (4) and (7). It remains to determine the functions A (¢), B() so that the conditions 
(3), (5) and (6) are also satisfied. 

If we substitute from Eq. (16) into Eq. (12) and make use of the condition (6), we 
find, after some manipulation, that, if the function (16) is such that the boundary 
condition (6) is satisfied then A (£), B(€) must be related through the equation 


(¢ 6 + 2vsinh — 6 cosh é 5) A(¢) = sinh’ ¢ SBE). 
We may therefore write 


2vs 
A= -(1-—»va'ett), B= -(1- ») £3 + 2sinh § 6 coshé 6 av?) 


sinh’ & 6 


to obtain the expression 


@ = —(1 — ved [ e*v t)[{2v + & coth — 6} cosh éf 
— {t 6 cosech’ — 6 + 2» cothé 6 + Ef} sinh Et] Jo(Ep) dé (18) 
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for the stress function ¢(r, z) occurring in Eqs. (11)-(15). Corresponding to this stress 


function we have the equations 


(1 — ve f* cosh &(6 — ¢) cosh §¢ 
= ET La — 2) PS _ gg 
Z 1 — 2y i | , sinh é 6 ” “ sinh’ & 6 
.. sinh &(6 — &) 
— ef sink & 3 [one dé, (19) 


21—v)yef sinh &¢ (€) J (Ep) 
eo — | sinh &(6 — ¢) + & 6 —>— — Ef cosh &(6 — £) at == dé (20) 
- de sinh & 6 sinh & 6 


by which to determine the displacement field. 
Further, if we substitute from Eq. (18) into Eq. (14) and put ¢ = 0, we find that 


: EA — ve “ 
(o.]. = ——— | (EVA (E 5) J (Ep) dé (21) 
L j (1 + v)(1 za 2v) d Jo pg s s s 
where 
ae — Sega .. £€6+sinhé 6 coshé 6 = 
K(é 6) = — {1 + AE 5)|] = € 6- no (22) 
sinh’ — 6 


If we substitute the conditions (3) and (5) into Eqs. (20) and (21) we find that the 


function ¥(£) must satisfy the dual integral equations 


| ell + HE d)]ve)J(es) dt = 8", O< p<, (23) 


| Y()Jo(te) dE=0, p>. (24) 

The second of the two boundary value problems can be treated in a similar way. 

If we substitute from Eq. (16) into Eq. (14) and make use of the condition (10) we find 
that, in this case, the relation between A(£) and B(é) is now 


(¢? 6° — Qvsinh’ ¢ 6)A + (E 6 + sinhé 6 coshé 5B = 0 


which may be written in the form 


(1 — 2)(1l+ )F (1 — 2)(1 + )F & 6 — Qvsinh’ é 6 
A) = = . x(&), Bg) = - © 9 ar 2 2—— (£) , 
E i £6 + sinhé 6 coshé 6 

The dual integral equations to be satisfied by x(£) are then 
| tM(é 8)x(é)Jo(to) dé=1, O<p<1, (25) 
| x()J(éo) de=0, p>, (26) 

where 

sinh’ § 6 — & & ™ 
= é (27) 


M(é 6) = 


£6+ coshé ésinhé 6 


5. Solution of the dual integral equations. In this section we shall derive the solution 
of the pair of dual integral equations (23), (24) using the method of Lebedev and Uflyand 
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[8]. Equations (25), (26) can be solved by a similar procedure. The essence of the Lebedev- 
Uflyand method is that it reduces the solution of a pair of dual integral equations to a 
Fredholm integral equation of the second kind which can then be solved by standard 
methods [10]. 

Since we are interested in values of 6 much greater than 1 we write the function 
H(é6) in the form 


H(é 6) = 2 _ (1 + 2mé d)e?”*’. (28) 


m=1 


Following Lebedev and Uflyand we write 
1 i 
v® = [ sin) dt, 0) = 0. (29) 


It is easily shown that the equation (24) is satisfied identically by this representation 
of ¥(£). Substituting from Eq. (29) into (23) we find that 


[ J (€p) dé [ £¢(t) sin (Et) dt + [ EH(E 5) J o(Ep) dé 


) 
al 


| o(t) sin &) dt = 8", O<p<i. (30) 


70 


Since (0 0 we see on integrating by parts that 
»1 1 
| £(t) sin (Et) dt = —o(1) cosé + [ ¢'(t) cos (Et) dt. 


Now 
f 0, 0<p< 
| J (Ep) cos (Et) dt = , ' (31) 


\(p — #)'”’, gag, 
so that, by a simple change in the orders of integration we find that the first term on 
the left side of Eq. (30) becomes 
ia ’ 


[ (p — t) "(0 at. 


70 


If in the second term on the left side of Eq. (30) we replace J,(é) by its integral repre- 
sentation 


(2/7) [ (p — u’)~*” cos (ué) du, 


we find that this term reduces to 


2° _du | _—— _ prep, — 2) 
VGuw i 1H*(u + 1) — H*(u — d]o(d dt, 


where 


H*(w) = 3 [ ” EH(& 8) sin (tw) dé. - 
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We thus obtain the equation 


e pe 
/ (p — t)7'°o6’(t) dt + (2/r 8) [ (p — wu)” du 
0 70 


al 


| o()K(u, )dt= 5", O<p<i1, (33) 


d 


where K,(u, ¢) = 3[H*(u + t) — H*(u — t)]. We now let 


p 


[ (p? — &)-9"(t) dt = fp); (34) 
0 


d 


> 


then it is well-known ((10], p. 39) that 


at 


2d , 1/2 
¢'(t) = - | v(f — v)7"*f(v) dv 


rdt Jo 
and since ¢(0) = 0 we may integrate to obtain the expression 
2° vf(v) dv . 
g(t) = = | a: (35) 
wTJo Vit —v) 


If we substitute this expression into Eq. (33), change the order of the integrations and 


write 


K(p,v) =v | (? —v)” dt | (p —u)’” &K,(u, t) du, o< 6 =<] (36) 


) 
we find that Eq. (33) reduces to the Fredholm equation 
a 7 el 
f(p) + 497? 5° f(v)K(p, v) dv = &", O < 5» < i. (37) 
By the classical theory of these equations ({10], p. 49) we can obtain a solution of 
this equation as a power series in (44% “6 -°) provided that 6 is sufficiently large for 
17°5° > M, where M is the bound on the kernel K(p, v). In this way we obtain the 
solution 


( “a ~ ae 
f(e) = 541 —40r? 5 | K(p, v) dv + (4x°* 5°)’ | K(p, v)K(v, w) dv dw + aa 3 
Since K(p, v) is of order 0 in 6 this means that 
nl 


f(e) = 641 — 40"? 5 | K(p, v) dvs + 0(67’). (38) 


v0 


To obtain the solution of Eq. (37) to this order we need only calculate 
ar al 


1 
[ K(p, v) dv = (po — w)"'” dw | t(6°K,(w, ¢)] dt. (39) 


v0 70 


Using the expansion (28) for H we find from Eq. (32) that 


H*(u) = 46? SS m {ull + (u/2m 87]? + ull + (u/2m 8)7)-*[3 — (u/2m 4)” ]} 
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so that 
K(w) = [ {e°K,(w, t)] dt = ps eh ~ kh +h = Wd, (40) 
where 
. t( + 0) it 3 3 ~ le 
I, 1 | ‘i+ | im +) ——. AT 7? = m* &'|tan™ {(w + 1)/2m 3} 
— tan™' (w/2m 6)] — 3m? & [1 + {(w + 1)/2m 6}*)"' 
‘ . tw — t) ) dt ee ‘es 9 
I, } i ll + {(w — )/2m a7) ~ m° S {tan {(u 1)/2m 6} 
— tan™' (w/2m 8)] + 3m? & [1 + {(w — 1)/2m 8; ]", 
, ff) tw + O13 — {(w + b/2m 4}*) dt , 
I : | il + iw +0 2m 6}")° = 2 {(w + 2)(w + 1)’ + 4wm’? 8} 


-{1 + [(w + 1)/2m 6)?}>? — (w* + 4wm’ 8°)[1 + w*/4 &m’)~"'), 


tw — O(3 — {(w — d/2m 4}") at ors 
[1 + {(w — t)/2m 5}7)3 = iw — 2)(w — 1)" + dwm* 6° 4 

-{1 + [(w — 1)/2m 6)?}>? — (w* + 4wm’ &)(1 + w?/4 8m’). 
Since the solution (38) is valid only to order 6° we replace J, , J. , 7; , 7, by their 


expansions in powers of 6 up to 6°. We find that 


9 « " 
cw) = 2 £3) — 3 6 2) \(5w* + 2) + 0(6* (41) 


where 
(3) = © m* = 1.202, ¢(5) = DY m™* = 1.0369. 
m1 m= 
The numerical values of the zeta-functions ¢(3), ¢(5) are taken from [11]. Substituting 
from Eqs. (40) and (41) into Eq. (38) we obtain the expression 
f(p) = 8° '[L — 4¢(3)/(Bmr 8°) + 3¢(5)(5p° + 4)/(20r 5°) + 0(6~*)]. (42) 


From Eqs. (35), (29) we find that 


y(—) = (2/n) [ sin &¢ at [ p(t — p)**f(p) dp 


70 


—(2/n)[fo D — $f, D*)E* sin E + 0(5"), (43) 


where 
fo = &'{1 — 4¢(3)/(3m 8°) + 35(5)/(Se 8°)}, fi = 35(5)/(4r 8) (44) 
and D denotes the operator d/dé. 


6. Shape of the deformed crack. We shall now proceed to find an approximate 
expression—correct to order 5 °—for the displacement w on z = 0,0 < p < 1. Putting 
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¢ = 0 in Eq. (20) we find that the normal component of the surface displacement is 
given by the equation 


21i-—»* f°, 
Wo(p) = —— ~e © | W(E) J o(Ep) dé 
l1—2Q Jo 


which by Eq. (29) and the result 


aX 


= 0, 
| J o(&p) sin (&t) dé = |‘ 
70 | 2 2\-1/2 


L(g — pgp) : t>p 


i <p; 


we find that 


al 


w(p) = 2(1 — v)*(1 — 2Qv)~"e ( — p)'o(t) at 
which, because of (35) can be written in the form 
al at 
wWo(p) = 4(1 — v)e/[r(1 — 2y)] | (( —p) dt | v(t — v)""*fv) dv. 


Substituting for f(v) from Eq. (42) and performing the integrations, we find that 


W,(p) = | Sis mo (sa —p)'” 


Tr | ] — Zyv) 


-{1 — [4¢(3)]/(me 6°) — [¢(5)(23 + 10p”)]/(30r 6°) + 0(877)} 


Using the numerical values of the zeta-functions we find that 


wo(p) = {4(1 — v)*e/[e O(1 — 2v)}}-(1 — p)'’*k(p, 4), (45) 


where “the correction term’”’ /(p, 6) is defined by the equation 


k(p, 6) = 1 — 0.510 6&* — 0.011 6-°(23 + 109°) + 0(67’). (46) 
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ON A PROBLEM OF MINIMUM WEIGHT DESIGN* 


By 
ZENON MROZ 
Institute of Basic Technical Research, Warsaw 


Summary. A problem of optimal design for perfectly plastic, isotropic structures is 
analyzed. It is shown that for such structures as plates or shells, an extremum of the 
volume, if it exists, is either a local maximum or a minimum. 

1. Condition for an extremum of the volume. Consider a region R of space that is 
bounded by a regular surface S. On the part S,; of S, let the non-vanishing surface 
tractions 7’; be prescribed, and on another part Sy of S, let the velocity u,; be required 
to vanish. It is assumed that S; and S, comprise the whole of the surface S. A rigid, 
perfectly plastic body B is to be designed to the following specifications: S; and Sy 
form parts of the surface S, of B; the remainder Sj; of Sz, is to be free from surface 
tractions, while S, is loaded in the prescribed manner and S, is rigidly supported; the 
body B is to be contained in R and should reach its load carrying capacity under the 
prescribed surface tractions; it is to have the minimum volume possible under these 
conditions. This problem has been studied by Drucker and Shield [4, 5]?; we shall briefly 
discuss it to lay the foundation for subsequent work. 

Consider a rigid, perfectly plastic body C that satisfies all conditions of the problem 
except that its volume V, need not represent a minimum. As this body is supposed to 
reach its load carrying capacity under the prescribed surface tractions, there exists a 
stress field o°, that is statically admissible for these loads and does nowhere exceed the 
yield point. Moreover, there exists a kinematically admissible velocity field uf the 
strain rates €°, of which are compatible with the stresses o{; and do not vanish identically. 
The principle of virtual velocities then furnishes the relationt 


[ i: a a / Tui dS; . (1) 


Denoting the surface of this body by S, = S7 + S, + S!, we consider a modification 
of the stress-free part S! of this surface, such, that the altered body remains at the 
yield point under the same surface tractions 7; . If we move S/ only outward (boundary 
E’ in Fig. 1), the load-carrying capacity, in general, will increase (cf. [15]). So, the varia- 
tion of S’ should be performed in such a way that a part of S/ is moved outside, the 
remainder inside of the body V, (boundary £ in Fig. 1). Further, we assume that there 
exists a stress field o* defined throughout the volume V* = V, + 6V of the modified 
body that is statically admissible for the prescribed surface tractions and does nowhere 
exceed the yield limit. By the first fundamental theorem of limit analysis ((14], p. 40), 
it follows from this assumption that the prescribed surface tractions cannot exceed the 
load carrying capacity of the modified body. We further assume that the velocity field 
u° is continued in a kinematically admissible manner into any volume that has been 

*Received Oct. 10, 1960. 

+Numbers in square brackets refer to the Bibliography at the end of the paper. 

tBody forces which do not alter the final conclusions, are not included here. 
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; 
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added to V, . It means that the field of strain rates «{,; derived from u{ satisfies every- 
where the incompressibility condition; moreover ¢;; are continuous functions within 
6V and on the surface S’ . Obviously, in general, there are many such extensions beyond 
the volume V, . If some additional kinematical constraints are imposed, the field of €; 
has to be compatible with them and the number of such extensions may be reduced, 
even to one possible field. (e.g. plates, shells etc. subject to thickness variation). In the 
volume 6V let of; denote a stress field that corresponds to the extended strain rate 
field by the theory of the plastic potential ({14], p. 15) excluding rigid regions, where 
o:; is not defined. Setting o* = of, + 60;; throughout the volume V*, and applying 
the principle of virtual velocities to the stresses o* and the strain rates «€{; , we obtain 
| o.¢;dV,+ | ai ;¢;; A(bV) + | 60,,¢;,dV* = | Tu, dS,. (2) 
Here, the power of plastic dissipation D = o*,e°; is completely specified by the strain 
rates ¢{; ({14], p. 37). 
Subtracting (1) from (2), we obtain: 


| D(e;;) d(6V) = —| be;,¢;; dV*. (3) 


The following theorem may now be proved:* Jf 0 < D (€;) = D = const. on S/ and 
D(é;) > D throughout V. , and if the velocity field u‘ can be so continued in a kinematically 
admissible manner beyond S’ that the strain rates of this continuation satisfy D(é.;) < D, 
the body C has the minimum volume that can be obtained if the given loads are not to exceed 
the load carrying capacity. Indeed, setting D(e:;) = D + AD, we write (3) in the form: 


D sv = -| éo,e5, dV* — | A Dd(sV). (4) 

By the principle of the maximum specific power of dissipation ({14], p. 37), we have 
éc,;€5; < 0, i.e. the first term on the right side of (4) is always positive and so is the 
second, since for points in the exterior of the body AD < 0, d(éV) > 0, while AD > 0, 
d(6V) < 0 for interior points. Thus, 6V > O whenever S’ is modified so that the pre- 


*An analogous theorem has been stated in Ref. [5]; it is not clear however, which parts of the body 
were assumed to vary. If the part S, may also undergo the variation, this theorem will no longer be true. 
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scribed surface tractions do not exceed the load carrying capacity of the modified body. 

A proof of this theorem may also be provided by means of the second theorem of 
limit analysis*. Indeed, considering the modified body V**, we assume that it has not 
yet reached or is just reaching its load carrying capacity. The kinematically admissible 
velocity field uS cannot therefore be unstable ({14], p. 40) and we have 


| D(é,) dV* > il Tui dS, . (5) 
Setting D = D + AD and comparing (1) and (5) one finds 
DiV>-— [ A Dd(éV) > 0, (6) 


which completes the proof. 

A special case of this theorem arises when D(e{;) = D = const. inside and outside 
of the body V, . Then AD = 0 and the right side of Eq. (4) contains only the first term. 
This case occurs, for instance, in dises, membranes or sandwich plates. If we consider 
a dise that is subjected to prescribed surface tractions 7’; on all of the cylindrical surface 
S, forming its ‘‘edge’’, only the lateral surfaces S/ , that is the thickness of the dise, 
may be varied to minimize the volume. Since stress and strain rate are constant across 
the thickness of the disc, the constancy of D on the lateral surfaces implies that 
De D = const everywhere. This design criterion for dises has been given in [4] 


7 


and [5]. 


Difficulties arise, however, when one deals with structures, such as plates and shells, 
for which D(€‘,) is greater in the exterior than in the interior of V, . The theorem proved 
above then is no longer applicable. Nevertheless, we shall try to choose the surface S/ 
so as to render D(e,) = D = const. on S’ . We shall see that this type of design has 
certain advantages. 

An infinitesimal modification of the surface S/ can be specified in terms of the normal 
displacement ng(/?) of the typical point P. Here, 7 is an infinitesimal constant, and 
¢(P) a continuous function of the position of the point P on S/ . If é¢,; is of the first 
order* (i.e. of the same order as n), then é0;,¢{; is negative and of the second or lower 
order according to whether the body V* is everywhere at the yield point or not. Thus, 
the first term on the right side of Eq. (4) is positive and of the second or lower order, 
while the second term is negative because now AD > 0, d(éV) > O in the exterior and 
AD < 0, d(é6V) < 0 in the interior of V, . This second term is likewise of the second 
order, since AD is a first order quantity in the vicinity of S’ . Thus, Eq. (4) implies the 
following: a body for which D(é.;) = D = const. on St! attains either a local minimum or 
a local maximum, depending on the sign of the expression on the right side of Eq. (4). 
This sign must therefore be studied in each individual case. It should be noted that, 
whereas the minimum is not necessarily analytical, the maximum must be analytical. 

It is interesting to note, that in the case considered last the second theorem of limit 

*The author is indebted to Professor W. Prager for this remark. In fact, the proof based on limit 
analysis theorems has been presented in Ref. [5). 

**The conditions that the function ¢(P) has to satisfy in order to ensure 60;; to be of the first order 
are not yet known in the general case. In the following, this problem will be discussed in detail for a 
circular plate, symmetrically loaded. 
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analysis does not provide sufficient information on the sign of 6V. Indeed, the inequality 
(6) now shows only that 6V is greater than a negative quantity of the second order. 
From this fact, no conclusions can be drawn regarding the sign or order of 6V. Equation 
(4), however, shows that we have an extremum of the volume. 

A somewhat different formulation of the problem of optimum design arises for 
bodies that are required to have a plane, axis, or center, of symmetry. Then any modi- 
fication of a part of the surface S’ entails the corresponding modification of the sym- 
metrically located part. If the prescribed surface tractions 7’; and the kinematically 
admissible velocity field u; need not have the symmetry properties of the body, the 
condition of constant D(«‘;) on S% must be replaced by the condition 


D(é,) + D”(é;) = D = const., (7) 


where D'’’ and D‘’ denote the values of D(e¢‘;) at a typical pair of symmetrically situated 
points of S/. 

2. Design of plates. Consider a solid plate of given plane-form that is subjected to 
prescribed transverse loads and supported in a given manner. A particular design of 
this plate is specified by the distribution of the plate thickness over the plane-form, 
and the discussion will be restricted to designs that reach their load carrying capacity 
under the given loads. The problem of finding the design of minimum volume has been 
thoroughly investigated for circular plates under rotationally symmetric loading [1, 3, 
6, 7] and some results have been obtained for other shapes or loads [2, 8, 9]. Drucker 
and Shield (4, 5) have established conditions under which the volume of a plate is station- 
ary with respect to neighboring designs. As we shall see, these conditions may lead to a 
local maximum as well as a local minimum. 

Let 2h. be the variable thickness of a design for which the given loads represent the 
load carrying capacity. Denote the rates of extension and shear at the lower surface of 
the plate by ¢'; , (¢, 7} = 1, 2), and assume that the specific power of dissipation has the 
constant value D at this surface: 


D(é.;) = D. (8) 


As is customary in plate theory, we shall assume the strain rates to vary linearly 
over the thickness of the plate. At the distance x, below the median plane of the plate, 
we therefore have the strain rates 


éij = (ts/hei;,  (t, 7 = 1, 2). (9) 


Since the dissipation function is homogeneous of the order one in the strain rates and 
does not change its value when the signs of all strain rate components are reversed, we 
have by (8) and (9) 


D(é;) = (| a3 |/h.) D. (10) 


I 


Equation (9) also represents the kinematically admissible extension of the strain rate 
field beyond the lateral surfaces of the plate that will be used in the following. 

The state of stress at a typical point of the fully plastic plate will be treated as plane, 
and the notation o'; , (7, 7 = 1, 2), will be used for the stresses in the lower half of the 
plate, those in the upper half having the same intensities but opposite signs. The same 
stresses are compatible with the extension (9) of the strain rate field. 
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If the variation of the plate thickness is denoted by 26h, and the element of area of 
the median plane by dA, the terms on the right side of (4) can be evaluated as follows: 


het+ibh 
[ b0.,¢5, dV* = 2 faa [be se8iCas/h.) des 
. 7 O 


(11) 
- | 50,,€°;[(he + dh)*/h,] dA, 
her+bh _ 
[ s Da(sV) = 2 / aA [ ((rs/h.) — 1] Dae, 
v Jhe (12) 
= D | (oh’/h,) aA. 
Equation (4) therefore furnishes 
6V = - [ (6h?/h.) dA — 1 D | 50;,€,;[(h. + 6h)*/h.] dA. (13) 


If the yield locus in the space with the rectangular Cartesian coordinates 
11 5 G2 , 012(2)' has continuously turning tangent plane, 60,;¢?; is negative and of the 
second order. The two terms on the right side of (13) are therefore of the second order 
and of opposite signs. This indicates that the volume of the original plate is stationary 
in comparison to neighboring plates that reach their load carrying capacity under the 
prescribed surface tractions, but no general statement can be made regarding the 
character of this stationary value (maximum, minimum, or saddle point). 

If, on the other hand, the yield locus has a plane face and the states of stress in the 
lower half of the plane are represented by interior points of this face, for the original as 
well as the modified plate, then ée,;¢;; = 0. The right-hand side of (13) therefore is 
negative and of the second order, so that the original plate represents an analytical 
local maximum of the volume. The same remark applies if the yield locus contains a 
straight segment and the states of stress in the lower half of the plate are represented 
by interior points of this segment for both the original and the modified plate. 

Finally, if the state of stress throughout the lower half of the plate is represented by 
points of an edge of the yield locus, ée;;¢;; will as a rule be negative and of the first 
order, so that the first term on the right side of (13) can be neglected. The original 
plate then represents a non-analytical local minimum of the volume. This statement, 
however, must be modified, if for each considered stress point o?, , 02. , o°,(2)' on an 
edge of the yield locus, the vector with the components e¢?, , €2> , €12(2)' happens to be 
normal to one of the faces intersecting in this edge. If then the corresponding stress 
point for the modified plate lies on this face, 60,;¢;; will be negative and of the second 
order; if, however, this stress point lies on the other face, 60;;¢;; will be negative and of 
the first order. In the first case, we have a stationary value, in the second a local minimum 


0 


of the volume. 
Equation (13) may also be rewritten in the following form: 


sv = —| (ah/h) aA - ‘ | 6M‘xi,[(h. + dh)/h,} dA, (14) 


where 61’, = 6c,,;h2 denotes the variation of bending moments due to the change of 
the stress components only and x‘; is the rate of curvature. Figure 2 represents the 
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Tresca yield condition referred to the principal bending moments. In view of the above, 
we conclude that for sides of the hexagon the optimal solution represents a maximum 
of the volume, while for the corners of the hexagon the plate will have either a minimal 
volume or its volume will correspond to a saddle point. The last case occurs, for instance, 
when x‘; coincides with one of the normals to the adjacent sides, say J at the corner A; 
the solution obtained will then represent a maximal volume for all states AB and will 
be a minimum with respect to all states AF, i.e. it will correspond to a saddle point. 

Example 1. Consider first a simply supported circular plate of the radius a under the 
uniformly distributed load g, and adopt the Tresca yield condition. On account of the 
rotational symmetry, the principal bending moments are the radial moment MW, and 
the circumferential moment M, ; the associated strain rates are the radial and circum- 
ferential rates of curvature x, and k, . 

Hopkins and Prager [1] treated this problem of optimum design assuming that the 
entire plate is at states of stress corresponding to the corner A of the hexagon (Fig. 2). 
For the plastic moment M, = o,h? (o, being yield stress in simple tension) they obtained 
the expressions 

M, = 1(@ — 7’), h. = 4(q/a)'"(a — r’)'”, (Fig. 3a). (15) 

The corresponding distribution of the rate of deflection w can be determined as 

follows. Starting from the condition D = const. we write 


a 2 > >. 172 
h, ; r dr adr (a —r’)'” 


c 


« Ma Me », dw a 
nae «x, + Myx, _ “is Ldw dw _ & a (16) 


where a = const. and & = 2a/(gc)’*. Integrating Eq. (16), one obtains 
» r 2 2\1/2 Fi 2 2\1/2 
w= a —aln-+(a@—r)” —aln {fa+(@ —r)”]/r} |. 
a 

It is easily verified that x, > 0 and x, > O within the plate, i.e. the solution really 
corresponds to the corner A. 

For the side AB we have x, = 0, i.e. d°w/dr* = 0 and w = A(r — a), hence x, = 
— 1/r dw/dr = A/r. The constancy of D yields then 


M.«,./h. = const. or h, = K,r, (Fig. 3b), (17) 
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(a) (b) 
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where K, is an arbitrary constant. If the side BC is considered, then x, = — x, ,w = 
A In a/r, and the condition D = const. leads to the solution h, = K.r’. Since for the 
states represented by CD a plate cannot deform, there is no solution for this region. 
For the side AB we have 


M,=0K%°=Cr, M,= (c os 2) OS, 


2/ 3 
Since 1, = 0 for r = a, C = q/2, M, = 0 within the plate, thus the stress state is 
represented by the point B. The volume of this plate is V = 0.946 2(q/co)' a® and is 
greater than that of a plate of constant thickness (V = 0.815 x(q ‘@)' a®). On the other 


hand it may be shown that no extremal solution exists for the state BC; thus plates corre- 
sponding to BC will have a greater volume than those corresponding to B, which there- 
fore represents a saddle point. 

An exact proof of the existence of a minimum for corners and a maximum for sides 
of the hexagon has been given in Ref. [7] in a somewhat different way. 

Example 2. Consider a circular plate symmetrically loaded and made of a material 
that obeys Mises’ yield condition. We then have 


M? —-M.M,+M,=M,, (M, + 6M‘)? — (M, + 6M‘)(M, + 642) 


(18) 
+ (M, + 6M’)? = M2, 
x, = (2M,-—M,), «, =M2M,—M,), (see [14)]), (19) 
and hence 
6M'x, + 6M’, = [8M“(2M, — M,) + 6M{(2M, — M.,)] (20) 


= —\(6M”? — 6M! 6Mi + 6MQ?). 
The quantities 61/’ and 6M{ can be determined from the yield condition (18) and the 
equation of equilibrium. Writing 


6h 


d | 6M.) — 8M. =0, 68M, = 2M, a + sM?, 8M, =2M,5*+ 6Mj, (21) 


dr 


we have 


d d bh bh 
6M!) — 6Mi = -2— (va, ) 2M, =: 22 
dr 8 OM?) dr a (22) 
Equations (18) and (22) determine 6M! and 6M? . From these equations it follows that 
6M’ and 6M will be of the same order as 6h, provided d(éh)/dr will be also of the same 
order. In view of (20) 64/7, x‘, is therefore of the second order and an extremum of the 


1 


volume is obtained. 
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Since the expression (14) represents the Weierstrass’ function for the volume of a 
plate, we can use 6M’ as the variation rather than 6h; note that 6M’ vanishes at the 
edge. We can then apply the known criterion for a minimum; setting M’ = dM,/dr, 
we have 

oF d°((M? — MM, + Md'r) 


= a 93) 
aM"? aM?’ (23) 


> 0 


, 


provided the Jacobi condition is satisfied. On performing the calculations, we obtain 
5M; + 4M,M, — 4M; > 0, (24) 


and hence 








Fic. 4 


lies within this region, it will correspond to the minimum of the volume. A problem of 
optimal design using Mises’ yield condition has been studied by Freiberger and Tekinalp 
[3]. It can be easily checked, that the solution obtained by them is contained within the 
angle P,OP, (Fig. 4). It is possible, however, that in the region P,OP, there is another 
solution, for which the volume attains its maximum or corresponds to a saddle point. 

If an extremal solution exists, which lies in the region P,OP, then there is no other 
solution in this region. Indeed, in the transition from one minimal solution to the other, 
we must pass the maximum, i.e. come into the region P,OP,; and then move back. This 
implies a contradiction, since we would have passed the same minimal state twice. The 
other minimal solution, if it exists, must therefore lie in the region OP;P, . Furthermore, 
since there is no upper bound for the voiume, if a maximum exists at all, some other 
minimal solution must also exist, unless we have a saddle point. 

3. Design of shells. Assuming that the thickness of a shell is to be altered sym- 
metrically with respect to the median surface, we cannot obtain D = const. on both 
surfaces. We should therefore apply the condition (7). Performing calculations similar 
to those presented for plates, we obtain the following results. When 


A = ojj(ej; + Kish.) + oi;(6; — «ish.) = const. (26) 


at every point of a shell, we have (neglecting orders higher than second) 


_ = o f prne : . = 
AsV = —-5 | Kio — 0,3) oh dA — | | be :sll, + 05) ae| dA, (27) 
=m JA JA he 
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where o,, , o,, denote the stresses at the two surfaces and ¢; , «{; are the strain rates 
and rates of curvature of the median surface. The first term in (27) is always positive, 
and the second, negative. Similarly as before, if 60;; ~ 6h, we conclude that condition 
(26) assures an extremum of the volume. Examples may also be provided which show 
that in some cases the extremal solutions are maxima. 

Problems of optimal design for shells have been studied in Refs. [10], [11], [12], [13]. 
In Ref. [12], the criterion (26) has been established as sufficient for a minimum of the 
volume, which is generally true for sandwich shells. 

If we assume that the thickness of the shell may be changed on one surface only, 
while the other surface is kept unaltered, we must take into consideration that the median 
surface is changed; the problem therefore requires separate treatment. If, for instance, 
D = const. on the surface that is changed, is smaller than D on the other surface, our 
shell will attain an absolute minimum according to the theorem proved above. 
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Theoretical physics in the twentieth century. Edited by M. Fierz and V. F. Weisskopf. 
Interscience Publishers, New York and London, 1960. x + 328 pp. $10.00. 


Pauli was considered by many physicists to have a deeper and keener insight into the fundamental 
problems of physics than anyone else in the past 50 years. This memorial volume to Wolfgang Pauli 
supplies the evidence for such a claim. There seems to be practically no field in physics to which Pauli 
has not contributed, and this account of his contributions written by experts makes extremely interesting 


reading. 
Roun TRUELL 


Special relativity. By W. Rindler. Interscience Publishers, Inc., New York, 1960. x + 186 
pp. $2.25. 

This is intended to be a concise introduction to the subject of special relatively and it is. It is a 
text that can be recommended for seniors and first year graduate students. The reviewer has found the 
problems and exercises from the first four chapters useful for seniors studying wave phenomena and 
radiation physics. The book assumes only that the reader is familiar with calculus and vector methods. 
Tensor theory is developed, to the extent needed, in an appendix. 

Roun TRUELL 


Annual review in automatic programming. Volume I: Papers read at the Working 
Conference on Automatic Programming of Digital Computers held at Brighton, 
1-3 April 1959. Edited by R. Goodman. Pergamon Press, New York, Oxford, London, 
Paris, 1960. xi + 300 pp. $10.00. 


Aside from papers on automatic coding systems for specific computers (mostly British) there are 
the following papers of more general character: Future trends in automatic programming (A. E. Glen- 
nie )—Some problems of a universal autocode (K. A. Redish)—The application of formula translation 
to automatic coding of ordinary differential equations (J. P. Cleave)—The philosophy of Programming 
(S. Gill)}—Automatie programming and business applications (G. Cushing)—Auto-programming for 
numerically controlled machine tools (J. E. Meggitt). 


Continuous geometry. By John von Neumann. Foreword by Israel Halperin. Princeton 
University Press, New Jersey, 1960. xi + 299 pp. $7.50. 


This is a carefully edited and corrected reprint of lecture notes multigraphed in 1935-7; it provides 
a definitive technical account of a fascinating generalization of projective geometry, invented by 
von Neumann. This generalization has applications to operators on Hilbert space, and hence it sheds 
some light on the foundations of quantum mechanics. 

No discussion of these applications is given in this book. However, as a display of brilliant algebraic 
techniques and incisive axiomatic analysis, by one of the most versatile pure and applied mathematicians 
of this century, the book is still extremely stimulating. In particular, it is one of the most penetrating 
studies available of the abstract algebra of subspaces of function space under intersection and union, 
and of the related algebra of semisimple and ‘‘regular’’ rings. 

GARRETT BIRKHOFF 
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HYDRODYNAMIC INERTIA COEFFICIENTS FOR A SLENDER BODY WITH 
A SAIL* 


By 
THEODORE R. GOODMAN AND RALPH GREIF** 


Allied Research Associates, Inc., Boston Massachusetts 


Introduction. In the course of studying the hydrodynamics of a submerged sub- 
marine it becomes necessary to calculate the hydrodynamic inertia coefficients of a 
body whose cross-section consists of a circle with one sail (fin). The simpler inertia 
coefficients for this shape appear in the open literature, although somewhat obscured 
by the presence of what, for this problem, would be considered extraneous fins. The 
inertia coefficient which gives rise to rolling moment due to roll velocity is not in the 
literature, and it is principally with the calculation of this coefficient that the present 
paper is concerned. 


Analysis. To determine the forces and moments on a rigid slender body moving 
through an infinite incompressible fluid where the motion of the fluid is irrotational 
and at rest far from the body, requires the evaluation of the six components of the sym- 
metric inertia tensor of the cross section of the body 


Nn. = Mm, = —p $ ¢, Ob ds j,k = 1, 2,3. (1) 
A on 
Refer to Lamb [1], Munk [2], Jones [3], Spreiter [4] and Bryson [5]. 

The cross section of the configuration is mapped conformally onto a circle in order 
to determine the velocity potentials for the low-speed two-dimensional motion of the 
cross section. The inertia coefficients may then be determined from Eq. (1). 

Consider a slender body of revolution of radius a(x;), bearing a sail of height 
t(v,) — a(x,), as shown in cross section in Fig. 1. Because of the plane of symmetry of 
the cross section, we see that: 


a) Mm; = 0. (2) 
That is, a motion in the x, direction cannot produce a rolling moment. In addition, 
b) my = 0. (3) 


That is, a motion in the x, direction cannot produce a force in the x, direction. 
The sail does not affect the inertia coefficient m.. since it lies along the streamlines 
of the flow. Therefore, the case of no sail suffices to obtain m.. . The result is 


Moo 


= 1. (4) 


pra” 


‘Received Oct. 20, 1960. This work was performed for the Office of Naval Research under Con- 
tract Nonr—2343(00) 


**Graduate student at Harvard University. 
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The expression for m,, may be taken directly from Summers [6] or Bryson [7]. The 


result is 
My; f a’ (t *1 (t\’ ft 
i = {t+ S(E-1)[(J 404 : (5) 
pra t 4t \a a a 
Equation (5) has been plotted in Fig. 2. 
t’ 
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Fic. 1. Mapping of a circle with a sail onto a circle (Eq. 9) 


The expression for m,,; may be taken directly from Bryson [7]. The result is* 
] 


(t/a — 1) mii. ji .) | 
(1 + 3 va * (t/a)? 


Miz (t/a — 1) 228 + m) (t/a + 1)" 
pra° T 16 (t/a)° * 4(t/a)'” 
where 
a. =x 8 om ois s (He = 1) 
Se OC CUMS 


Equation (6) has also been plotted in Fig. 2. 


*An error in sign has been corrected which makes my; positive. 


(6) 
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The apparent additional moment of inertia, m3; , for a body of revolution with a 
sail is not given in the literature and will now be obtained. 

[It is required to find a solution, ¢;(x, y), to the two-dimensional Laplace equation 
in the z = x + iy plane satisfying the boundary conditions 

on the sail 


0b; 7” 


= (2-0) Le = S e< ze < f, (7) 
OY 

on the body 
0 z ° > 
~ (r, coSa,r, sina) |,,-, = 0, (3) 
or, 


where (x, y) are a set of right handed, cartesian coordinates and the longitudinal di- 
mension is measured positive downstream from the body nose. Note that we have as- 
sumed unit angular velocity about this axis. 

The conformal transformation 


(1 +) —-2g9=¢+H4, (9) 
Z § 
where 
(1 + a) (10) 
¢ — 
q 4 i 
maps the cross section of the body in the 2(= x + ty) plane (Fig. 1) on the 
circle | ¢ | = g in the ¢ plane, the sail appearing as an are extending from — to 4 , 
where 
2 6a — t — (a 
cos 6) = es ee (11) 


g 2a + t + (a°/t) 
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On the sail the mapping function given by Eq. (9) and (10) reduces to 
az = [g(1 + cos 6) — a] + {g(1 + cos 6)-[g(1 + cos 6) — 2a}}'”. (12) 


The boundary conditions (7) and (8) transform to 


0o3 ' 2L sin 2 ‘ 

er r a a 27 _ 0 a) (le 

ar (89) lene = 7 — Cate <0<86 13) 
*) 1 ! 

~ (rr 9) Inne = 0 |O0|>H. (14) 


$ 


The problem has thus been reduced to an exterior Neumann problem for a circle and 
the solution for ¢, is given by 


alt 2 a 
o3(7; , 0) = g [ a In| 1 — <2 cos (@—y)+ (2) | dy. (15) 
0 [ 


9 Vane 
an. or; rr=g ry | 


On the circle ¢ = ge’’ the solution for the potential ¢, is given by 


a2 4 ; ~ | 
$3(9g, 6) = g | = In | 2 sin : = ¥ dy. (16) 


WT . ) or; re=g _— 


All integrals are to be evaluated in the Cauchy sense. Equations (13) and (14) are substi- 
tuted into (16), and z is eliminated by virtue of Eq. (12). Then, after an integration by 
parts to eliminate the logarithm there is obtained: 


2 ate f : : 6, 1/2 
¢;(9, 0) = —£ sin 6 | 4 cos ¥ (2 cos” — cos 2) (cos? ¥ — cos oe) (17) 


2 2 A ! > 2 6 ’ 
+ 2 cos’ ¥ (cos r - oe = f cos” y _ cos’ 5 dy, 


where cos 6, is defined in Eq. (11). Carrying out the integration it is found that 


¢3(g, 0) = ~# sin #\sin 0) + (0) + 2)(1 — cos 6) + cos 6) 
(18) 
6 || 
P ta 
(cos @ — cos >){I -F cos 6) | tan 3 + si 
sin 6 : 65 6 | 
tan _ tan 5 
Now, from Eq. (1) 
. Og; 
ae ae : ds (19 
M33 p $ 3 an ) 
and on the body surface 
ae 
Os ds = —x dx — y dy. (20) 


on 
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But Eq. (20) equals zero on the body surface, | z | = a; hence, the integration need 
only be confined to the sail. The resulting integral is quite tedious. The result is 


33 9 Pm 
e Lau(n + a, — 8n + 11’) 
7 pa rn’ | 2 

T 


+ 69 — 10n)[n(1 — »)]'? + 5 (14 — 259 + 159’ — 5n°)[n(1 — »)J'? ~— (21) 


bo 


2 


. oa 9 , 
+(1— |* O-igtes +4) + (2 — 20 on’) |, 


where 


3 
II 
Ss 
Il 
nn, 


1+ ) 
a 
Equation (21) has been plotted in Fig. 2. 

As a check on this result, it will be noted that the value of m,, is given by Miles for 
the case of a body of revolution bearing two fins of equal length which are mounted 
180 degrees apart [8]. For very small fins, the value of m;; given by Eq. (21) for the 
single-finned body should be one-half of the value given by Miles for the two-finned 
body, since, if the fins are small, they will not mutually interfere. Carrying out the 
limiting process for both cases bears out this relationship. An additional check may be 
made by proceeding to the limit as the body shrinks to zero. This case is given by Bryson 
[9], and it can be shown that Eq. (21) yields the correct limit.* 


REFERENCES 


1. H. Lamb, Hydrodynamics, 6th ed., Dover, New York, 1932, Chap. VI 

. M. M. Munk, The aerodynamic forces on airship hulls, N. A. C. A. Report No. 184, 1924 

R. T. Jones, Properties of low aspect ratio pointed wings at speeds below and above the speed of sound, 

N. A. C. A. Report No. 835, 1946 

4. J. R. Spreiter, The aerodynamic forces on slender plane and cruciform wing and body combinations, 
N. A. C. A. Report No. 962, 1950 

5. A. E. Bryson, Stability derivatives for a slender missile with application to a wing-body-vertical-tail 
configuration, J. Aeronaut. Sci. 20, No. 5, 297 (May 1953) 

6. R. G. Summers, On determining the additional apparent mass of a wing-body-vertical-tail cross section, 
Readers’ Forum, J. Aeronaut. Sci. 20, No. 12, 856 (Dec. 1953) 

7. A. E. Bryson, Comment on the stability derivatives of a wing-body-vertical-tail configuration, Readers’ 
Forum, J. Aeronaut. Sci. 21, No. 1, 59 (Jan. 1954) 

8. J. W. Miles, The potential theory of unsteady supersonic flow, Cambridge University Press, 1959, 
p. 171 

9. A. E. Bryson, Evaluation of the inertia coefficients of the cross section of a slender body, Reader’s 
Forum, J. Aeronaut. Sci. 21, No. 6, 424 (June 1954) 


*The authors are indebted to the referee for pointing out this limiting case. 
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Classical mechanics. By H. C. Corben and P. Stehle. John Wiley & Sons, Inc., New 
York, 1960. xi + 389 pp. $12.00. 


The first edition of this widely used text was reviewed on p. 334 of volume 9 of this Quarterly. 
The present edition contains many changes, the most important of which are as follows. A chapter on 
time-dependent forces and non-conservative motion has been added containing sections on the inverted 
pendulum, rocket motion, atmospheric drag, the Poynting-Robertson effect, and the damped oscillator. 
The presentation of the Hamilton-Jacobi theory has been made more concise and a chapter on special 
applications of this theory has been added containing sections on non-central forces, spin motion, 
variational principles in rocket motion, the Boltzmann and Navier-Stokes equations. A chapter on 
continuous media and fields has been added and the chapter on orbits of particles in high energy accel- 
erators has been considerably expanded. The new edition has four appendices dealing with Riemannian 
geometry, linear vector spaces, group theory and molecular vibration, and quaternions and Pauli 
matrices; the second of these replaces a chapter of the earlier edition. 

W. PRAGER 


Theoretical hydrodynamics. By L. M. Milne-Thomson. Fourth edition. The Macmillan 
Co., New York, 1960. xxviii + 660 pp. $11.00. 


According to a statement in the preface, this edition differs from the preceding one primarily by 
the addition of articles concerned with the following subject: the formulas of Plemilj, flows under 
gravity with a free surface, an exact treatment of the surface wave of constant form and the resulting 
“exact linearized theory,’ and comparison theorems including Serrin’s ‘“‘under-over’’ theorem. 


Theory of thermal stresses. By B. A. Boley and J. H. Weiner. John Wiley & Sons, Inc., 
New York, 1960. xv + 586 pp. $15.50. 


Interest in thermal stresses and their importance in engineering science have increased rapidly 
during the last decade. The book presents the extensive material in a clear manner. Written in a broad 
style, it is divided into two main parts. In the first part, the mechanical and thermodynamical founda- 
tions of the theory are laid. Although the laws of thermodynamics are introduced in a very general 
form, their application to the thermal stress problem is restricted from the outset to small deformations 
and to linear elastic and viscoelastic media. A number of uniqueness proofs for various cases is given 
in considerable detail. The—in general negligible—influence of inertia effects and of the coupling 
between the heat equation and the stress-strain equations is thoroughly and clearly discussed. Illustra- 
tive examples are given. Of course, most of the particular problems treated are of the uncoupled, quasi- 
static type. 

A long chapter—more than one hundred pages—is devoted to a presentation of the theory of heat 
conduction. In the opinion of the reviewer this inclusion makes the book unnecessarily bulky. The 
excellent treatise by Carslaw and Jaeger will have to be consulted anyway by everyone interested in 
thermal stress and, hence, in heat conduction. 

The second part of the book is, apart from a few cross-references, completely independent of the 
first. The basic equations are developed once more; this time, however, from a more intuitive point of 
view. First, only purely elastic systems are considered, coupling effects and, with two exceptions, 
inertia effects being neglected. Beams, built-up structures and plates are treated in detail. A chapter 
on thermoelastic stability is included. The final three chapters deal with inelastic systems and present 
a very good introduction to linear and nonlinear viscoelastic and to plastic stress analysis. 

It seems that the book, which to a considerable extent is based on original work of the authors and 
their coworkers, has been written with two different types of readers in mind: the scientist who wants 


(Continued on p. 173) 
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— NOTES — 


A NEW INTERPRETATION OF THE 
PLASTIC MINIMUM PRINCIPLES* 


By PHILIP G. HODGE, JR. (Jllinois Institute of Technology, Chicago, Illinois) 


In the theory of elasticity, the principle of minimum potential energy states that 


i = [ (0. de.,) av -/ Tu;dS — [ Fu; dV (1) 
JV « ST JV 


is a minimum for the actual state among all kinematically admissible states. A mathe- 
matically similar principle for elastic/plastic materials states that 
W(eZ , ¢) =~ W(oe’, e’) > 0,7 (2) 
where the functional W is defined for any independent states of stress o,; and strain ¢* by 
Wo ,€,) = [ oe* dV — [ Teu* dS — i Fu* dV (3) 
Jv Jsr Vv 
and primes denote differentiation with respect to t. We shall consistently use the asterisk 
to denote an arbitrary kinematically admissible state and unmarked symbols to denote 
the actual state. Precise definitions of terms, proofs of (2), and credits for original refer- 
ences may be found in any of several books such as [1, 2, 3, 4, 5). 

Despite the familiarity of (2), the principle has never received a name. This is doubt- 
less due to the fact that (2) has the unwieldy dimensions of rate of power and it is difficult 
to find any physical significance for W(o’, «’). The purpose of the present note is to 
show that the principle embodied in (2) may be given either of the following two in- 
terpretations. 


Theorem. Among all kinematically admissible rate states the actual rate state minimizes 
the time variations of 


, 


(1) the potential dissipation function W(c, €’); 
(2) the pseudo-potential” energy Tl. 
The proofs of these two assertions rest on the principle of virtual work and related 


principles involving rates. In the present notation 
W(c, e,') = We, igo © (4) 


where the superscript (¢) stands for any number of time differentiations. Also, we note 
that since the existing state is always known and only the rates are varied, only rate 
quantities need be distinguished by an asterisk. Therefore, we can write 


’ ox ° , ’ es 
ly — UW’ = Wo, ef) — We, ) = 0 (5) 
and 
” a ’ ” r/ ” i Ph on la" ’ 
I,’ — ll” = We, &{) — Wo, &’) + Wlef, ef) — Wo’, e’) > 0. (6) 
*Received May 9, 1960. This investigation was supported by the office of Naval Research. 
‘Numbers in brackets refer to the list of references collected at the end of the paper. 
?Since fo;;de;; may be partly dissipated for a plastic material, II as defined by (1) is no longer a 
potential energy. 
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Equation (5) is an immediate consequence of (4), and the inequality (6) follows from 
(5) and (2). The first statement of the theorem is thus verified and the second state- 
ment becomes evident when II(¢ + At) is expanded in a Taylor’s series in At about the 
time ?. 

If the complementary energy is defined by 


I. = | ( €;; dc.) dV — [ Tu; ( 
“ V * Sp 


Woo , €) =| o° es ~J Tout (8) 
V Sp 


and a functional W, by 


then the elastic principle of complementary energy states that II, is a minimum for the 
actual state among all statically admissible states and the analogous plastic principle 
states that 


W (a5 , 6) — Wo’, e’) > 0. (9) 
Just as for the first principle one can easily prove two consequences of (9). 


Theorem. Among all statically admissible rate states the actual rate stale minimizes 


the time variations of 


(1) the complementary dissipation function W(a’, e): 
(2) the complementary energy II. . 
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ON THE SIMULTANEOUS DIAGONALIZATION 
OF TWO SEMI-DEFINITE MATRICES* 


By ROBERT W. NEWCOMB (University of California, Berkeley) 


1. Introduction. The use of congruency transformations for simultaneously diago- 
nalizing two symmetric matrices, one of which is definite, is well known. One merely 
diagonalizes the definite matrix to (plus or minus) unity. This is then followed by an 
orthogonal transformation which diagonalizes the other matrix while preserving the 
unit matrix already obtained [1]. If, instead of being definite, one matrix is semi-definite, 


*Received July 11, 1960. 
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this method fails. However, if both matrices are semi-definite, this standard procedure 
can be extended. 

2. Diagonalization. In the following let the superscript ¢ denote matrix transposition 
and 1, denote the unit matrix of order r. Further let the rth order zero matrix be denoted 
by 0, . The main result is then the following theorem. 


Theorem: Let A and B be nzn real, symmetric, positive semi-definite matrices. 
Then there exists a real non-singular matrix 7’ and real diagonal matrices A, and B, , 
[see Eqs. (3) & (8)], such that 


II 


A=T'AT 
B = T'B,T 


Proof: Let A have rank a and B rank b and assume that b > a. We first find a 
real, non-singular 7, such that 


(1) 


A=T.A,T, 


(2) 
B =T" BT. 
where 
A, = diagonal [1, , 0,-.]. (3) 
If any of the last n — a diagonal elements of B’ are zero, the corresponding entire 
row and column of B’ are zero, since B’ is semi-definite. For the last n — a diagonal 


elements of B’ which are nonzero, we can reduce the remaining nondiagonal elements 
in these rows and columns to zero. We must do this by always adding the diagonal 
element to the off diagonal element in order to preserve A, . We can then write 


A = T'TIA,T,T, 


(4) 
B = T'T'B"'T,T, 
where 
B, 0 
B’ =| | In| | (5) 
: eon 


Here 8 > 0 is defined as the rank of B, . We now diagonalize B, by an orthogonal trans- 
formation 7, and put 


ae fei ll (6) 
ee Bice 
Now let 
f = 7,7,7,. (7) 


Then Eq. (1) results with 
B, = diagonal [A, , --- , As , Ou-g , 1s—s , On—o-0+8] (8) 
where A; > 0,2 = 1, --- , 8B. 
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By observing that neither the “sign” of A, nor that of By enters into the proof, we 
see that we can diagonalize two semi-definite matrices (possibly of opposite sign). We 
can also easily extend the theorem to Hermitian matrices. Thus let a superscript asterisk 
denote complex conjugation and let A and B be complex Hermitian, positive semi- 
definite matrices. Proceeding as above, but using complex 7, , 7; and unitary 7, , 
we can write 


(9) 


where A, and B, are as in Eqs. (3) and (8). 

3. Applications. The above theorem is necessary for the synthesis of networks which 
are passive or active at a point (to be published, for the basic concepts see [2]). It also 
can be used to advantage in the synthesis of two element kind networks, as well as in 
studying equivalent networks (see pp. 96 and 142 of [3]). Its use in studying the vibraticns 
of systems satisfying Lagrange’s equations should also be apparent. 
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AN UPPER BOUND ON RIGHT HALF PLANE ZEROS* 
By DOV HAZONY (Case Institute of Technology ) 


Abstract. An upper bound is placed on the number of right half plane zeros of 
functions of the type Z — m/n. Z and m/n are RLC and LC driving point impedance 
functions respectively. In addition, it is shown that if ReZ > 0 on j axis, the number 
of right half plane zeros is determined precisely. 

Introduction. In problems of control and network synthesis, it may be necessary 
to determine the number of right half plane zeros of certain impedance functions. In 
control problems, zeros in the right half plane may cause instability while in synthesis 
they may require active networks. In this paper an upper bound is placed on the number 
of these zeros of the class of functions Z — m/n and Z — n/m. These terms are defined 


below. 


Lemma. 


Given: I. Z is prf (an RLC driving point impedance function). 
II. m-+nisa Hurwitz polynomial, of degree d, of the complex variable 


S; m is an even and n is an odd function of S. 
*Received July 27, 1960. This paper, although based on the work sponsored by the U.S. Air Force, 
Cambridge Research Center, Bedford, Mass. Contract No. AF 19(604)3887, has not been approved or 


disapproved by that agency. 
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IIIa. The real part of Z, ReZ, may have some isolated zeros on the j 
axis but it cannot be identically zero along any segment of the j 
axis. (Part 1). 
or IIIb. ReZ > 0 on j axis. (Part 2). 


Then mz—n and m—nz have, at most, d zeros each in the right half S plane(Part 1) 
and exactly d zeros if IIIb is exposed (Part 2). 


Part 1. Proof. Let ¢ be given by 


ee 


¢ = (mZ — n)/(m — nZ) (1) 


€7)-er0te) 2 
J m-+n/\¢ + 1 


Since Z is prf (by 1) and |(m — n)/(m + n)| = {l on the axis (by II) 


irom W hich follows 


f— - 3 <1 onthe j axis. (3) 
Let 
(Z — 1)/(Z +1) — 


- = and .. |w!<1= onjaxis. 4 
(m — n)/(m + n) t+ 1 a y] (4) 


To prove the lemma, it will suffice to show that ¢ has as many pole zero pairs in the right 
half plane as w has poles, since w may not have more than d poles in the right half plane. 
From (4) the following can be verified 

Regt >0O for |w] <1, 

Reg =0 for |w|=1, (5) 

Re¢é <0 for |w|>1. 
Let D be any S plane domain such that (Fig. 1) 

Ww | — 


| along the boundaries of D, 
w|> 1 inside D. (6) 


Then, by the maximum modulus theorem, w has at least one pole inside D. Let ¢ be 
given by | ¢ | e'*. Then by (5) and (6) 


a | = 90° along the boundary of D, 


a | > 90° inside D, ( 


~I 
— 


a | < 90° outside D. 


Now, since the 90° phase locus must start at a zero and end at a pole, it follows that all 
poles and zeros of ¢ lie on the boundary of D and occur in pairs. Further, any phase locus 
starts at a zero and ends at a pole. Then, inside D, ¢ has as many 180° phase loci as the 
number of pole zero pairs it has on the boundary of D. Along each 180° phase loci the 
magnitude of ¢ increases monotenically* from 0 to ©. Then, at one point along this 


*This is due to the orthogonal properties of magnitude and phase plots. 
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Iw| = 2 


Iw) <1 











Fic. 1 A plot of w in the S plane. 


line, £ = — 1 and w has a pole. It follows that ¢ has as many pole zero pairs on the 
boundary of D as the number of poles w has inside. In addition, the following properties 
of w in the right half S plane may be ascertained: 


w has d poles at most. (by (4)) 

If w has poles it also has at least one domain D. 

This D must lie exclusively in the right half plane. 

(by postulate III and because | w | < 1 on the j axis, by (4)) (8) 


It follows that ¢ has d poles and d zeros in the right half S plane at most. This completes 
the proof (Part 1). 
Part 2. Rouche’s theorem [1] will be used to prove Part 2. 


Rouche’s theorem. If P(s) and Q(s) are analytic interior to a simple closed Jordan 
curve C and if they are continuous on C and 


| P(s) | < | Q(s)| on C (9) 
then the function F(s) = P(s) + Q(s) has the same number of zeros interior to C as 
does Q(s). 

Let Q(s) = 4 (n — m)/(m + n), Q(s) has d right half plane zeros and 
| Q(s) | = } on j axis (by II). 


and P(s)= 4 (Z — 1)/(Z7 + 1). | P(s) | < 4 onj axis as ReZ > 0 there 
from which it follows that F(s) = (nZ — m)/(m + n) (Z + 1) and therefore also nZ — m 
have d right half plane zeros. Similarly, by letting Q(s) = 3 (m — n)/(m +n), mZ — n 
has d right half plane zeros. 

This completes the proof. 


Applications. Example 1. In the following network (Fig. 2) Z = prf, m/n = a 
reactance function (m + n is of degree d). Then, the transfer function, 


_ F, _ nZ—m 
Zu = E, |\1,<0 ~ nZ+m (10) 


has, at most, d right half plane zeros. 
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Fia. 2 A lattice in which E,! 


nZ—m 


E,\1, -0 nZ +m 





Example 2. In the following network [2] (Fig. 3) ¢ is given by (Z is prf) 


(1 a Pg bZ(a) —_ aZ( DY; _@ Z(a)Z(b) b? —_ a” bs | 
} ab bZ(b) — aZ(a) 


{( nm Ss bz) — aZ{e)) _¢ Z v-@ i 
ab bZ(a) — aZ(b) ~ ab bZ(a) — aZ(b) 


b2(a)-a2(b) 
b-a ( 
=, 


:- 5 ~ _ (b2(a)-a2(b))? 


, i a 
. b2(a)-a2(b (b° -a~)(bZ(b)-aZ(a)) 


ab bZ(b) — ada) | 




















Fic. 3 Physical realization of an impedance function, Z, in terms of ¢. 


where a and b are positive real or complex conjugates with a positive real part. ¢ is 
constructed to have the factor (s — a) (s — b) common in both the numerator and 
denominator. Further, [2] Re ¢ > 0 on the ; axis. To complete the proof of the realiza- 
bility of ¢ it is necessary to show that no additional right half plane zeros and poles 
are present. By the above lemma, both numerator and denominator have only two 


right half plane zeros each. This completes the proof. 
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THE DISTRIBUTION OF QUADRATIC FUNCTIONALS* 
By JOHN B. THOMAS (Princeton University) ann EUGENE WONG (IBM Research Center ) 


1. Introduction. The purpose of this note is to consider the joint probability 
distribution of quadratic functionals of the form 


Yn = [ ys 2j(r)hin’(r)a,(7) dr, (1) 


where the x;(7) are correlated Gaussian processes and the h{?’(r) are given functions. 
The probability distribution of a single functional of the form of (1) has been treated 
by another method Siegert [1] for the case where the z,(7) are stationary. 

The solution to the problem described reduces to the solution of a homogeneous 
matrix integral equation. This equation shows explicitly the transformation functions 
hi’ and the covariance functions of the 2;(7). 

2. Analysis. Let the functionals under consideration be denoted by 


> 


Ym = | X(r)-hi™ (7) X(17) dr, m=1,2,--- ,p. (2) 


Here X(r) is an n-element column vector whose elements 2; are correlated Gaussian 


processes with zero means, and h’”’(7) is ann XK n symmetric matrix [2] whose elements 
h ‘x’ are given functions. In addition let X(7) be written formally as the series 


X(r) = DY ag,(7), (3) 


where the a, are random variables and the ¢,(r) are non-random vector functions. 
With the substitution of this series, (2) becomes 


Ym = > > a,a, | (7) -h”’(r)d,(7) dr. (4) 
1 k=1 v0 


The characteristic function for the joint probability distribution of the y,, is given by 
( Pp } 
. A ot ° = 
F(m ,°*: , 7) = E\exp a Mmm f (5) 
\ m=1 
where i = (—1)'”’ and E{ } denotes the expectation of the bracketed quantity. From 
(4), the characteristic function becomes 


. aa, | (7) -w(7)g,(7) ar}, (6) 


a 
Fim , +++, 7) = Ejexpi >» 
\ I7=1 k=l /70 


where the order of summation has been interchanged and the weighting function matrix 
w(r) has been defined by 


w(r) + Zz njh(r). (7) 
7=1 
Let the elements of the covariance function matrix of the z;(r) be given by 
Ry(r, 0) = E{x;(7)ax(0)}, (8) 


*Received Aug. 12, 1960. 
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and let the ¢,,(r) be the eigenfunctions of the homogeneous matrix integral equation 


t 


R(7, a)w(a)o(e) do 


Ad(7). (9) 


Then it is shown in the next section that the ¢,,(7) are orthonormal with respect to the 
weighting function w(r), Le., 


af 


| (7) -w(r)d,(7) dr 


II 


Sars (10) 


and the a, are uncorrelated; hence the a, are independent Gaussian variables with zero 
means and variances \, . The characteristic function can be written from (6) and (9) as 


Fim, , +++ 5) = kl exp l > ai} = Il (1 — 2id,)7'”. (11) 


Equation (11) is the solution to the problem provided that the series representation 
(3) is valid and that the eigenvalues A, can be found from (9). These points will now be 
discussed. 

3. Series representation. It will be shown that the vector eigenfunctions of (9) 
which correspond to different eigenvalues are orthogonal with respect to the weighting 
function matrix w(r). This proof parallels the standard proof for the scalar case [3]. 

or two distinet eigenvalues \, and \; , expression (9) can be used to form 


| at al 
| o,(7)-w(r)o,(7) dr = " R(7, c)w(c)d,(c) -w(7)¢,(7) do dr, (12) 
and 
| ,(7)-w(r)o;(7) dr = : o,(7r)-w(r)R(7, o)w(o)d,(o) dr de. (13) 


The matrix identity 
MA-B = A-M,B, 


where the subscript 7 is used to denote the transpose of the n X n matrix M and where 
A and B are n-element column vectors, may be applied to the right side of (12) to yield 


. [ o.(c)-wr(o)R7(r, c)w(r);(7) do dr. (14) 


T):W(7T)d(7) dt = 
| ¢ vu(r)d,(7r) d val 


Now, if o and 7 are interchanged in (14) and if it is noted that w;(r) = w(r) and that 
Rr(o, 7 Ri(r, oc), the result is 


[ o,(7)-w(r)o;(7) dr = : | o,(7)-w(r)R(7, c)wlo)d;(c) dr do. (15) 
* ) k AL) 


A comparison of (13) and (15) shows that, for A; # & , 


[ d1(7)-w(n\6,() dr = 0 (16) 


and, with proper normalization, (10) follows. 
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That the a, of (3) are uncorrelated as a consequence of (9) may be shown by forming 
E{a,a;}. From (10), 


E{a,a;} = BA / | [ * [X (1) -w( abel DIX(@) -w0(0)6,(0)] dr ao}. (17) 
The integrand of (17) may be rearranged to yield 
E{a,a;} = [ i E{ X(1)X 7(c) }w(o)¢;(0) ae | (4.00) dr. (18) 
Since E{X(r)X7(c)} = R(r, oc), it follows from (9) that 
E{a,a;} = ; [9.42 (rel dr, (19) 
and, from (10), that 
E{aa;} = ; 6 . (20) 


Thus the a, are uncorrelated random variables with variances \, . The means of the 
a, are each zero, i.e. 
a 
Eia,| = | E{ X(r)}-w(r)¢,(7) dr = 0, (21) 
Jo 
and the linearity of (10) makes the a, Gaussian. 
The series expansion (3) is a generalization of the Karhunen-Loeve representation 
[4, 5] and reduces to that representation for the case of one-dimension and for w(x) = 1. 
As in the one-dimensional case, it can be shown [6, 7] that the ¢; exist and that the 
components of the vector series (3) converge mean square to the components of the 
vector process X(r). In addition, if R(r, a) is a positive definite kernel, then the ¢; are 
complete and the \; have at most a finite degeneracy. In these proofs a generalization 
of Mercer’s theorem is required to represent the elements of the matrix R(r, ¢). 
4. The single-functional case. For the special case where only a single functional 
y is considered, (9) can be written as 


nt 


| R(r, a)h(c)o(c) do = vd(r), (22) 


and the characteristic function (11) becomes 


F(n) = [] (1 — 2im,)°’”, (23) 
k=1 
where ny, = A, . Now if the matrix g(r, yu; n) is defined by 
2 (Ll) (m) 
ide (T)de” (nu) 
a(t, a3) & > ee 24) 
Jim Ty BK; n) > 1 ne Zinn, , ( 1) 


then, by direct evaluation, 


g(r, 3) — in [ R(r, a)h(c)g(o, w; n) do = R(t, yp). (25) 


J ( 
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For the stationary case (25) is identical, except for a change in notation, with Eq. (79) 
of Ref. [1]. In terms of the matrix g(r, u; 7) the characteristic function can be written 


F(n) = exp i / dé [ Tr{g(o, 0; )h(o)] do, (26) 
where 77r[ | is the trace of the matrix. 
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FURTHER EXTENSIONS OF SCHUSTER’S INTEGRAL* 
By E. T. KORNHAUSER** (H. H. Wills Physics Laboratory, University of Bristol) 


The integral, 


[= [ (C’(x) + S*(x)] dz, 
where C(x) and S(x) are Fresnel integrals defined by 


C(z) = [ cos f° dt, 


S(x) = [ sin f° dt, 


was conjectured by Schuster’ to have the value (1/8)*. Proof that I does in fact have 
this value was given by Hardy’ and more elegantly by Ingham’*. More recently Bateman* 
has extended Ingham’s treatment to evaluate integrals of the form 


[ C(a)C(az) dex, [ C(x) S(ax) de, 


*Received September 1, 1960. 

**On leave from Brown University. 
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etc. An additional generalization, which is of interest in the theory of imperfect dif- 
fraction gratings’, involves the inclusion of additive rather than multiplicative constants 
in the argument of the Fresnel functions. 


Consider the integral 


I, = | C(x + a) dx = | C(x) dx = [xC*(x) |? + 2 | x cos x C(x) dx 


70 va Ja 


a @ 


—aC’(a) + [sin 2°C(x) |? + sin x cosa dx 


va 


II 


= —aC’(a) — sin a’C(a) + 4 2°‘ S(2'”’a). 


In like manner 
2 


i. = | S’(2 + a) dx = —aS*(a) + cos a’S(a) — 42°'7S(2'7a). 


Furthermore, 


he f (Ce+o+ Set ola=1+h, 


= —alC*(a) + S*(a)] + cos a’S(a) — sin aC(a). 


Changing the sign of a yields 


I, = | [(C"(e# — a) + ot a)| dx 
= a[C*(—a) + S*(—a)] + cosaS(—a) — sina C(—a), 
and since C(— a) = (x/2)? — C(a) and similarly for S(— a), 
I; + I, = (/2)'*(cos a® — sin a’) + 2a(m/2)'”* [(x/2)'* — Cla) — S(a)]. 


Now consider the integral, 


I, = | (Cia + a)C(x — a) + S(x + a)S(ax — a) dx. 


‘ 


Successive integrations by parts yield 


1/2 2 ° 2 ¥/ 2 
I; = 3(x/2)’(ecosa — sina)+a | [C(x + a) cos (x — a) 


7“) 
— C(x — a) cos (x + a)? + S(x + a) sin (x — a)? — S(x — a) sin (x + a)'| ae, 


which may be written 


) ) ‘ ) =. . : ; . ' 
I, = 39 2)'(cos a’ — sina’) +a | 5 Cia + a)C(x — a) + SQ + a)S(x — a)} dx 
9 oa : 


Thus it is clear that J, obeys the differential equation, 


dl, 2 > . 2 ’ 
1 = 4(4/2)'(cos a” — sina’) = F(a) 
da ° 


I; —c¢ 
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which may be solved by quadrature, 


I,=«a | u"F(u) du = 4(r/2)'”*{cos a’ — sina’ — 2a[C(a) + S(a)]}. 


Finally, the most significant integral from the point of view of the physical problem 


n= | {(C(w +a) — Ce - a)? + [S@ + a) — SQ — a)]*} de 
= J],+ I], — 2/1; = za, 


a remarkably simple result. 


NOTE ON THE POINCARE BOUNDARY-VALUE PROBLEM* 


By E. E. JONES (University of Nottingham, England) 


1. This note is concerned with the solution of a modified form of the Poincaré 
boundary-value problem [1]. It is required to solve the Poisson differential equation 
Vo f(x, y) for ,(x, y) defined in S; , the region enclosed by the circle C of equation 
/), such that on C 


a. (2 r+ 


0g; 


0 , 
k i + l1—~ + m¢; = g(x,y), (1) 
on Os " 
where /, 1, m are constants, f(x, y) is prescribed in S, , g(x, y) is prescribed on C, and 


d/dn, 8/ds denote differentiations along the inward normal and positive tangential 
directions respectively to C. 

It is assumed that ¢, = ¢;9 + ®; , where ¢,o(z, y) isa particular integral of the Poisson 
equation reflecting all the singularities of the complete solution ¢,; , and (2, y) is 
harmonic in S, , being together with its first partial derivatives single-valued and con- 
tinuous in S, . It is thus possible to write &, = re W,(z), where W,(z) is a regular func- 


tion of zin S, . If z = ae’ = ¢, 2 = a’/g on C, and by definition 
I in _, | Adin 
h(¢) = =e = k a + an + MP? io ’ (2) 
. or a 006 
then 
1 ¢ (+ 2)h(0) 
hia) = 5 [ SEP ay, (3) 
271 Joc : 2 


is the Schwarz integral representation of a function h,(z) regular in S; with a real part 
equal to h(¢) on C—for this to be so it is necessary for h(¢) to satisfy the Lipschitz 
condition on C [2]. The boundary condition (1) then takes the form 


ili a ) 
wil dW (¢) ~ wie — his) = 0, (4) 


la dt 


*Received September 30, 1959. 
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where a = a/(k — il). When re(ma) < 0, then a unique solution exists for W,;(z) in 
S,; in the form 
1 
Wie) =a [ hi(et)t-"*~ dt. (5) 


D 

This form of expression has an advantage over that given in [1], since it also includes 
the conditions for the existence of a solution. Thus when m = 0 the equation (5) de- 
termines W,(z) to within an arbitrary constant, provided the integral is uniformly 
convergent, which necessitates h;(0) = 0. It thus follows that since the integral of (5) 
is uniformly convergent for re(ma) < 0, and for m = 0 if h;(0) = 0, then W,(z) isa 
regular function of z in S,; , as required, and W‘(z) can be obtained by differentiation 
under the integral sign [3]. 

Repeated integration by parts extends the range to all positive values of re(ma), 
including ma a positive integer—in this latter case solutions exist but are not uniquely 
determined [1]. 

2. The particular integral ¢,) can be determined by integrating the Poisson differential 
equation in the form 


and the solution can be written as 
Gio = Yo + re 0, (2), 
where ;,(z) is a function of z, regular outside C, having prescribed singularities in S, 
only. It follows that (2) can be written as 
sf ae i # 3) ( 2. a (a 
h(g) = h¥(g) + rej —= ON] + m2,\—]?- 
\ a \$ $ 
But h(¢) = {re h;(z)}c, then evidently 


hz) = hX@ — e. a,(“) + mo,(~) , 
az zZ 2 


~ 


which is a regular function of z in S; , as required. On substitution into (5), and inte- 


= a,(“) > a. [ a( er dt. (6) 
a Z a J; z 


This result determines W,(z) when the singularities of 2;(z) are prescribed, and extends 
the results of Ludford et al. [4]. 

A similar type of reasoning leads to the conclusion that it is possible to add the 
real part of any function of z regular in S, to ¢;o without changing the value of ¢, . 

3. For the external problem with ¢, defined in S, , the region | z | > a, and satisfying 
a boundary condition similar to (1) with 0/dn denoting differentiation along the normal 
to C drawn into S, , the solution analogous to (6) is 


- 2 Snbx? 1 2 : 
W.2) = W*(2) + : a,(“) + = [ ny dt. (7) 


a Z 


grating by parts, then readily 


Wiz) = W*@) + 


N 


Here 


W*(2) = -a | h*(et)t"*— dt, 
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where h*(z) is regular in S, , including infinity, such that on C 


. 09%, l d¢*, 
bh ¥* (+ . = * 
re h*(z) = g(x, y) — k a 2 mo*, , 


and 2,(z) is regular in S, , having singularities only in 8S, . 
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ON MATRIX DIFFERENTIAL EQUATIONS* 
By JOHN JONES, JR., (/nstitute of Technology, Air University) 


The purpose of this note is to obtain a necessary condition and a sufficient condition, 
of an algebraic nature, for the matrix differential equation 


AX" + XB = C, (1) 


where all matrices considered are n-square, to have a single-valued solution. Capital 
letters denote matrices and prime denotes the derivative with respect to z. The elements 
of A, B, C belong to the polynomial domain &[z] of the field $ of real numbers. I is the 
identity matrix. 


Theorem 1. If a solution matrix X of Eq. (1) exists, then the following pair of matrices 
[A C [A 0 
| i. J (2) 
10 BI 10 B 


are equivalent. 
Proof. Clearly 


[f -X\(A c){7 -x") _[A ¢ -Ax" seas fd in 0 


= | (3) 
lo rjilo Bio «sj to B J 0 B 


and so the matrices of (2) are equivalent. 


Theorem 2. If the matrices of (2) are similar and there exist non-singular constant 
matrices P, Q such that PAP™', QBQ™' are diagonal matrices exhibiting the invariant 
factors a; ,71 = 1, 2, --- , a, a: ul b;,j = 1,2, --- , 8 of A, B, respectively, along the 


*Received October 12, 1959; revised manuscript received August 23, 1960. 
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main diagonal, then there exists a solution X of Eq. (1) for-— 7 <a<r<b<+o 
except at the zeros of the invariant factors of A. 
Since the matrices of (2) are similar, 


P 0\(A cl(P* 0) _ {PAP* PCQ")_ » 

|| | = | + 
lo Qgilo BiLo @Q- 0 QBQ"] - 
P O\;A OjfF 0 ect 


| 
| 
| 
} | 


_ far? 8 | 
0 gio Billo @} | 0 @BQ"| 


and the matrices J/, N are similar. Form the equation 
DU" + UF =G, (5) 
where D = PAP", F = QBQ", G = PCQ"™. Consider Eq. (5) element-wise, namely, 
4; 8) + uiihis = Gi: ; (7@,j7 = 1,2,---,n), (6) 


where u/) ¢ U"",u;; e U, d;; ¢ D, f;; « F, and g;; ¢ G. 

Next to show that Eq. (6) always has a solution. W. E. Roth* has shown that if 
the matrices of (2) are equivalent, then for the elements g;; of G: (i) g;; , for 1 <7 <a, 
1 < ;} <@,isa multiple of the greatest common factor of d;; and f;; ; (ii) g;; for 1 <7 < a, 
B <j <n, is a multiple of d;, ; (ili) g;; , fora <i < n, 1 <j < B, isa multiple of 
fi: Cv) 9i;,fora <i<n,8 <j <2, is identically zero. Thus g,; cannot be different 
from zero when both d,; , f;; are identically zero, and in each of the four cases above, 


(6) has a solution. 
Let U = PX Q", then U” = PX"Q™' (PXQ"')”, so Eq. (5) may be written as 


(PAP™")(PX"Q"') + (PXQ"')(QBQ"') = PCQ". (7) 


Multiplying (8) on the left by P™' and on the right by Q, we have Eq. (1); thus X 


P™ UQ isa solution of Eq. (1 


FURTHER PROPERTIES OF CERTAIN CLASSES OF 
TRANSFER FUNCTIONS: II** 


By A. H. ZEMANIAN (College of Engineering, New York Unive rsity ) 


This note is a sequel to a previously published paper [1]. The notation and termi- 
nology used here is the same as before. The purpose of this note is to point out a conse- 
quence of some previously published results [1, 2], which is immediately applicable to 
rational transfer functions that have no poles in the right-half plane and have at least 
twice as many poles as zeros. Such transfer functions arise quite commonly in physical 


systems. 





*W. E. Roth, The Equation AX — YB = C, and AX — XB = C in Matrices, Proc. Am. Math. 
Soc. 3, 392-396 (1952). 
**Received October 5, 1960. 
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Let F(s) be such a transfer function. It may be written as follows, where m > 1. 


F(s) = Ge + Gns8" $F +++ + ao _ NO), 
8b bg" + ee) +b = DGB) 





(1) 


The symbol f(t) will denote the inverse Laplace transform of F(s) that is zero for negative 
t. (This is the unit impulse response of the system.) 


Theorem: Let F(s) be given by (1) and let D(s) have no roots in the right-half s-plane. 
For m odd, letn < (m — 1)/2 and, for m even, letn < m/2. Then 


SE waliodie Galt la, | ¢"~? lait 
: +--+ +7 + . (2) 


()|< - am€ "l= 
j (m—-n — Dit (m — n)! (m — 2)! (m — 1)! 


Proof. Let D(s) be a Hurwitz polynomial and let m > 2, at first. By Theorem 3 of 
[2], 1/D(s) isa subclass m function. Hence, by Theorem 2 of [1], s*/D(s) is a class (m — yp) 
function where yu is any integer in the range, 0 < uw < (m — 1)/2, if m is odd and in the 
range, 0 < » < m/2, if m is even. Under these restrictions on u, the inverse Laplace 
transform of s*/D(s) is the wth derivative of the inverse Laplace transform of 1/D(s) 
Letting £* denote the inverse Laplace transform, Theorem 3 of [1] implies that 


| a ao 
D(s) —(m—n-—1)! 


Combining this result with (1), (2) is obtained. 
If m |, the hypothesis implies that m = 1 and n = 0. In this case, (2) becomes 
f(t) | < | a | and this result is easily established. 

The proof is completed by noting that the conclusion still holds when the restriction 
that D(s) is a Hurwitz polynomial is replaced by the weaker restriction that D(s) has 
no roots in the right-half s-plane. This is so since the unit impulse response is a con- 
tinuous function of the pole positions of the corresponding transfer function. 

The principal advantage of this theorem over the previously published results is 
the following. In the previous cases one always had to establish whether a transfer func- 
tion is a member of some class of functions before a particular bound could be placed 
on the corresponding unit impulse response. In this case one must establish that the 
degrees of N(s) and D(s) are appropriate and that D(s) has no poles in the right-half 
s-plane. The restriction on m and n can be checked by inspection while in many practical 
cases it. is known beforehand that D(s) has no roots with positive, real parts. For instance, 
if the system under consideration is a lumped, linear, fixed, finite, and passive one, the 
transfer function will be rational and have no poles in the right-half s-plane. 
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ON THE DEFINITION OF STRESS RATE’ 
By E. F. MASUR (The University of Michigan) 


The concept of stress rate has been the subject of several current papers. Prager [1] 
has introduced a discussion of some of the proposed definitions from a Cartesian spatial 
point of view and, in reviewing their suitability for use in the constitutive equations 
of plasticity, has arrived at a decided preference for Jaumann’s [2] definition. Naghdi 
and Wainwright [3] have generalized the concept of tensor rate in such a way as to 
include the definitions mentioned by Prager as special cases. Convected coordinates 
have been employed by Sedov [4] in a discussion of Prager’s results. In continuing the 
approach adopted by Sedov we investigate in the present note the connection between 
the various stress rates and the concept of a convected derivative introduced by 
Oldroyd [5]. Additional arguments in favor of Jaumann’s definition are established; more- 
over, a modification of the latter is proposed in connection with the constitutive equa- 
tions of plasticity. 

In order to make this note reasonably self-contained, a brief resumé of Oldroyd’s 
derivation is given in what follows. Let x; be the instantaneous Cartesian coordinates 
of a particle * and let a* represent a convected coordinate system defined by x; = 2,(a", t) 
in such a way that the material coordinates a“ of a particle are time independent. Let 
o°*, the contravariant components of stress relative to the material coordinate system 
in its instantaneous position, be considered functions of the material coordinates a* 
and of the time ¢t. The partial derivative of o** with respect to ¢ is called “‘the convected 
derivative’ by Oldroyd. The components D,c“” so defined form the contravariant com- 
ponents of a tensor, since the functional relationship between the assumed convected 
coordinate system and some other coordinate system which is also convected is inde- 
pendent of time. 

We are therefore justified in defining a tensor o’”’ by 


ol a D,o (1) 


and its components relative to the assumed Cartesian spatial system*® by 


enh eae (2) 

Since 
Ci; = 2,085.90" (3) 

and 
Dw, = 0; ., (4) 


where v,; represents the Cartesian components of the velocity vector, we obtain by 
differentiation of Eq. (3) Oldroyd's definition 


ed — D,0;; ine: Vi kT jk ee Vj KO ik 7 (5) 


‘Received November 21, 1960; revised manuscript received December 27, 1960. This note is part 
of a research program sponsored by the National Science Foundation. 

An extension to non-Cartesian spatial systems is straightforward, but contributes nothing essential 
to the purpose of this discussion. 

*Greek indices refer to the convected coordinate system and Latin indices to the fixed Cartesian 
system. 
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Similarly it is possible to define a tensor o’’ whose covariant components are given 
by D,o., . Note that the tensor so obtained is in general not the same as the one defined 
in Eq. (1) (that is, o/5and g.,gs,0””* are in general different), because the metric tensor 
Jag iS time dependent, hence convected time differentiation and lowering and raising 
of indices are not commutable operations’. If the Cartesian components of o”’ are ob- 
tained by a process similar to the one employed above, we are led to the definitions of 
Cotter-Rivlin [8]. Furthermore, if the same process is applied to the Kirchhoff tensor, 
then the use of the contravariant components leads to Truesdell’s [7] definition, while 
an analogous derivation can be carried out without difficulty in connection with the 
covariant components. 

The four stress rate definitions introduced thus far are rejected by Prager on the 
basis that their vanishing does not make the stress invariants stationary; consequently 
their use in the constitutive equations of plasticity leads to complications. This be- 
comes immediately apparent in a convected coordinate system since the use of either 
contravariant or covariant stress components necessitates the introduction of the time 
dependent metric tensor in the invariants; the latter are stationary in the presence of 
stationary stress components only in the case of locally rigid body motion. 

This argument for rejection can be strengthened further by the (not unrelated) 
consideration that the establishment of the constitutive equations of plasticity may 
make it necessary to separate the stress rate tensor into its isotropic and deviatoric 
contributions. Moreover, it is reasonable to demand that the vanishing of the stress 
rates make the deviatoric stress components stationary. 

This is not true for any of the definitions introduced thus far. For example, if the 
contravariant components of the stress deviator, which are defined by 


ap as ny 8 4 
Ss =g — ig” Iyg° (6) 


are used for the definition of the stress deviator rate” 


gr? = Ds? (7) 


, 


then it follows that «’”’ = 0 does not necessarily imply s’’”’ = 0. Conversely, with the 


introduction of the strain rate tensor e by means of 


Zeus = Va.p + Usa = D.Gas (8) 
the vanishing of s’’”’ implies 
rnb ‘ 1/8 d <8 
7 aor” 62 — 2a2(ef — 46 65). 


It is seen that for this case o’”’ is not isotropic, as it should be, unless « is isotropic too. 
The objections raised against the stress rate definitions introduced above become 


‘It is equally valid to refer to o’’”’ and to o”’ (or to any suitably normalized combination thereof) 
as a “‘convected derivative’’ of o; any rate tensor so defined clearly satisfies the requirement of being 
insensitive to rigid body motion. The conceptual significance of this ambiguity has apparently not been 
fully explored by Oldroyd [5], [6] or Truesdell [7], but is implicit in [4]. 

‘Note that this defines the rate of the stress deviator, which is not equal to the deviator of the 
stress rate. 
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invalid with the use of the mixed components. We therefore define the unsymmetric 


tensors *c’ and **o’ by means of 


ol? = Dio," (9) 
+e" = Do’. (10) 

and the symmetric tensor® 
ao” = 4(*o’ + **o’). (11) 


The Cartesian components of this tensor are given by 
of, = Diya; — pio jp — We; Oi (12) 
in which 

= 4(v;,; — 0;,;) (13) 


are the components of the skew-symmetric vorticity tensor. 


It is seen that o’ is Jaumann’s definition; in other words, the artifice of introducing 
a local rigid coordinate system moving with the particle as employed by Jaumann and 
again recently by Oldroyd [6] is equivalent to the use of the symmetrized convected 
time derivative of the mixed components. Note also that o’ is the mean of o”’ and o’” 
An analogous definition applied to vectors leads to a rate definition which is in agree- 
ment with elementary concepts. 

Under certain circumstances, the use of a weighted tensor in the constitutive equa- 
tions may be more convenient. For example, Drucker’s argument [9] establishing the 
convexity of the plastic loading function may lead to the conclusion that this function 
should be expressed in terms of the Kirchhoff tensor components, since the latter govern 
the rate of work per fixed unit of mass (rather than that per time dependent unit of 
volume). The equivalent Cartesian stress rate, referred to the instantaneous system, 
is given by 
Oj = Dio, — weit je — Wei HVE RC;; (14) 


which is acceptable in the light of the preceding discussion. 

It may be of some interest to discuss briefly the question of the definition of the 
strain rate. Offhand, it appears as if similar arguments should make it necessary to 
define the latter also in terms of the mixed components. This is actually not the case 
since the definition of the strain tensor 7,, itself employs the metric tensor and, con- 
sequently, some of the resulting relationships become somewhat simpler. In fact, with 


the strain deviator defined by 
has = —} ” (15 
lap = Nap 3779 Yas 2) 
then its rate, in view of Eq. (8), takes the form 


Cag = Dihas = (€as — 4€27gas) — 3(nr€a3 — 1” €yrGas)- (16) 


‘Of course the symmetry of o’ is a consequence of that of o. Similarly Eq. (11) defines an anti- 
symmetric tensor if the original tensor is also antisymmetric. In particular, if the operation is performed 
on w itself, then (as is plausible from elementary considerations) the material time derivative coincides 


with the convected derivative. This fact is useful in the determination of higher order rates. 
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In this case the vanishing of the strain rate components €,, implies also that the deviator 
components h,, are stationary. Moreover, a limited inverse is also true; that is, the 
vanishing of e,, implies that the strain rate is isotropic. 
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STEADY SPHEROIDAL VORTICES—MORE EXACT SOLUTIONS TO THE 
NAVIER-STOKES EQUATION* 


By VIVIAN O'BRIEN (Applied Physics Laboratory, The Johns Hopkins University) 


Abstract. The vorticity equation, the curl of the Navier-Stokes equation, is con- 
sidered in ellipsoidal coordinates. The steady spheroidal vortex solutions are demon- 
strated as examples of a class of exact flow solutions characterized by a simple linear 
vorticity distribution. 

Introduction. The Navier-Stokes equation expresses the conservation of momentum 
in a form suitable for continuum fluid mechanics. For incompressible fluid flow the 
system of the continuity equation plus the Navier-Stokes equation has a solution that 
is completely determined, in principle. However, solutions that satisfy the equations 
exactly are rare due to the non-linear convective term. One such solution that should 
be included in the meager list of exact solutions is the Hill-Hadamard spherical vortex. 

The spherical vortex is an example of a class of steady axi-symmetric flow solutions 
consisting of a rotational part where the vorticity is proportional to the distance from 
the axis and an irrotational part. Other examples of this class will be given, in particular, 
the steady oblate spheroidal vortices which appear to be physically significant. 

Hill-Hadamard spherical vortex. Hill {1] showed in 1894 that the stream function 

y = sin’ 0(Ar’ + Br’) (1) 


*Received December 5, 1960. Supported by the Bureau of Naval Weapons, Department of the 
Navy, under Contract NOrd 7386. 
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satisfied the inviscid condition for steady motion (Lamb [2], p. 245) 


ay 3° l dy 
A es ee 
2 eke) DOD 


= Dy = of(y), (2) 


where (x, @, ¢) is the cylindrical polar coordinate system. The relation to the spherical 
polar coordinate system (r, 6, ¢) is given by x = r cos 6,@ = rsin @. Hill’s solution takes 
f(w) as a constant = —10B. 

Hadamard [3] (and Rybezynsky [4]) in 1911, while discussing the viscous motion 
of a fluid sphere within another fluid, used the same stream function (1). Stokes [5] 


had given 


y = sin’ (Ar + Br‘ + Cr + Dr) (3 
as a solution to the steady viscous vorticity equation 

fr oa 10? 

\<5 +5 — = Spy = Dy =0. (4) 

Ox” T kom wiul v v 


The physical requirement of finite velocities as r > 0 eliminated the two lower powers of r. 

Thus it has been known for a long time that this stream function (1) independently 
satisfies a condition derived from the Navier-Stokes equation under the assumption 
of purely inviscid flow (2) and the condition derived for purely viscous flow (4). It 
follows, though never previously shown to the author’s knowledge, that the same 
stream function satisfies the complete Navier-Stokes equation. Therefore, it is an exact 


solution. 
The velocity field follows from the derivatives of y. q = ui, + vi; where 
, 1 dy 1 ay 
= =—— — (i= — 
ada’ D Ox 


The pressure field is found by integration of the momentum equation written in the 
“vorticity form”’ 


ra} f ] 
4 a xe=-V(2+ 04404) 1 Xe (5) 
with the vorticity vector o = V X q = wi, = (1/aDy)i, . Whence 


Il 


P4o4-@¢@ = 10B{Ao + Boz? + ow) — Que} 
p 


Il 


10B{(Ar’? + Br’) sin’ @ — 2vr cos 6}. 
For a spherical vortex y = 0 at r = a. This requires A = — Ba’. The Hill-Hadamard 
spherical vortex is one of the class of flow solutions consisting of a rotational part where 
the vorticity w is proportional to m and an irrotational part. Spherical vortices have 
been demonstrated in real fluids by Spells [6] and Savic [7]. 

Poiseuille flow. A more familiar example of this class of exact flow solutions with 
a linear vorticity distribution is the well-known Poiseuille flow. This is often quoted as 
an exact solution, but the vorticity distribution receives little attention. Its stream 
function can be written (a = pipe radius) 

y= Av(a -@). (6) 

This can be rewritten into a rotational part and an irrotational part as 


( 6, oil ers 
y =< 64 r* sin” or + {atc ”? (cos 6) + 203" * (cos 0}. 


(vorticity w = 12 Ar sin 86) 
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The irrotational part is given in terms of the Gegenbauer polynomials Cy'’*(cos @) 
after Savic [7] rather than the more usual Legendre polynomial form for reasons that 
will become apparent in the next section. The validity of this solution for real fluids in 
laminar flow is well-established. 

Exact solution in oblate ellipsoidal coordinates. It will be demonstrated that the 
steady spheroidal vortex is another example of the same class of exact flow solutions. 
For this purpose the pertinent equations are first transformed to the oblate ellipsoidal 
coordinate system (n, 8, ¢) by means of the relations 


x = sinh 7 cos 8 


cosh 7 sin B cos @ (7) 


II 


y 


| 


z = cosh nsin sing (@ = cosh nsin 8) 
The continuity equation becomes (see Goldstein [8], p. 114) 
0 | . . 2 g\1/2,, 7 a) . . 12 2 9\1/2, 7 7) 
: [cosh 7 sin B(sinh » + cos B)’~u,] + ry: [cosh n sin B(sinh » + cos B)’~vs] = O (8) 
On dl 
for q ui. + Sake 
Define a current stream function y by 

(cosh nsin B)' dy 


u.= -— TG ; _—— 
(sinh* » + cos’ 8)'’* a8 


(9) 
(cosh nsin B) ' Oy 
(sinh* » + cos” B)'”* an 


vs = 


, 


to identically satisfy (8). The steady vorticity equation can then be written 
(see [8], p. 115) 


- (cosh nsin B)' fay, Dy) 2 Dy _ a(y, cosh 7 sin 8) | ~ ots, (10) 
sinh” 7 + cos 8) | (ny, B) cosh 7 sin 8 O(n, 8) 
where 
= ' 9° a 9° 4) 
D sinh” » + cos £) (2; — tanh " 5y + 5a — cot B 2). 

Let 

Y = y, (irrotational) + y, (rotational). (11) 

The irrotational part ¥, , where 
] 


represents a trivial solution to (10) but it is convenient to discuss this first. 


Assume ¥, = f(n)g(8) and Dy, = 0 is separable for n(n — 1) > 0 into the pair of 
ordinary differential equations 


{’’(n) — tanh nf’(n) — n(n — 1)f(n) = 0; g’’(8) — cot Bg’(8) + n(n — 1)g(8) = 0. 
By applying the transformation 7 sinh » = z to the former and cos 8 = z to the latter, 
both equations reduce to a Gegenbauer differential equation (see Magnus and Ober- 
hettinger [9], p. 76) 

(1 — 2°)P’(z) + n(n — 1)P(2) = 0. (12) 


A solution for this is P = C,'*(z), the Gegenbauer polynomial of order n 


and degree (— 3 


> 


). 
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Thus a general solution y, valid for regions containing the origin is 
y= 2. AC,” (eos B)Cz'”? (i sinh 7). (13) 


The n = 0, 1 terms have been dropped because they lead to velocity singularities 
along the axis. The Gegenbauer functions of the second kind which also satisfy Eq. (12) 
lead to infinite velocities as 8 — m/2 or 7 — O and thus do not appear as part of the 
physical flow solution. However, the Gegenbauer functions of the second kind with 
arguments a function of » would be important for other flow problems than the type 
considered here. 

Consider now a rotational part 

¥. = Bsin’ B cosh’ n (cos B — sinh 7). (14) 
This makes both the inviscid part of the vorticity equation and the viscous part identi- 
cally zero. w2 = Dy, can be transformed by cos 8 = uy, 7 sinh n = z into 
‘ile — , 7 2 ' er . 
(u —z) I —w)=3-—- (0 — 2) = JB — pw) — 2) t+). 
Ou Oz 
Carrying out the indicated operation, this equals 
—10B(1 — w)(1 — 2’) —10B sin’ B cosh’ 7 


= KW. 


I 


(Or w@ = KW.) 
The inviscid part of the vorticity equation (10) gives 


| OY. Ir OY. Aye | | 

. | nea ea e ae. | 

(cosh n sin 8)" | On 0B i On 0B 
eS ss 3K — 2ko \ 
(sinh” 7 + cos 8) | a , ‘ | 
0D dD 0D OD | | 

On Of | | On 0B |j 


which is identically zero. The viscous part vD°y, is also satisfied because the term equiva- 
lent to Dy, appears in the general solution y, for n = 2 or Dy, = 0. 

The oblate spheroidal vortex. The exact solution y = y, + ¥. with the boundary 
conditions y = 0 on the ellipsoid 7 = 7 and that there be stagnation points on the 


central axis yields: 
‘ - : 
18. 2 2 . 92 ee cos B — et ae , 
y = Csin’ 8 cosh n\sinh n — sinh 7 + ~>>s— (sinh 1 — sinh 7)?. (15) 
sinh” mo J 


This is a spheroidal vortex. The velocity follows directly from the derivatives of y and 
the pressure field can be obtained by integration of the momentum equation in the 
vorticity form. The component momentum equations in this coordinate system are 


Jogi , 172 O [Dp , ° 
—vgo = —(sinh’ y + cos Bp)” ; (: + 2+ 3[u, + il) 
On \p 
eo r 1/2 Ow 
+ y (sinh 7 + cos 8) w cot B + =], 
op 
° 2 2 1/2 a) Pp 1 2 2 
u,w = —(sinh’ 7 + cos 8) 98 + 2+ 3[u, + v5] 
08 \p 


° > 2 > ra] 
—yv (sinh? » + cos’ B)~’ *(« tanh n+  ). 
on 
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The pressure field is 
P + 2+ 3(u; + v5) = xf —w + 2 cos Bsinh 7}, 
p 


where x = —10C is the proportionality factor between the vorticity w and the axial 
distance @. 
Transforming back to the cylindrical coordinate system 
Y = Cw (cosh’ » — w — 2° coth’ m), (16) 
where y = 0 on the oblate ellipsoid 


z* (sinh 99)"' + @° (cosh’ 4)7' = 1. 


An example is shown in I'ig. 1 for an ellipsoid with fineness ratio = tanh np) =.5(m)=0.55). 
The velocity components parallel and normal to the axis are 


u = —2C (cosh’ "o — 2a — 2x coth No); 
v = —2C(ro coth’ no). 
In the axial plane, x 0, note a circle of zero velocity. This center of circulation occurs 


at the same fixed ratio of the semi-major radius, (2)~'’*, for all values of m . 


tanh 9, = 0.5 





stream- 
lines 
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Fig. 1 
An example of a spheroidal vortex 
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Real vortices often have boundaries that can be approximated by oblate spheroids. 
A comparison of some real flows with theory is given in another place, O’Brien [10]. 
The measured centers of circulation of quasi-steady vortices over a fairly wide Reynolds 
number range compare reasonably well with the theoretical value given here. 

A more general viscous solution. Neglecting the inviscid part of the steady vorticity 
equation, a rather general solution, ¥. , can be determined from the relation Dy, = y, . 
Although the equation is not separable, the rotational solution in oblate ellipsoidal 
coordinates (valid for the region containing 7 = 0) has been found to be 


v2 = > B,C,” (cos 8)C;'*(i sinh )(cos’ 8 — sinh? 7). (17) 


But for n > 3 this stream function does not simultaneously satisfy the inviscid part. 
The prolate spheroidal vortex. In similar fashion the problem can be reformulated 
in prolate ellipsoidal coordinates. The stream function for a prolate spheroidal vortex is 


yY = Co (sinh’ & — ow — 2’ tanh’ é,). (18) 
In analogy to the oblate cases, the center of circulation is at the same fixed ratio, (2)~'’*, 
of the midsection radius for all & . 

Conclusion. It is not known whether there are other members of this class of flow 
solutions. The second order Laplacian equation that governs“ perfect” fluid flow is separa- 
ble in just eleven coordinate systems. The fourth order purely viscous vorticity equation 
(4) is not even separable in the two coordinate systems considered here (oblate ellipsoidal 
and prolate ellipsoidal). Assume, by use of another coordinate system, a new solution 
to this linear equation is found. If it has vorticity proportional to the axial distance, 


@, it will be another exact solution. 
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A NOTE ON UNIQUENESS IN CLASSICAL ELASTODYNAMICS* 
By M. E. GURTIN AND ELI STERNBERG (Brown University) 


Introduction. The conventional uniqueness theorems of classical elasticity theory, 
first established by Kirchhoff [1] for the equilibrium case and later extended to elasto- 
dynamics by Neumann [2], rest on the assumption that the elastic constants satisfy 
inequalities which are necessary and sufficient for the positive definiteness of the strain- 
energy density. The theorems referred to thus presuppose that 


up > 0, —-l<o<}, (1) 


if uw and o@ respectively denote the shear modulus and Poisson’s ratio of the material. 

Actual elastic materials obeying the linear isotropic stress-strain law conform to (1) 
and, in fact, fail to exhibit negative Poisson ratios. Nevertheless, the question as to 
whether conditions (1) are necessary for uniqueness, is of obvious theoretical interest. 
It is the purpose of the present note to show, with limitation to the first boundary- 
value problem of elastodynamics (surface displacements prescribed) and for bounded 
domains, that the inequalities (1) may be relaxed without loss of uniqueness. Specifi- 
cally, we find that in this instance (1) may be replaced by the less stringent assumption 


uw > 0, —2 <a < }, L<<e«ct @, (2) 
or, equivalently, by 


gat>e, gu Are, (3) 
oo (1 — 2c)p — 
in which ¢, and c, are the velocities of equivoluminal and irrotational waves in a medium 
of unlimited extent, whereas p stands for the positive mass density.” Accordingly, the 
generalized inequalities (2) admit a simple intrinsic interpretation: they require that 
the two basic wave-velocities of the elastic solid be real. 

The foregoing extension of Neumann's theorem, which will be proved subsequently, 
is the dynamic counterpart of a known extension of the uniqueness theorem in elasto- 
statics. The solution to the first boundary-value problem in the equilibrium theory is 
unique (under suitable regularity assumptions) if, instead of (1), one stipulates merely 
that 


un ~ 0, —-~e~ <o<}, l<ao<om, (4) 


The relaxed inequalities (4) are evidently met provided the two wave velocities c, and c, 
are real and non-zero. For a bounded domain, this generalization of Kirchhoff’s 
uniqueness theorem is implicit already in observations due to Kelvin [3]; it may also 
be inferred from available results appropriate to the Dirichlet problem for strongly 


*Received December 15, 1960. The results communicated in this paper were obtained in the course 
of an investigation conducted under Contract Nonr-562(25) of Brown University with the Office of 
Naval Research in W ashington, D. C. 

Numbers in brackets refer to the list of publications at the end of this paper. 

*Actually, the inequalities (3) are sufficient for uniqueness in the present circumstances even if 
p is negative. 
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elliptic systems of partial differential equations,* since (4) hold if and only if the dis- 
placement equations of equilibrium are of the strongly elliptic type. As shown by Duffin 
and Noll [6],* the inequalities (4) assure the uniqueness of the solution to the first prob- 
lem of elastostatics also for unbounded exterior domains, provided the displacements 
are required to vanish uniformly at infinity. Finally, it is clear from counter-examples 
due to Ericksen and Toupin [8] and Ericksen [9] that (4), in the equilibrium case, cannot 
be weakened further without loss of uniqueness. 

A generalization of the uniqueness theorem in elastodynamics. The displacement 
equations of motion in the linear theory of homogeneous and isotropic elastic solids 
take the form 

du 


uV*u + = VV uti= pa, ’ 
oer = ot 


es 
— 


in which u(z, ¢) is the displacement vector’ and f(z, ¢) is the body-force intensity vector, 
the single argument x representing the triplet of cartesian coordinates and ¢ denoting 
the time. Equation (5), which must hold throughout the region of space D occupied 
by the medium, is subject to the initial conditions 


du 


. =vi(z) i (6 
At |e o, v(x) in D, )) 


u(r, 0) = u(z) 
where u(x) and v°(z) are the prescribed initial distributions of displacement and velocity. 
In the first problem of elastodynamics the accompanying boundary conditions are 
characterized by 

u(z, ?) = u*(z, 4) on B O<t<o), (7) 


if B is the boundary of D, while u*(z, ¢) are the given surface displacements. 

We now prove the following theorem. Let D be a bounded regular region of space® 
with the boundary B and suppose that (3) hold. Then there exists at most one u(x, t), twice 
continuously differentiable for (x, , 22,23) in D+ B(O <t < @), which satisfies (5) in 
D(O < t < &), and meets (6), (7). 

In view of the linearity of (5), (6), (7) it is evidently sufficient to show that 


2 aut i 
uVut - c — VV-u= p i in D (O<t< &), (8) 

1 — 20 ol 

together with 
0 . 

u(z,0) = — =0 in D, (9) 

OF | ¢z.0) 
u(z,/) = 0 on B (0 < t < @&), (10) 


imply u(z, 1) = 0in DO <t < @), 


Let K(t) be the kinetic energy of the entire body, so that 


K(t) = § [ || av, (11) 


3See Browder [4] and Morrey [5]. 

4See also [7]. 

5Letters in boldface designate vectors. 

*The term ‘‘regular region of space’’ is used in the sense of Kellogg [10]. 
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and define G(t) through 
. : mes x1 —o 
G(t) = ° | k(V-u)y + (V Xu) dV, c= a =). (12) 
/D 


2 1 — 2¢ 


From the divergence theorem, with the aid of (8), (11), (12), and elementary vector- 


analytic identities, follows 


»| E wVut x (7 x w | ndA = on + &, (13) 
where n is the outward unit normal of B. 
According to (13) and (10), 
K(t) + Git) =c O<t< ~), (14) 
in which c is a constant. But by (11), (12), (9), 
K(0O) = GO) = 0, (15) 
whence, from (14) and continuity, 
Ki) +G=0 O<t<.). (16) 


Since «x > 0 by virtue of (3), both K(8 and G(é) are either non-negative or non-positive, 
so that (16) yields 


K(t) = 0, G(t) = 0 O<t< oo). (17) 
On the other hand, the first of (17) in conjunction with (11) leads to 


= =0 in D (<t<.2), (18) 
Cc 
from which, because of (9), u(x, t) = u(z, 0) = 0 in DO < t < @). This completes 


the proot. 


REFERENCES 


1. G. Kirchhoff, Ueber das Gleichgewicht und die Bewegung eines unendlich diinnen elastischen Stabes, 
J. fiir die reine und angewandte Mathematik 56, 285 (1859) 

2. F. Neumann, Vorlesungen tiber die Theorie der Elasticitaét der festen Kérper und des Lichtathers, 
B. G. Teubner, Leipzig, 1885 

3. W. Thomson, On the reflexion and refraction of light, Phil. Mag. 26, 414 (1888) 
F. E. Browder, Strongly elliptic systems of differential equations, Ann. Math. Studies 33, 15 (1954) 

5. C. B. Morrey, Jr., Second order elliptic systems of differential equations, Ann. Math. Studies 33, 101 

1954 

6. R. J. Duffin and W. Noll, On exterior boundary value problems in elasticity., Arch. Ratl. Mech. and 
Anal. 2, 2, 191 (1958 

7. M. E. Gurtin and E. Sternberg, On the first boundary value problem of linear elastostatics, Arch. Rath. 
Mech. Anal. 6, 3, 177 (1960) 

8. J. L. Ericksen and R. A. Toupin, Implications of Hadamard’s condition for elastic stability with 
respect to uniqueness theorems, Can. J. Math. 8, 432 (1956) 

9. J. L. Ericksen, On the Dirichlet problem for linear differential equations, Proc. Am. Math. Soc. 8, 
3, 521 (1957) 


10. O. D. Kellogg, Foundations of potential theory, J. Springer, Berlin, 1929 











172 BOOK REVIEWS [Vol. XIX, No. 2 


BOOK REVIEWS 


(Continued from p. 142 


to be informed about the present state of the art, and the practicing engineer who has to solve a given 
concrete problem. The latter will find a wealth of useful material in the form of tables, graphs, formulas 
and references. To the first the limitation of the general chapters to linearized theory will probably be 
somewhat disappointing. For both groups, however, the book will prove to be an indispensable tool. 
H. ParkKus 


Stochastic population models in ecology and epidemiology. By M.S. Bartlett. John Wiley & 
Sons, Inc., New York, and Methuen & Co. Ltd., London, 1960. x + 90 pp. $2.00. 


This monograph discusses models for biological populations in ecological and epidemiological 
contexts where a full stochastic formulation is advisable. It is noteworthy that even such simplified 
models as are presented imply complex consequences. First various models to determine the frequency 
distributions for single species are presented; the Poisson model is the simplest of these; birth-and- 
death processes and the problem of extinction are next discussed. Growth and interaction, and com- 
petition, are complications leading to more complicated models. In the concluding chapters, epidemiologi- 
cal models are considered for the vicissitudes of animal or human populations subject to infection by 
some disease transmitted from individual to individual. The restriction to a closed population is relaxed 
to permit recurrent epidemics and in the last chapter, epidemiological problems are noticed in which 
the spatial or topographical factor is relevant. 

The work will appeal to biometricians with the necessary background, to theoretically minded 
ecologists and other interested research workers in biology and medicine. 

W. FREIBERGER 


Automatic language translation. By Anthony G. Oettinger, with a Foreword by Joshua 
Whatmough. Harvard University Press, Cambridge, Mass., 1960. xix + 380 
pp. $10.00. 


About a dozen projects in mechanical language translation are being currently conducted in the 
United States alone; the great majority of them have, as their ultimate aim, a reasonably accurate 
and smooth translation from Russian into English, mostly in scientific and technical fields. This is also 
the case with the research being conducted at the Computation Laboratory at Harvard University of 
which Mr. Oettinger is the principal investigator and this book a thorough and informative progress 
report. 

The book under review here deals with only one of the major problems which have to be resolved 
in machine translation research: the compilation of an adequate bilingual dictionary in a form which 
would make it readily accessible to the computer (stored on magnetic tape in the Harvard project), 
and the devising of programming procedures through which the words of a running text could be 
accurately matched with entries in this dictionary, thus achieving one of the essential steps of transla- 
tion, the finding of lexically equivalent or near-equivalent elements in two languages. Once this step has 
been mastered, a rough approximation of the original text in another language, a so-called word-by-word 
translation, can be achieved. Word-by-word translations are of some use but are basically unsatisfactory 
because they contain no information about the grammatical properties of words or their syntactic 
relationship (except what may be communicated through lexical characteristics or the preservation 
of the original order of items), not even speaking of their complete disregard of such major problems 
as multiple meanings of words, idiomatic and other non-equivalent constructions, correct word-order, 
etc. Nevertheless, the compilation and operation of an automatic bilingual dictionary is a prerequisite 
to research in the other essential aspects of translation by data-processing machines, be it the mechanical 
analysis of the grammar and syntax of the source langugage, or the generating of grammatically and 
syntactically acceptable strings in the language into which the translating is being done. Oettinger 
discusses these other aspects only in general terms, pointing out that the Harvard project has not 
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progressed beyond the initial stages in such research. However, the fact that this book deals primarily 
with the compilation and operation of an automatic dictionary, possibly one of the more easily man- 
ageable of the major steps in mechanical translation, does not detract from its importance. In spite of 
the occasional optimistic newspaper reports, the time when a good translation can be obtained entirely 
by automatic means is not exactly around the corner. Thus, the thorough treatment which Oettinger 
gives to the work already accomplished at Harvard and his sober discussion of the limitations and 
difficulties which still must be overcome in the future will be appreciated by those researchers interested 
in automatic linguistic analysis and translation, who are aware of all the complexities of this field of 
research. 

The book is divided into ten chapters, which are interspersed with numerous tables, diagrams, 
appendices, and reproductions of actual machine-produced lists. Some chapters also contain special 
bibliographical information in addition to the sizeable bibliography at the end of the volume. Although 
the author tried to make the discussion of information-processing machines and their programming 
as general as possible, the book is unavoidably strongly oriented toward the Univac I computer which 
the Harvard project has been using, and some of the special characteristics as well as the limitations 
of this machine have determined the approach to individual problems. 

The first chapter is a general introduction to the organization and programming of data-processing 
machines, with the Univae system serving to illustrate the individual points. The chapter begins with 
a very elementary outline of the principles of such machines, obviously intended for the complete 
beginner in the field, but the speed with which the subsequent complex information is introduced is 
such that one may doubt whether a reader who has never studied a computer manual would be able 
to follow the explanation for more than a few pages. The first chapter also contains several sample 
problems and a number of exercises (without solutions supplied). The exercises, in the later chapters, 
get progressively rarer until, toward the end of the book, their occasional appearance in the midst of 
the discussion cannot but impress the reader as a distinct afterthought. I found this apparent hesitation 
about the format of the book, and the sporadic and not very effective attempts to give it a bit of a 
textbook character, rather distracting. 

The second chapter, which deals with the structure of signs and pays particular attention to the 
properties of sets, is well organized and informative but still comprehensible to the non-specialist. 
It also contains a description of the Cyrillic Unityper, a specially adapted typewriter used for transcribing 
characters onto magnetic tape. The machine serves, in lower-case, for encoding Russian characters 
while, in upper-case, the usual layout of Roman letters is maintained. Chapter three discusses flow 
charts and automatic programming, and includes an illustrative flow chart and programming information 
for the transformation of the computer representation of Russian characters into one of the commonly 
used transliteration systems of Russian through letters of the English alphabet. Chapter four investi- 
gates the various aspects of a translation process and states the assumptions which make mechanical 
translation a realistic, even if complex, possibility. The bulk of the book, chapters five to nine, are 
devoted to the compilation and operation of a dictionary. The Harvard Automatic Dictionary contained, 
at the time this book was written, some 22,000 stem entries (it has apparently been enlarged to 30,000 
since then). One stem entry does not necessarily correspond to a word in the conventional sense; one 
word may be represented by one or more stem entries, depending on the degree and complexity of the 
inflection which it undergoes. However, the Harvard Automatic Dictionary is not a “paradigm dic- 
tionary’ which would list all the inflected forms of a word separately, but rather something of a ‘‘canoni- 
cal dictionary’, containing only those stems which are necessary for the successful matching of the 
strings in actual text and, of course, a list of affixes, so that strings of the text can be factored before 
they are matched with the dictionary entries. The dictionary is compiled in several steps, some of 
which are automatic. The appropriate words are selected for inclusion from conventional dictionaries 
and other sources, and then manually classified into their various word-classes (nouns, verbs, ete.) and 
inflectional sub-classes. This classificatory information is included when the words are transcribed onto 
magnetic tape. The computer, supplied with a list of inflectional endings and a program for applying 
them, generates the complete paradigm of each word (if that word can indeed be inflected); after this 
has been accomplished, all the forms now in the list are automatically factored into stem and affix, 
duplications eliminated, and the final list of distinct stems entered in the dictionary, together with the 
English equivalents which have been assigned to them by a human translator. The running text which 
is to be analyzed is processed in a similar fashion, with words factored into stems and aff xes. Since the 
same algorithm is used in the generation of stems which become entry keys in the dictionary and in 
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the operation to factor strings of the analyzed text, there is a guarantee that a matching stem will be 
found for each word of the text, provided, of course, that such a word had been originally processed for 
inclusion in the dictionary. 

The final format of an entry in the automatic dictionary consists of 30 machine words (of 12 charac- 
ters per word) which include, besides the Russian stem and one or more of its English equivalents, the 
class designation of the Russian item and various other information, such as the specification of verbal 
aspect or syntactic government by verbs, which are intended for use in the subsequent automatic 
grammatical analysis. 

Mr. Oettinger also gives an actual example of a word-by-word translation (which he calls a “trot’’), 
containing all the possible English meanings included in the dictionary for every Russian word. Through 
manual post-editing which supplies some grammatical information and selects the appropriate English 
alternative, the “‘trot’’ can result in a useful, even if inelegant, translation. 

The last chapter, with the descriptive title ‘From dictionary to translation’’, discusses the future 
plans of the Harvard project which will hopefully lead to an acceptable mechanical translation of 
Russian texts. 

The book is well organized, with some minor exceptions. The numerous tables and diagrams, which 
are scattered throughout, are referred to in the text by their numbers rather than by pages on which 
they occur; this results in a tedious search when they are mentioned in various places in the book. 
The index which completes the volume leaves something to be desired. 

As a whole, this book is a competent and informative contribution to a complex field of research. 
As mentioned previously, it deals, in a detailed manner, only with one of the essential aspects of auto- 
matic translation, the lexical one, and it is careful to avoid any implications that the other major steps 
which still must be tackled are about to be resolved. After reading the book, one is left with the impres- 
sion that research in automatic translation has made considerable strides in the last few years but that, 
all the same, the day when one could have a secretary copy a Russian article onto magnetic tape and 
get, in a few minutes, a good English translation from the computer is not yet within sight. 

The foreword by Joshua Whatmough, which is promisingly announced on the title page and on the 
dust-jacket, is only a brief statement of somewhat whimsically expressed praise for the book. 

Henry Kucera 


A guide to mathematical tables. By A. V. Lebedev and R. M. Fedorova. English edition 
prepared from the Russian by D. G. Fry. Pergamon Press, New York, 1960. xlvi + 
586 pp. $15.00. 


The tables listed are arranged under the following headings: Powers, rational and algebraic func- 
tions—Trigonometric functions. Various values connected with the wide and the sphere—Exponential 
and hyperbolic functions—Common and natural logarithms—Factorials, Euler integrals and related 
functions—Sine and cosine integrals, exponential and logarithmic integrals and related functions— 
Probability integrals, probability distribution functions, and related functions—Elliptic integrals and 
functions—Legendre functions and polynomials—Cylinder functions—Certain special functions and 
integrals—Solutions of certain equations—Sums and values connected with finite differences—Mathe- 
matical constants—-Primes, factors and products, quotients and fractions. While the text has been 
translated all tabular matter has been directly reproduced from the Russian original. This involves 
minor deviations from the customary notations (e. g. the use of tg for tan) but has the great advantage 
of excluding typographical errors that might have arisen if this material had been reset. 


Stationary processes and prediction theory. By Harry Furstenberg. Princeton University 
Press, New Jersey, 1960. 283 pp. $5.00. 


The question of the “predictability” of a sequence is examined in this study. Specifically, given a 
sequence of observations --: , &—n), «++ , &—1), &(0) (possibly nonnumerical), the notion of prediction 
of future possible values is considered. The discussion is highly motivated by the statistical theory of 
prediction (see the work of A. N. Kolmogorov, N. Wiener and their followers) but is formulated in 
rather abstract language. The sequence of observations can be thought of as the realization of a random 
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process in the past and present. If the realization is consistent with stationarity of a potential random 
mechanism generating it, a prediction in terms of a probability distribution on possible futures is pro- 
posed. The presentation rests heavily on the theory of commutative C*-algebras. 

M. RosENBLATT 


The potential theory of unsteady supersonic flow. By J. W. Miles. Cambridge University 
Press, New York, 1959. xii + 220 pp. $8.50. 


This monograph is a purely theoretical one, essentially on the unsteady aerodynamics of slender 
bodies in supersonic flow. The treatment is limited to frictionless potential flows and, almost entirely, 
to linearized flows. Within these bounds of subject matter the author has prepared a comprehensive 
and well-organized monograph. In the book are evident not only the results of the author’s own extensive 
research in the field but also his thorough knowledge of the literature. 

The first three chapters introduce the subject, by presenting the basic linearized equations and 
discussing their transformations. Chapter 4, treating reduction of the equations to those for steady 
flow, includes the reductions of Magnaradze and of Garner. Succeeding chapters treat in detail two- 
dimensional problems, simple planforms, the rectangular wing, the quadrilateral wing, slender wing 
methods, the delta wing, the low aspect ratio rectangular wing, and non-planar bodies. A final chapter 
on nonlinear problems and an appendix on reverse flow theorems close the monograph. A bibliography 
and three indices are provided. 

This excellent monograph is one of a series of monographs on quite specialized topics in mechanics 
and applied mathematics. The availability in organized treatments of such material is of great value to 
those not already specialists in a field. Few indeed could easily get an understanding of the state of the 
art in the theory of unsteady supersonic flow from sources other than Miles’ treatment. The editors 
are to be commended. 

Watiace D. Hayes 


Digital computers and nuclear reactor calculations. By Ward C. Sangren. John Wiley & 
Sons, Inc., New York, London, 1960. xi + 208 pp. $8.50. 


This book talls naturally into two parts: the first four chapters constitute an introduction to digital 
computers, programming, numerical analysis and reactor problems, and the last four discuss in some 
detail a particular reactor problem, namely fission-product poisoning, as a typical example of the 
application of computers to reactor physics. 

The style is leisurely and most readable, the argument on the whole elementary. The book serves 
as an excellent survey for those without prior knowledge or specialised interest, or as an introduction 
with methods for generating random numbers obeying given probability distributions. In “Fourier 
to more advanced works, such as ‘‘Numerical methods for nuclear reactor calculations” by G. I. Marchuk 
(see Quart. Appl. Math., Vol. 18, p. 36), for anyone intending to work in the field. 

WALTER FREIBERGER 


Contributions to probability and statistics. Essays in honor of Harold Hotelling. Edited 
by I. Olkin, 8. G. Ghurye, W. Hoeffding, W. G. Madow, H. B. Mann. Stanford 
University Press, California, 1960. x + 517 pp. $6.50. 

This book is a collection of forty-two papers on probability and statistics, mostly mathematical 
statistics, dedicated to Harold Hotelling on his sixty-fifth birthday. Most of the papers are results of 


new research rather than expository papers. 
G. F, NEWELL 
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